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Introduction 


Cet ouvrage est destiné au public des étudiants de second cycle de sciences 
économiques et dé géstion, mais aussi plus largeruent à tous ceux qui souhaitent 
acquérir une formation de base en économétnio, Il s'agit donc d'un manuel dont 
l'ambition &st de présenter de façon claire et accessible au plus grand nombre, un 
ensemble de techmiques dont l'usage est devenu courant dans les milieux 
professionnels. C’est aussi un ouvrage qui privilégie les applications par rapport à 
la héoric purc. On x trouvera les techniques lés plus usuelles illustrées par des 
exemples d'applications, Mais l'économétrie est une discipline vivante dans 
laquelle existent une recherche active et un flux régulier de nouveaux 
développements. Rendre complé de tous ces travaux dans touté leur richesse 
dépasse largement le cadre de Cet ouvrage. ]l existe des revues. des manuels et des 
ouvrages qui leur sont consäcrés el auxquels nous renvoyons le lecteur qui 
souhaitera dépasser le contenu volontairement limité de ce manuel. Susciter, par 
cette première approche, le désir d'un approfandissement ultérieur où d'un effort 
pour suivre l'évolution de la discipline fait partie intégrante des objectifs de cet 
ouvrage. 


L'essentiel de ses dévelappements est tiré d'un cours enseigné sous différentes 
variantes de contenu et de niveau aux universités d'Orléans, d'Abidjan, de 
Niames, de Cotonou, de Batare et dérnièrermçnt en matttise à l'université du 
Maitic, 


L'économétrie est une branche de la science Écemomique consistant à établir des 
lois ou à vérifier des hypochèses à partir de données chiffrées tirées de la réalité, 
Les étudiants de première année apprennent que pour différentes raisons, il ost 
impossible d'utiliser en scientes sociales la méthode expérimentale, comme ecla 
se pratique dans les sciences de la nature. Cette limitarion poserait des problèmes 
redoutables concernant la validité des développements de Ja théorie économique 
s'il n'existait un substitut de l'expérimentation consistant à en confronter les 
résultats aux donnéés tirées de l'observation de la céalité, Dans cette démarche 
essentielle, l’économétrie fournit les méthodes qui permettront de téster les 
hypothèses où de quantifier tes relations. Elle à par conséquent un caractère 
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instrumental et constitue le complément indispensable de l'analyse économique. 
Dé ce fait il est à peu près impossible de faire de li recherche en sciences 
économiques sans $e trouver devant la nécessité de lire ou de réaliser des travaux 
d'économétnie à un moment ou à un autre. C'est la raison pour laquelle, duns tous 
les pays, la lormalion dés éconvrristes suppose l'acquisition de ces techniques. 

Mais l'économéric pormet également de quantifier cs relations entre variahles 
économiques, ainsi par axcrmple de chiffrer la réponse de l'investissement en 
logement à une variation du taux d'intérêt. Pour cela, l'écanamiste part d'une 
analyse théorique des déterminants de l'investissement en logement, identifiant 
les variables explicatives et Le type de relation qu'elles sont susceptibles d'avoir 
avec la variable expliquée. Ces éléments analytiques vont être traduits en une 
équation mathématique mettant en relalion l'investissement cn Logement et sus 
principales variables explicalives. En donnant aux variables les valeurs 
enrogisirées dans lt réalité, l'économétne pormet à là fois de tester Ia validité du 
modèle et d'en éhitlrer Les paramètres. Elle rend passible li modélisation, c'est-à- 
dire la représentation de phénomènes économiques, nécessairement complexes, 
par un ensemibfe cohèrent de relations mathématiques quantifiées. À son tour là 
modélisation constienc gn enrichissernent important, sous Ta forme d'une aide à la 
décision, en rendant possible la simulation, l'optimisation ou la prévision. La 
siroulation consisté à calculer les conséquences d'un ensemble d'hypothèses 
réprésentécs par des valeurs dus variables explicatives, La prévision étend cet 
uxtrcice à des périudes futures. L'optimisation adopte I& chemin inverse et a pour 
but de déterminer les valeurs optimales de éertaines variables en fonction 
d'objectifs fixés. La plupart des grands centres de décision, qu'il s'agisse des 
banques, des agences gouvernementales, des institutions intemationales et des 
grandes entreprises utilisent régulièrement ces techniques. 1.cs constructeurs 
d'avion wilisont des mudèles éconamétriques pour analyser la demande actuelle 
et future, Les Ministères des Finances, de l'Économie où du Plan construisent des 
modèles maëro-écencrmiques pour la prévision et pour l'évaluation des politiques 
économiques. Les institutions non gouvernementales d'aide au développement 
construisent el aclualisent régulièrement des modèles des grandes régions du 
monde et des principales économies nationales. [i en résulte que pour travaillec 
dans cos institutions, il est indispensable de pouvoir manier les techniques 


Introduction 3 


économétriques souvent à haut mvcau, Îl en résulte aussi que pour pouvoir 
travailler avec ces institutions, il faut maîtriser ce langage commun. 

Il est vrai que depuis quelques années la mode ést au libéralisme et à des 
politiques économiques aussi peu aclives que possible. Les marchés sont supposés 
tout régler dès lurs qu'ils ne sont pas entravés. Dans ces conditions, on peut douter 
de l'utilité des modèles macro-économiques, d'autant plus qu'ils se sont révélés 
peu performants, pour né pas dire plus, dans la période récente. De fait, l’un des 
débouchés traditionnels de la modélisalion, qui avait largement contribué à son 
développement dans lé domaine de La macro-éconumie et exercé un effet d'appel 
sur l’écunométrie, se trouve actuellernent quelque peu délaissé. On peut cependant 
pensér que les idées évolueront, en fonction des résultats de l’expérience, vers une 
plus grande autonomie reconnue à la politique économique. De leur côté les 
téchrniques économétriques ne cessent de progresser et la qualité des modèles 
s'améliore, De nouvelles furmes structurelles comme les modèles à correction 
d'erreur ont permis simultanément d’allégér lès équations et d'améliorer leur 
qualité dans de nombreux cas, en particulier en macro-éconemie. L'utilisation des 
modèles de séries temporelles développés par Box et Jenkins ou leur association à 
des modèles structurels permettent, de leur côté, d'améliorer sensiblement Jes 
qualités en prévision. De nouvelles méthodes d'estimation permettent de tenir 
compie des relations pouvant exister entre Les chocs aléatoires subis par plusieurs 
variables etc. Tout ceci amène à penser que le recours à la modélisation et à 
l'économétrie devrait connaître de nouveaux développements dans l'avenir. Les 
applications à l’entreprise ét au marché en général, malgré les problèmes 
particuliers qu'ellés posent pourraient bien représenter la partie la plus active de 
ces futurs développernents. 

Par ses méthodes comme par la nature des informations qu'elle utilise, 
l'économétrie fail partie des disciplines quantitatives. On sait que depuis fort 
longtemps la question du rôle et de la place de ces techniques en sciences 
économiques fait l’objet de discussions acharnées. Si la possibilité d'utiliser la 
puissance de l'instrument mathématique en raison du caractèré quantifié des 
variables constitue un avantage indéniable, le caractère réducteur de la 
formalisation, la nécessité de postuler des comportements simplistes {les seuls que 
l’on puisse traiter facilement}, le recours systématique à l'abstraction sont autant 
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d'aspocts qui prêtent à la critique, Le simple fait que lors d’un des grands débats 
méthodologiques des années soixante on ait pu forger le concept dé “réalisme 
déplacé" {misplaced concretness}, en analyse économique. montre biën que pour 
les tenants les plus extrêmes de l’abstraction lorsque la réalité est en contradiction 
avec la théarie, c'est la réalité qui a tort. Il y aurait beaucoup à dire sur le rapport 
de la théorie économique à son objet, sur le rôle éclairant ou mystificateur dé ja 
formalisation. Force est dé constater que cétlie dernière constitue sousent La seule 
originalité de nombreux travaux et dissimule plus qu’elle ne supplée l'insuffisance 
théorique. Un exercice salutaire consiste à reprendre les conclusions et à Les 
confronter directement aux hypothèses de départ. On mesure micux, alors, ce 
gu’apportént les étapes intermédiaires du raisonnement. 


L'’économétric échappe largement à ces débats malgré son caractère formalisé. En 
effet, elle s'inscrit dans une démarche inductive qui s'efforce de découvrir à partir 
de phénomènes particuliers, les lois générales susccptibles de les avoir engendrés. 
À l'opposé d'une généralisation hâtive de modèles a prior, l'économérrie 
s'appuie sur la réalité qu'elle permet de micux déchiffrer, Elle se nourrit de cette 
confrontation pérmanente Entre les hypothèses théoriques à la réalité telte qu'elle 
peut être représentée par les données. Dès lors la recherche en économéirie, dont 
le caractère formalisé est évident, n'a d'autre but que de permettre une meilleure 
appréhension de la réalité. En perfectionnant les méthodes d'estimation, cn 
précisant leurs conditions de validité, elle offre la possibilité d'obtenir des 
résuitats plus précis ct plus fiables, Ainsi Les développements récents de la théorie 
de la coïntégration ont permis de micux cerner lés conditions de validité des 
procédures d'estimation courantes et d'éliminer les corrélations fortuites. 
L'utilisation généralisée des tests de racine unitaire qui en résulte augmente la 
fiabilité des résultats numériques. C'est pourquoi même si la matière de ces 
travaux est sauvent d'un accès difficile, leurs conclusions et les procédures qui 
peuvent on découler sont d'un grand intérêt. La formalisation n’est pas vaine, elle 
est au contraire le mcilleur chemin vers le concret. À partir d’un petit nombre de 
lois, comme les principales lois de probabilité. l'économétnit a construit un corps 
de savoirs essentiellement tourné vers tes applications. 


Comme dans d'autres disciplines, les problèmes qui se posent au niveau de ces 
applications renvoient, généralement, à la théorie en précisant par exemple la 
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distnibution de probabilité d’un type d'estimateur dans des conditions données. Là 
encure, la cointégration fournit un bon exemple puisqu'après avoir moniré que 
dans certaines conditions les distributions n'étaient plus standards. les 
économètres ont calculé les vraies disribetions ou leurs approximations. 


La connaissance des procédures, si lle est importante. n’est jamais suffisante. 11 
faut également maîtriser {eur conditions d'emploi et leurs limites. Le détour 
théorique, que l'on à voulu aussi simple que possible, n'est donc pas inutile. 


Chapitre I. L'estimation d’une équation 
unique 


Estimer une équation. c'est déterminer la valeur dus paramètres de la foncuon qui 
va relier Une variable à ses déterminants présumés. 

Verte t@ry +4 (l) 

Les paramètres à estimer sant @1.....45 

Si on était dans un univers délerministe, la loi en queslion serait parfaitement 

vérifiée, mais ve n'esl pas Le vas pour plusieurs raisons : 

— Les phénomènes économiques sont caructérisés par l'interdépendance entre de 
nombreux éléments, te qui entraîne que Les variables explicatives susceptibles 
d'exercer une influence sur la variublt expliquée sont très nombreuses {on ne 
Les retient pas taut£s). Mais. l'elfel des vanables qui ont été ornises explique 
qu'il y ait des écarts entre Îa réalité observés 1 le résultat de La fonction. 

— Les phénomènes économiques sont mis en œuvre par les individus qui n'ont 
eux mêmes pas un comportement détetené. Cest donc leur libre arbitre qui 
fair que parfois ilR n'agissent pas comme prévu el donc où n'obtient pas Le 
résultat escompté, 

Les lois Économiques sont donc seulement vraies en moyenne. Elles ont Ie 

caractère de lois statistiques. 

On a choisi de traduire cet éçart par une vatiahle aléatoire noté ut appelée terne 

d'erruur. 

Compte tenu de l'existence de ec terme, la loi qui va gouverner yt n'est plus une 

loi mathématique mais une loi statistique, èt. 1] faudra sc donner une distibulion 

Pour CaractéTIsEr U4. 

L'estimation des cacfficients 4, , a, lournit tes valeurs qui seront elles-mêmes 

des variables aléatoires ayant une distribution de probabilité que Fon devra 

Spécifier selon la distribution de u, 
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Qn dispose d'un ensemble d'observations porlant sur La variable expliquée et sur 
Les variables explicatives. Sachant que la telation à un caractère probahiliste, il 
nest pas possibte de déterminer de façon unique les valeurs vraies des paramètres, 


On peut simplement déterminer des estimations et connaissant la distrihotion de 
probabilité de cés éstimations, remonter vers des intervalles de confiance pour les 
paramètres. 

Les paramètres estimés sont désignés par dd, et, les connaissant, on pourra 
remonter aux valeurs vraies des parantètres @1,%5....a@7. 

Pour que celte opération d'inférence statistique Soit possible i faut qu'un cértain 
nombre d'hypothèses soient vénfées. À défaut, on aura (toujours un résultat 


numérique, mmais le résultat &sc sans valeur car il ne permet pas de porter un 
jugement sur les vraies valeurs des coefficients recherchés. 


Les hypothèses en question portent sur Les variables ainsi que sur le 1erme 
d'erreur. 

Soient ÿ1...yr les valeurs observées de y et #11 àX17,..x à ar les valeurs 
observées des variables explicarives, On peut écrire Ÿ et X sous la forme : 


Xu . . CE © | 


kru=|l. XCF,n = 
JT 
LAIE 7 
{variables explicatives} 
Les coctficients peuvent aussi s'écrire : 
u] 
al 
AUUD= | de même : LT, u= 
car] 
CFA 


On a alors 
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[@) v=xA+u 


forme matricietle de l'équation {L} 


al 
Le vecteur des estimaleurs sera AG.) = 
.i 
äi 
d'où : 
A La 
Y=XA 
Las A 
Y-Y=L 


U'est le térme d'erreur -> inconnu 
Ü est le résidu — connu. 
On distingue 7 hypothèses : 
l- Les variables sont observées sans erreur. (L'écouomètre s'‘interdit de 
s'imtniscer dans {a collecte des données}. 
On peut cependant avoir une petite idée sur ta fiabilité des chiffres qu'on 
utdise. Je degré de précision de la méthode d'estimation est dans unx certaine 
mgsure imité par a qualité des données. Il ne sert à cien d'utiliser des 
méthodes sophistiquécs sur des donnécs elles-mêmes douteuses. Leur 
estimation grossière peut suffire pour danner une idée du phénomène étudié. 


2- Le terme d'erreur a une espérance marhéniatique nulle 


: (up) = 0 
c'est-à-dire qu'on se trompe tantôt dans un sens tantôt de l'autre. L'erreur n'est 
pas sYstémalique. 

3- us à une variance a qui çst constante el finie. C'est ce qu'on appelle 
l'homoscédasticité c'est-à-dire que la marge d'erreur est la même en début ct 
en fin de période vu bien qu'elle est identique en tout point de l'échantillon. 
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C'est une hypothèse parfois contestable dans la mesure où on 4 des raisons de 
penser que la variance peul augmenter ou diminuer sur l'ensemble des 
observaiions. 


Cette hypothèse peut être tusléc ct quand on à repéré l'hétéroscédasticité on peut 
corriger ét obtenir néanmuins de bons estimatcurs. 
4- Les termes d'erreur sont mdépendants les uns dés autres, dans ke temps. 

Etuçu, j= @ sr # r 

Eux, }= oè daté 
Cette hypothèse traduit Le fait que l'erreur n'est pas systématique. En cffet. l'erreur 
en t+ 1 est mdépendante de l'erreur en L 
Cette hypothèse peut, elle aussi, être restée, Quand HA n'est pas vérifiée, on parle 
alors d'auto corrélation, ce qui eraduit souvent un problème de spécification, En 
effet, une Spécification défuctueuse est de nature à entraîner une crreur 
systématique : 
—. Onission d'une variahle importante 
— Choix erroné de la fonction (linéuire à ka place d'exponentielle par exemple}. 
Les hypothèses sur l'homoscédastivité et l'indépendance peuvent se résumer en 
écrivant : 


avec U vecieur colonne et L' sun transposé. I étant la matrice unité. 


E (U C'à est la matrice des vanançes covaciances des termes d'erreur. C'est done 
une matrice {T'X'T. 

Efayn).. Etuur j ] 

E { 4; } 


EUU")= 
(TT) 


EXuru ) Et(kyur) 


L'estinarion d'une équation unique il 


Sous l'hypothèse H4 d'indépendance des 16rmes d'erreur, si 2 terimes d'erreur sont 
indépendants E (uj uj} = 0 (à # j). Par conséquent les tormes situés en dehors de la 
première diagonale sont nuls. 


Sous l'hypothèse H3 les terines restants sent égaux à o? 


2 0 
1 
= Oukirer) 
0 oi 
Donc 


Leuvy=or | 
3- Les variables explicatives x,...x; ne sont pas linéairement dépendantes. 


Sinot, on parle de collinéarité des variables explicatives. On peut Ja déceler et 
y remédier le plus souvent. 


6- Les variables explicatives sont indépendantes du terme d'erreur. Ce n'est pas 
vérifié dans plusieurs cas : 

. Dans les modèles à équations simultanées certaines variables explicatives 
peuvent ciles-mêmes être des variables expliquées dans une autre équation 
et donc être des variables aléatoires. 

. Quand on travaille sous l'hypothèse d'anticipation rationnelle, puisque les 
variables anticipées sont supnosées résulter de l'application de la théorie 
économique et sont par conséquent assimilables à des variables 
endogènes. 

. Lorsque l’on travaille avec des variables explicatives connues avec érreur 
{remise en cause de HI). 

Dans tous ces cas on utilise la méthode des variables instrumentales. 


7- Les termes d'erreur suivent une loi normale. 


C'est grâce à cette hypothèse que l'inférence statistique peut sc réaliser car 
c'est elle qui va préciser la distribution statistique des estimatcurs à. Elle joue 
donc un rôle essentiel bien qu'elle soit difficile à vérifier. 
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TI existe des tests de normalité mais ils ne sont pas applicables sur les termes 
d'erreurs, qui par définition demeurent inconnus, 

Sachant que les termes d'erreur ont une loi de probabilité qui cst normale, les 

estimateurs c'est-à-dire les à suivent aussi une loi normale et par conséquent les 


â; - Eté; : : : 
grandeurs 2 Lée suivent une Loi normale centrée réduite. {— N(0,1)) 
à 


Mais, comme on ne connaît pas la variance du terme d'erreur ofi dait l'estimer. La 
loi normale est alors remplacée par une loi de Student Fisher : 


d; — E(à;} 


— Student(T — 1) 


di 


Cette loi permet l'inférence statistique. On peur revenir sur les vrais valeurs à, à 
l'aide d’an test simple 

HO:a; =0 
Si À, cst un estimateur sans biais, Efà;} = a, ct donc sous HO : Lia;} = 0 et 


Gi, Sudent{T —1) 


Cu, 


Soit tx = Valeur critique de la table de Student 


"a 
LL 


= 


D 
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. À; ; Sr ; 
Si —cst compris entre - «y êt ta alors on accepte IG, Le cocfficient ai n'est 
Ga 
donc pas signficanvement différent ile Q. 4 l'exiérieur de l'intervalle on rejette 


HG. 


E ou <a fl : 
Tü Fa 

+ et X; est biert une variable explicative. 
Rappel 
La plupart des propriétés qui ont été démontrées en économétrie sont relatives à 
des variables stationnaires. 


SECTION 1: L'ESTIMATION PAR LES MOINDRES CARRÉS 
ORDINAIRES DANS LE CAS GÉNÉRAL 


Soit L: modèle Y = Y A - 1 1} 
(TJ) (OU (TR 


À esl un estitüateur de A. donc 
Y= XÀ (2) 


eV. Ÿ =Ü-Y-XÂ (3 
Le principe de la méthode des moindres cartés consiste à trouver les valeurs des 
estimatears qui minimisent [a sommic des carrés des résidus. 
A- Les estimateurs simples 
On part du calcul de : 


MIN ZA? c'est-à-dire MIN EU 

Ü=Y:xÀ d'après (3) 

Par définition le transposé de Üest C'=(Y-XAY=Y'-ÀÂX (4) 
Donc, C'Ü = €". XIV -XA) (5j 


Ü'Ü =YY-Y XA-A XY+A X LA (6) 


RE il ils GI 
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Tous Les éléments du produit sont des scalaires. Sachant qu'un scalaire est égal à 
son transposé Y'XÀ = Â'X'Y et le résultat peut s'écrire indifféremment. 


UU=YY-2F XA+A X'XA 
1 D bar 
ou 
A À # La ñ 
D'U=YY-2A X'Y+A'N'XA (7} 
«LD LAS IL) 
On cherche le minimum de celte expression par rapport à À. Cela revient à 
annuler ja dérivée de C'Ù par rapport à À : 


MINGO D) = EE 0 
d A 
QUE 3 xry 42 XÀ = 0 
gA 


d'où X'Y = X'XA. Si X'X est une matrice régulière, adtnettant une 
inverse. À S'obtient en multipliant l'équation précédente par (X'XY L : 


d'où 
| À ={X'XYISY (8) 
Remarque 
Dans le calcul précédent. on a : 
XX Te 
da A 
9 À X'XÀ 2x x 
gÀ 


Ce résultat se vérifie aisément cn explicitant les éléments du calcul. Appetons 
&,%2..0, les éléments du produit X°Y. 
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@] 
(9) Soit X'F= 
41) 
œ 
En multipliant X'Y par À', il vient : 


a 
À'CX' 1) =...) 


a; 
= â @ + 50... 410 (10) 
La dérivation de cette expression par rapport à À revient à dériver par rapport. À 
chacun des éléments de À, sait &... Àj 
da + dde +... +de) | 2 | 
| 
dû i 
a. . 
BA CAT) 
d4 
dé: +5) +, + æ, } ns 


dû 


€; : 


Le résultat est un vecteur (l, 1) composé des termes Gz à O4, c'est-à-dire X'Y lui- 
même. 


À 
De mème, appelons Bjles éléments du produit X'X de format (LA). Par 
construction XX est une matrice SyMÉtrIquE, 


| Économétrie appliquée 
By + - - Bi 
{li) X'X= | | avec B;=B; Viet j 
Bu + 8, 
En multipliant à gauche par À’, on obtient un vecteur Ligne à { éléments : 


AUX X = (fn + da +. db Br +. deBar +. +4Br) (12) 
Puis en multipliant à droite par À, il vient 

Léh Au1 + Sa faute di ht Dé + +6 Pr + Ga fort #8 Br à 

(1%) qui ext un scslaire, 


La dérivée de À'X'X À par rapport à À s'obtient en dérivant par rapport à âi… 


[aa X'XÀ 
TS = 2Biû +242 +.—2âlhr 


DA XXA 
24 


| ee” = 23 + 2dB5:+...- 2818 
L dây 


AB + 42Bas +... + 
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Après avoir calculé À, nous pouvons examiner les propriétés de cet estimateur. 


a) À est un estimateur sans biais si E(U)=0 (H2) 
On sait que ÀÂ={XX)!X'Y 

et que Y = XA +U 

En teraplaçant Y par sa valeur : 
A=(XX)!X'(XA+U) (14) 

Â=(XX) XXA+(XXYIXU (15) 
Â=A+(XXYIXU (16) 

d'où E (À) = A + (XX)! X'E(U) 

donc E(Â) = A si E(U) = 0 

Sous H2, À est un estimateur sans biais de A. 
b} Variance de À : 
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Ea variance du vecteur À est donnée par la maltice des variances covariances 


obtenue en développant en lignes et en colonnes : 
v)= u- EA)XA- roy | 


[ Var Cove + : va | 


vtÀ)= 
1) 


Vs - | 
COV (à) to ee Ve | 
On sait que À = (XX)! X'Y et que E(Â}= A. 
Donc À - EÂ)= (XX)! X'{XA + U)- A 
À-HÂ)= (XX)! XU 
et[Â-EA) =UX(XX)t (17) 


(X'XY! matrice symétrique par construction est égale à sa transposée. 
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D'où V{Aj = EfX'XY I X'UL' XX XI IT 18) 

Comme les éléments en X sont des éléments certains, 4 peut encore écrire : 
V(A)=(XXTIX ECL") XEXXT! (19) 

Or, sous Les hypochèses HT ét F4, on à É(EU" = o?1.c étant un scalaire et Fa 
matrice unité de dimensions (T, ‘Ti. 


Donc V{Â) =o2iX x) X'AX{X X) ! (20) 
soit finalement 


VGA) + oc XX) | 

Remarque ; La simplification n'est possible que si HA et H4 sou vérifiées er 
donc si en peut écnire LCL") = gi 

On peut montrer que sons les hypothèses qui ont été posées, l'estimatcur À des 
moitdres carrés ordinaires (MCO) est un estimateur à variance mimimale, dans la 
classe des estimateurs Sans biais. 

Si maintenant EtUU'] # oi. mais E{(UL'";} est connue, par exemplé E{UU") = Y, 
il cst passible de retrouver los propriétés précédentes comme on le verra dans la 
suilt dé ce chapitre, à La section 5. 


Le calcul de V{Â) fait appel à la vanñance du terme d'erreur, oi. 


. . TC : ? . 
Celle-ci est inconnue. On peut cepeodant utiBser un sstimateur de &% construit à 


L 3 ; 3 
parûr de la varranve du terme résiduel Gi. 


B- Estimations sous contrainte ct multiplicateurs de Lagrange : 
Les estimatcurs qui ont té présentés jusqu'ici présentent la propriété de minimiser 
la somme des carrés des écasts entre valeurs calculées ct valeurs observées de la 
variable expliquée y, On peut leur imposer de satisfaire également d'autres 
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propriétés. On obtient alors des estimaleurs sous contrainte qui minimisent 
toujours la somme des carrés des écarts, tout en satisfaisant certaines contraintes. 
Ils sont normalement “moins bons” que les estimateurs sans contrainte et ceci 
pose la question de savoir si la ou les contraintes sont acceptables compte tenu du 
modèle et des données utilisées. 


Exemple de contraintes : on peut imposer que les rendements d'échelle d'un 
processus de production soient constants. Dans ce cas si les coefficients à estimer 
représentent les élasticités de lu production par rapport aux facteurs de production, 
cela se traduira par la contrainte Ya; = 1. 


On supposera dans la suite que les contraintes sont de forme linéaire. 

Les restrictions de formé hinéaire peuvent être écrires sous forme matricielle celles 

que : EL ASh 

— À = vecteur des estimateurs sous containte. C'est un vecteur colonne à 
léléments. 

Pour Le distinguer du vecteur À des MCO on l'écrira À 


-  H = matrice des restrictions de format (q x 4), "dq" étant le nombre de 
restrictions. 


—  hesr done un vecteur colonne de format {q, 1}. 
Exemple : avec 4 coefficients et 2 cuntraintes 
t=4 

qa=2 

& +3 =] 

Lada + Haba = 5 


On à donc : 
ä 
k 1 0 © -[] 
0 0 &; “le az 
da 
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Dans ce cas, on nc peut plus calculer les cstimatieuts selon la méthode qui à été 
précédemment utilisée, on les calcule en utilisant un Lagrangien qui s'écrit : 


L=(y -À' XUY— XA)+24 (HA h) (22) 


carré des 
On choisit pour des raisons de commodité d'écrire le multiplicateur de Lagrange 


sous la forme + 2 À. 


L=FY-YXA-À XY+AX XA+21' (HÀ-h) (23) 

Tous les termes de cette équation sont des scalairts ot X' est donc par conséquent 
un sector ligne à q éléments. 

On peur utiliser la propriété qu'un scalaire est égal à sa transposéc pour Y'XÀ et 
À (HÀ - h). 

Il vient : 

L=YY-2AXY+AXXA+2(ÂH'-hX (24) 

L'application de fa méthode des MC conduit at calcul des conditions du premier 
ordre de l'équation de Lagrange : 

2.05 

24 

dL : 

—=0 (26 

Sa nu 


He _ JXYIIXXA+IHA=O 


dÂÀ 
dE 
À 
En prémultipliant par H GX}! l'équation (27) 
HIXXYI XY 4 HIXX)TN'XA + HOXXYIHA=U (25) 

Mais (XXI X'Y n'est autre que l'estimateur des moindres carrés ordinaires sans 
contrainte, À. De ce fait, on peut écrire : 

-HÂ+HA+H(XX)IHA=C (30) 


=AH'-h=0 (28) 
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Soit 

HÂ-HÂ=H(XXY IH (31) 

Mais d'après la contrainte, HA =h 

donc 

HÂ-h=H(X'xY' HA 

soit en raultipliant à gauche pac [H €X'XY HT! 

À=IH(XXY HT IHAÀ -h] (32) 

Cette valeur de À peut être reportée dans Les canditions de premier ordre : 
-XY+XXA+HIHIXX} HI [NA -h}=0 (33) 

En prémultipliant par (XX) 1, il vient 

(XXYUXY + OX'XTUXXA + (XX) HITKX'X5! HI][HÂ-hn]=Ù (34) 
À +A+(XX) IT [HXXF HI[HÂ -h=0 (35) 

et 

A=Â-(XX)! HIH(XX)! HIIHÂ-h} (36) 

On constate que l'estimaicur sous contrainte ne peut pas coïncider avec 


l'estimateur des MCO. Cela revient à dire que l'imposition de contraintes dégrade 
la qualité des estimateurs, 


On peut tester la validité des contraintes qui ont été imposées sur les estimateurs. 
C'est le principe des tests de résiriction à prion. 


2. est un vecteur ayant autant d'éléments que de restrictions qui mesure la perte ou 
le coût associé aux restrictions à priori. 


Il y a un mulüiplicateur par contrainte dent l'interprétation est identique à celle des 
multiplicattues de Lagrange habituels. 1L mesurc le coût en termes de fonction 
objectifs, donc ici en termes d'augmentation de Ja somine des carrés des résidus, 
dc Ja contrainte à faquelle il est associé. 


Il existe toute une série de lests sur la base des propriétés des multiplicateurs de 
Lagrange. Ces tests uni l'avantage de n'utiliser qu'une seule régression alors que 
deux seraient normalement nécessaires pour un {est de restrictions a priori. 
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Soient Z4,? la somme des carrés des résidus de l'équation sans contrainte et 
Es? la somme des carrés des résidus de J'équation avec contrainte. Par 
définition 24? > Et. Un test de restriction a priori classique passe par la 
comparaison des variances résiduelles, 


Eh-ER? 


Le rapport : (47) Fe= uit ne loi de Fisher {gT—-1} 
DD 7e 
T-1 


En posant : 
HO : Îcs contraintes ne modifient pas significativement la vatiance résiduelle, Le 
modèle avec contraintes est acceptable, 


H1 : Les contraintes modifient significativement la variance résiduelle et no sont 
pas compatibles avec les données. Le modèle avec contraintes est rejeté. 


SiF empirique < Ft valeur critique de la table de Fisher, on conclut HO. 

Si F empirique > Fa on conclut H1. 

Mais on sait par ailleurs que la loi de l'isher F {n, m) tend vers la loi du É 
centrée lorsque m augmente. Le rad centré {ou x° réduit... où ‘augmenttd x L 


2 
en anglais} étant égal à Xn 7. : 


2 
F{n,m) sf qd #1 — =. Soit encoren F{n, m)—. y? 


En régressant les résidus de la régression sous HO sur les contraintes, on obtient 
un R? qui mesure le pouvoir explicatif de [a régression, done ici La part de la 
variance résiduelle {sous HG) imputable aux contraintes. En utilisant l'expression 
de Fe donnée dans (37), il vient : 


Fe = T-i 


R? soit{T-DR?=q Fe. 
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Pour T grand cette dernière expression peut être approchée par FR? a gauche et 
par * à droite. 


Le test de restriction à priori peut donc se ramener au test du 1° pratiqué sur le 
R2 de la régression auxiliaire. 

Soit xè la valeur critique de la table. la règle de décision est : 

x? empirique < xÈ on cunelut HO. 

x? empirique > 1È on conclut H1. 


Les tests sur les multiplicaleurs de Lagrange, encore appelés LM tests dans La 
littérature peuvent être appliqués à un ensemble très diversifié des problèmes en 
fonction de la nature des contraintes figurant dans l'équation (32}. Ce sont en fait 
des cacres dans lesquels on peut l'aire rentrer toutes sortes d'hypothèses. On en 
verra un exemple dans la section 3, 

Dans ta réalité, les conditions nécessaires pour l'application de la méthode des 
moindres carrés ordinaires ne seront remplics qu'exccptionnellement. La situation 
la plus courante est eclle où unç où plusieurs hyparhèses ne sont pas vérifiées. H 
est essentiel], tout d'abord, de s'en apercevoir puis d'anpliquer les procédurcs de 
correetlon appropriées qui permettront aux estimateurs des moindres Carrés qui 
sont à La fois parmi les plus simples à calculer et les plus robustes de Conserver 
Keurs propriétés. 


Les sections qui suivent sont consacrées aux principaux problèmes susceptibles de 
se poser lorsqu'on cherche à estimer une équation unique. Le problème de la 
stationnarité des variables, qui est en lui-même un problème essentiel], susceprihle 
de transformer compiètement l'approche de l'économétrie sera traité plus en détail 
aux chapitres à ct 4. 
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SECTION 2: LA COLLINÉARITÉ DES VARIABLES 
EXPLICATIVES 


À 


Dans cette situation, les variables explicatives ne sont pas linéairement 

indépendantes. Cette collinéarité peut être stricte", se traduisant par une relation 

“déterministe" du type : 4j xj,+ a xj = k 

Cette collinéarité peut être seulement “approchée”, auquel cas la relation n'est plus 

détérministe mis plus large. 

Dans tous les cas c'est l'hypothèse HS qui ést prise en défaut, (in cammencera par 

exposer les conséquences de la coilinéanté. puis on indiquers les techniques 

permettant de la mettre ex évidence, avant d'indiquer quelques transformations 

usuelles. 

Conséquences de la collinéarité 

En cas de collinéarité stricle 

C'est une situation cxtrémement rare ot qui correspond généralement à une erreur 

de conception du modèle au de définition des variables. 

Dans ce ças l'ordinateur affiche le message CRROR éar (X'XT ! n'existe pas, En 

effet, la matrice X aura au moins deux colonnes linéairernent dépendantes, X'X 

aura au moins deux lignes et deux colonnes linéatremunt dépendantes et sera ne 

matrice singulière {det = Q). [l est donc impossible d'en calculer l'inverse. 

En vus de collinéarité approchée 

Dans ce cas, la matrice X'°K n'est pas singulière, L'ile tendra vers la singularité 

c'est-à-dire que son déterminant tendra vers (. Ceci va impliquer que les valeurs 

numériques de {X'X)! seront très grandes puisque le déteratinant tend vers 0. 

I! en résulte les conséquences suivantes : 

— Les cstimateurs tendront à avoir des variances élevées 

= On risque d'avoir des estimateurs “instables”, c'est-à-dire que si on change 
légèrement les données, (Exemple : modification de fa période), Je résultat 
nütnérique est susceptible de changer fortement. 
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B- Comment faire apparaître ka présence d'éventuelles coilinéarités 
entre les variables explicatives ? 
La méthode la plus simple consisté à calculer Le tableau des corrélations partielles 
entre les variables explicatives. 


Comme ce tableau est symétrique, on en donne généralement là partic supérieure. 


Ce tableau ne donne pus une réponse claire car Le caefficient varie entre -| ét 1. 
Flus on se rapproche de -1 ou de 1, et plus le risque de collinéarité est fort. [] faut 
donc se donner des règles de décision en tixam des valeurs critiques. 


Ex: SiR < 0,7 — pas de problème 

SiR>0,9 > présomption de collinéarité. 
Mais, de tableau ne garantit pus contre toutes Les formes de collinéarité, 1 
renseigne sur la collinéarité éventuelle entre deux variables mais pas sut la 
collinéanté entre une variable et un groupe de plusieurs variables, 
On ne sai don rien d'une situation du type : 


Cl = EE je + Cuxps + À 


La variable i est en relation avec la somme des variables j et h mais ni avec l'une 
ni avec l'autre des variables, prises indépendamment. 

IT peut donc y avair des situations aù it n'y à pas de cursélation apparent entre 
variables explicatives mais où il ÿ 8 quand même callinéanité. Ces situations sont 
heureusement assez exceptionnelles, 
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De façon yénérale, dans toute situation de collinéarité on rémargue une 
mcohérence enlre les différents Lèsls qui permettent d'apprécier la qualilé d'une 
TÉrCSSIon. 

Souvent. il ÿ a une bonne qualité d'ensemble c'est-à-dire un R? élevé) associé à 
des coctficients qui ne sünl pas significatiés, La raison est que cs variances sont. 
1rès élevées &t donc Ios tests ne vafident pas Les cocffisients. 

Symétriquement, la contradiction éventuelle entre une bonne qualité d'ensemble 
de l'ajustement ct le caractère non signilicuif de ecrtains coefficients de 
régression peut constituer un indice révélant la présence de cullinéarité. 
L'estiation en présence de collinéarité 

« En ças de collinéarité stricte, c'est souvent lié à un problème de définition de 
variables ou à l'introduction maladroite de va:cables muettes {dummies), 

Les variables muettes sont des variables qui mesurent l'effet particulier lié à une 
périude. Elles prennent La valeur ÔÜ où la valeur L. Lorsque l'observalion 
correspont à la pérniude. elles valent 1. Elles prennent la valeur O0 dans les autres 
CAS. 

Exempie 

Si on s'intéresse à la consommation d'électricité en données trimestrielles (y,). On 
va expliquer éètle consommation par le revenu R,. Mais on sait aussi que ja 
censommalion d'électricité est influencée par le climat. Ja température, la 
luminosité, la longucur du jour, etc. 

Des variables mueties traduisent l'influence des saisons sur la vonsommulion. 

À pour l'automne, 

H pour l'hiver. 

P pour le printemps. 

E pour l'été. 


Un modèle très simple cxpliquera la consommation du trimestre par le revenu 
correspondant ec l'effet saisonnier. 


=ap+éR, Has A +aill+asP+asE+u,  (H 


Hal rence Gé autns se bles iniélies 
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On suppose qu'an prend la période IV 1988 -— III 1993. 


La première observation correspond à l'Automneé 


À = | si Automne (A = Ô sinon) 
H=0û 
P2<0 


H P 
QE 
| | 1 0 DA à 
6 6 |  Varisbies 
j (ll | 0 1 muettes 
ü Le, 0 saisonnières 
| Nora 
Ÿ | : (6! 
| ts 0 | | 
DER AR EIRE NI 


On remarque qu'à chaque période an à À +1T+E + P 21 Ÿ,. Comme la 
régression (1} comporte un terme constant, 4 matrice X contient une colonne tte 
termes égaux à 1. Par conséquent, la somme des quatre colonnes constituées par 
les variables À, H, P et E ei éyale à la céionne correspondant au terme constant. 
Les colonnes de X ne sont plus indépendantes. 


On a une relation de collinéarité stricte entre les vaciables et (X'X} ! ne peut être 
calculée, 


Par conséquent quand on introduit des vañables saisonnières suus cette 
forme, il faut laisser 40 minimum un degré de liberté. 
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Autre exemple : on s'intéresse aux performances d'un échantillon d'entreprise et 
on veut tenir compte de l'effet pussible du cype d'activité. Les entreprises sont 
classées dans l'un des 3 secteurs : 

. biens d'équipement: 
. biens de consommation; 
. biens d'investisstment. 

er une variable mucite mesure l'effet du secteur, On retrouve exactement le même 

problème, si l'équation estimée comporte un terme constant. 

Les solutions 

1. Dans le cas de collinéarité stricte 

— Supprimer l'une des variables “dummies" qui correspond à L'effet Le plus 
faible. ce qui revient à choisir une variable comme référence et transformer les 
autres de sorte qu'elles mesurent les écarts par rapporl À cette Variable de 
référence, 

Si on prend À comme référence, on a : 

Se = do + à Re + ap CA+TIHPHE) + 143 - 42) H + (ag - 13) P + (as - as) E 2) 
et! y, =cot+aR +(@-)H+(as-a)P+{as-@iE 

On mesure les effets différenticls de H, P, E par rapport à A qui sert de 
référence. 

2. Dans les cas de collinéarité approchée, on dispose de toute une panoplie de 

solutions possibles qui passent par des transformations de l'équation. 

On peut faire la somme dé 2? variables (quand ça à un sens} 

Autre possibilité : normer l'équation c'est-à-dire diviser toutes les variables par 

l'une d'elle. Cela permet d'éliminer le trend c'est-à-dire la tendance commune aux 

différentes variables qui peut être à l'origine de la collinéarité. 
Exemple 
Q=ag+aiYi+aietr (1 
C'est une vieille formulation de la fonction de consommation dans 
laquelle Lé revenu et la consommation des ménages suivent les évolutions 
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très proches en données chronologiques. On peut La transformer en 
divisant par ÿ, Où Ci; 
Dr +a + a SL (4) 
Yr  Y Yr 
La signification des coefficients est celle du modèle du départ( 1). 


On peut aussi normer par iL revenu une fonction d'investissément ét la variable 
expliquée sera un taux d'investissement 1/Y, 


En normant par Le capital on obtient 1c taux d'accumulation T/K:. 


On peut aussi transformer les variables en taux de croissance ce qui supprime les 
effets tendanciels. Mais ce n'est possible que si la Loi qu'en cherche à estimer peut 
être assimilée à une loi à élasticité constante : 


Exemple : 

— Pour les fonctions de commerce extérieur ; 

= Pour les fonctions d'emploi. 

On peut aussi utiliser certaines propriétés de la fonction : 
Soit La fonctian de production Cabb-Dougias 

Di= AUTRE, (5) 

log Dr=log A+aiog £, +Blog K, (6) 


Marc Guillaume donne l'exemple suivant{2) estimé à partis de la fonction de 
production linéarisée : 


log C4 =0,28logL, +0,52logX, +C A°=0,984 (7} 
Les écants-1ypes des cocfficients sont indiqués entre parenthèses. 


C'est un exemple typique de collinéarité : le pouvoir explicatif est bon mais 
aucune des variables n'est significative. 


ME Remarque : si Le terme d'erreur de l'équation de départ est homoscédastique, la 
transformation proposée est susceptible d'introduire de l'hétéroscédasticité dans (4). 


(2} Marc Guillaame “Modèles économiques” Presses utivérsitaires de France 1971 page 99. 
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Dans le cas particulier ci-dessus il y a une ahsurdité supplémentaire. Si on somrne 
a + B on trouve Ü,80 (cas de rendements fortement décroissants), C£ qui est 
inacceptable, compte tenu de ce que l'on sait de l'évolution économique de Ta 
période sur laquelle porte l'estimation. 
De plus, très souvent pour estimer les fonction de production on se base sur les 
parts distributives. Si les rendements ne sant pas constants cela ne marche pas car 
c'est en contradiction avec l'élaboration des données, 
En posant par hypothèse que Les rendements sont constants, tl est possible de 
transformer l'éguation qui devient : 

K 
Leg Æ 2(1- og +e.  (8} 


à us K 
Les estimations donnent : log É = 026108 +6 (9) 
PE + 
‘&=0,74 


i 
18 = 0,26 
La transtormation de Ja fonction permet danc d'éliminer le problème. 
Mais parfois, aucune des transformalions précédentes ne donne de buns résultats. 
On peut alors utiliser une méthode d'estimation par balayage qui va Consister à se 
donner des valeurs plausibles pour l'un des coefficient afin de caleuler les autres. 
Exemple 
Yt=aikittaxittas (10 
On se donne a, et aux valeurs de a, vont être associées celles de az et de 
43. 
On va calculer une nouvelle variable : v, —a,x4, 
Yu = 2x, +4, et on régresse pour obtenir 3 et 3. 
On fait çe travail pour un ensemble de valeurs de à, ct on retient celle des 


valeurs de à, qui minimise la somme des carrés dés résidus ainsi que [es 
Valeurs de 2 et de a; qui lui sont associées, 


Mais cette procédure par balayage pose plusieurs problèmes : 
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— Le choix du parainèur : est-ce à] wu a3 ? 

Ça dépend de la signification des variables et de leur importance respective, 

— Los différents critères qui caractérisent te réenltat de la répression ne sont pas 
toujours convergents, | peut arriver que là somme des carrés des résidus 
dininue tandis qu'un cocffiient cesse d'être significatif. Ceta suppose que l'on 
définisse des critères permeltant de sélectionner la solution la plus 
satisfaisante, 


SECTION 3 : L'AUTOCORRÉLATION DES TERMES D'ERREUR 
C'est une situation dans laquelle les termes d'erreurs ne sont pas indépendants au 
coufs du temps. 

On va introduire une hypothèse selon Iaquelle le terme d'erreur dépend 
Üinéairement de sa propre valeur passée. Cette hypothèse se traduit par : 
ut= Put] +vt(l} 
P 2 coefficient d'autncorrélation 
vi= nouveau terme d'erreur qui satisfait aux propriétés habituelles 
c'est-à-dire 
Étvi=0 VYr 
Variance de (V)= of 
Le terme d'erreur de l'équation comprend donc une partie systématique el unc 
partie aléatoire. 
v£ est une partie purement aléatoire 
8%. est unc partie systématique 
{É} est une relation possible mais n'est pas forcément la seule car : 
— Ellsest linéaire 
— C'est une autocorrélation du premier ardre car ctle ne dépend que de la 
pénude antérieure, 


{autocorrélarion d'ordre 2 : ug = P] Ut-] + P; 4-2 + vt) 
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Les différentes formes que peut prendre la relation d'autocorrélation peuvent 
Être testécs. 

L'analyse du signe des résidus peut donner une première information sut 

l'existence où non d'un éventuelle auto-carrélation : plus on observera de 

régularité au niveau des signes d'erreur, plus on vérifiers l'équation {1}, 

— Si on a une alternance +/-, on risque d'avoir une auto-corrélation négative 
c'est-à-dire ? < 0. 

En ettet cela traduira lé fait qu'à une erreur positive succède une erreur négative et 

inversement. 

— Sion a une succession de résidus du même signe suivie par une succession de 
termes de signe opposé du type : +++ ---. ++ +, une erreur est suivie d'une 
erreur de même signe, ce sera le signe d'une auto-corrélation positive P > 0 

L'auto-corrélation est un problème stastisique qui peut aussi renvuyer à la 

spécification de l'équation car dans certains vas, elle trouve son crigine dans une 

spécification inadaptée. 

Supposons que la vraie fonction sait de type parabolique Ÿ £t qu'on l'ait estimé 

par une fonction linéaire. 


résidus 
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La fonction estimée, Ÿ, va d'abord sous-estimer la variahle expliquée, v, puis la 
surestimer, puis SOUs-estimer à nouveau. On observera une succession de résidus 
positifs, puis négatifs et ensuite positifs, 


Done, toute situation d'auto corrélation appelle à une réflexion sur le modèle. (Si 
il y à un cûté systématique dans l'erreur c'est qu'il ÿ a une cause porinanente). 


-  Lne distribution non systémarique des signes des résidus est compatible avec 
une absence d'auto-corrélation. Ce n'est évidemment pas un critère très précis 
et il existe, heureusement, des tests plus performants. 


À- Le diagnostic de l'auto-corrélation 


On peut tout simpiément régresser à, sur à, êt tester le coefficient de régression 
de l'équation { 1) : 


ût= Püi+v 
H0:/=0 
HL:P +0 


Ce n'est pas une méthode très performante. On lui préfère Le test de Durbin- 
Watson (DW). 


— Test de DUÜURBIN WATSON 


T 
A < + 
> ch) 
2 
T 


DE = EE — 


DAC 


I 


F7 
Du) 
L 
T 
Say 
l 


Pour ua grand nombre d'observation c'est-à-dire quand T — + 0, 


ou 


DWw- 
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on monire que : 
lbw 2.20 | 


En cffcr : 


Ecs 2 premiers rapports tendent vers 1 car numératcurs et dénominateurs ne 
différent que par un ere, 43 OU Uy respeetivement. 


Le 3° rapport tend vers À l'estimateur des MCO dans f'équation de définition (1) 
F 


DETTE 


Puisque : = 


T 
À 
> 
[ 


Cetie propriété permet d'établir & domaine de variation du test de Durhin 
WATSON 
Pour p =1 DW -»0 
Pour p =: DW - 4 
Pour p = 0 DW - »2 qui cst la valeur centrale de ce test. 
Plus de tust est faible, plus on à de chance d'uvoir f > G (c'est-à-dire de 
l'auto-corrélation positive) 
Les valeurs critiques du test de D'W tenamt compte du nombre d'observation et du 
nombre de variables exphrçatives om 66 calculées, Le problème est que ces 


valeurs critiques ne sont pas uniques. | y à en fait un seuil inférieur et un seuil 
supérieur, Entre les deux, i[ v a une zone d'incertitude. 
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La table de DW est sculement calculée pour les auto-corrélations positives. Les 
autres sc déduiront par symétrie. 11 cst donc prudent de retracer loutes tes 
situations possibics. 


D p>0 dl ? du p-0 4-du? 4-dl p<«04 
—————#——#>——+—+— 


«] = seuil inférieur 
du = seuil supérieur 


Remarque : La zone d'incertitude diminue quand le nombre d'observations 
augrnente et elle augmente quand le nombre de variables explicatives augmente. 


Dans les tables courantes, elle est maximum pour les valeurs suivantes : 
T=l5 5 variables explicatives 
à] = 0,56 
du=2.?l 


—  Îl existe une autre table qui à été calculée par THEIL et NAGAR avec un 
sæuil unique {ce qui élimine là zone d'incertitude). On peut donc s'y reporter 
en cas de doute. Toutefois cètte table n'est valable que sous certaines 
conditions. 


—  Qn peut tester des formes plus complexes d'auto-corrélation à l'aide du test de 

Breush - Gudfrey qui appartient à la catégorie des tests sur les muitiplicateurs 
de Tagrange (LM test) présentée dans la section 1. La régression auxiliaire est 
ici constituée par les résidus de l'équation avec contrainte qui sont régressés 
sur les variables explicatives et Les résidus décalés. Elle permet de tester la 
significativité globale de l'ensemble des coefficients d'auto-corrélation. 
Le modèle sans contrainte est celui qui contient des auto-corrétations. Le 
modèle avec contraintes est celui qui n'en contient pas, puisque cette absence 
est équivalente à poser que les coefficients d’auto<orrélation correspondants 
sont nuls. 

Exemple 


Soit le modèle général Y = XA + U sur [lequel on désire tester une auto-corrélation 
d'ordre 3. 


36 Économétrie appliquée 


(2) 4 = Pi, y Pate _z + Pau a+ 
Pi 

et P=|PM 
Pa 

Le modèle sans contrainte est : 

(3) F=XA-Zp+4Vv 

où ds (a. dr) 

Le modèle avec contrainte est : 

(4) Y=XArU 

Les cuntraintes sont : 


PA1=0 fP2=0 p;,=0 ou de façon équivalente h: =0 
(3,1 


Les hypothèses sont donc : 

HG: p =0 Le modèle contraint est acceptable, 

H1:p #0 Le modèle contraint n'est pas acccptahle, 

Soit À Le vecteur des résidus du modèie avec contrainte, modèle correspondant à 
l'équation (4). 

La régression auxiliaire consiste à régresser : 


{5} Ê = xA+ Z2p+V où Z est la matrice des résidus décalés de l'équation (4}. 
Soit R? le coetficient de détermination de cette régression. I! mesure le pouvoir 
explicatif des résidus décalés. 

Le test consiste à calculer TR? et à comparer avec Fe pour le seuil dé confiance 
retenu. 

Plus les résidus sont corrélés entre cux, plus Le R? de la régression a des chances 
d'être élevé er donc plus on a de chance de dépasser les seuils critiques du y? {car 
plus TR? sera élevé). 
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Si TR? > 2 — autu- corrélation — H1 
Si TR? < xà — indépendance — HO 
B- Les conséquences de l'auto-corrélation 
Les conséquences de l'auto-carrélation sont : 
- E{UU'}# 0°/ 


En effet, les lermes d'erreurs n'étant plus indépendants au cours du temps la 
matrice des variances covariances ne peut plus être une matrice diagonale : 
Fuur) #OSsitat 


On traduit l'auto-corrélation par 
H =pu;_j+v (D 
Cette relation permet de calculer les covariances entre les termes d'erreur : 
Etuu,. )= Eltpu,.r +v,)4, 1} (6) 
= pos + Eva 1) 


2 
= 26, 


Frs =pos (7). 
| 


Ur = Pl 2 Fr) 
u=p[pu. tv. i]rs (8) 


2 ; 
= OU + Pr + 


Less. 


Elu 2)= E[(pfu-2 +pvir + vou. 2] (9) 


2,2 L 
sp"o, +DE(v,_ 14, 9 )+E(vu, 
ru P Ne 7 2) Ga 2) 


Donc 


| Etqmy)=poi | 
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soit 


ie 
Eu.) Ha | 
Remarque : Les termes en v sont postérieurs aux termes en u et ne les 
contiennent pos donc les deux variables sont indépendantes. 
Si ef v avart des termes en y antérieurs aux termes en u on n'auroi plus Île 
droit de conclure à l'indépendance car quand on développe uz selon la 
formule : 
u, "x DÜPu,_3 +t.)+v (8) 
= p? ; , 

Hu =D" 2 PV FYe 
u, CONTIENT V,.} 
On constate que up} contient Vis 
u, se "souvient" des u,; or danx Ces up;it y a les v, de la formule de 
définition ue = Puf + Ve 
— les une sit pas indépendants des v,; mais les v, sont indépendants 
des lt-j- 

Exprimons u, en fonction des v,, voi … 

u= put 

Uri s D U2 ve (10) 

Ue2 = Plra + Vis 


We P (pue + Un + 
ue= POP (Puce + Veai+ Ven Er 
= p° Ua + p? vez t Pvir+ ve UD 
En substituant de proche en proche, an voit qu'il est possible d'exprimer u, en 
fonction des termes vt décalés et des puissances successives de p. 


| He — Zp'v,; 
Si Efv,)=0 vr alors Ffu,}=0 Yt 
Les termes d'erreuts restent donc d'espérance nulle. 


Conséquences sur la matrice des variances covarianees E(L'U') : 
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. 


2 


Tia COvHH COV Hylir | 
EÇUU" )= . à (12) 
? 
COVU,H] Tir 


La matrice n'est plus diagonale, mais les termes situés sur la diagonale restent 
égaux entre eux tant que l'hypothèse d'homoscédasticité reste vraie, On peut dont 


à ; ; ; L ; 
mettre ojen facteur puisque tous les éléments de 11 matrice s'expriment en 


: ? : 
fonction de aÿet des puissances de p. 


1 p p2… PUR 
d-2 
PT PP …p : 
E(UU }= ©? = o°Ww 
PU ’ 
Re pr 
CRE CR OT 
avec W= 
T-L | 


Cette matrice W fait que l'estimation par les MCO n'ese plus correcte puisque 
E(UL'"}# ci £. De plus, il faut distinguer deux vas de figure : 

—  Équation sans variable endogène décalée 11). 

—  Équation avec variable éndogène décalée {2}. 


Dans un modèle à une équation, le terme vartable endogène désigne la variable 
expliquée. On donnera dans Le chapitre suisant une interprétation plus complète, 


Ainsi l'équation vf = 4jX1, + &xa + LU, ne contiènt pas de variable endogène 
décalée. 
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L'équation yt= ax, + zx + Pire + U, contient y.1, variable endogène décalée. 
+ Équafion sans endogène décalée 

Le biais résultant de l'application directe des MCO porté uniquement sur la 
variance des Cstimateurs. c'est-à-dire que les estimatcurs cux même demeurent 
sans biais : ECÂ) = A 

Preuve 

La matrice des variances covariance n'intervient pas dans le calcul de vés 


estimmateurs et donc ils n'ont pas à être affectés. Par contre la variance de À va être 
biaisée. 


(À) = EX OLXUU XX x] 


avec E(UL")= o2W 
Ti vient : 


V{A)= o CX'XTUX W XX XY 


alors que l&s MCQ calculeront : 


V(A}= o2(X' 41! 


Quet est Le sens du biais introduit ? 
— En cas d'autv-corrélation positive c'est-à-dire 9 > 0 
L'ordinateur utilisera La matrice ? au lieu de la matrice W 
Si p > 0 alors on va généralement sous-estimer la variance V(À) 
Cela impligue un risque d'accepter pour Yalable un coefficient qui n'est pas 
significatif (On surestime effectivement les tests de Student}. 
—  Encas d'auto-carrélation négatives c'est-à-dire p < 0. 


Les termes nan diagonaux de Ÿ sont tantôt positils tantôt négatifs et la 
conclusion est moins évidente, 
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Remarque : Le fait que les estimateurs soient sans biais EtÂ} = A ne doit pas être 
interprété comme un signe d'exactitudk mais simplement comme le fait qu'ils ont 
autant de Chance d'être biaisés dans un sens que dans l'autre. 
Ea particulier. quand on estime une Équation sanx endogène décalé par les MCO 
ct par une procédure de correction de l'anto-corrélation an ne trouve pas les 
mêmes résuMats numériques. Cette différence est normale et pravient de Ja 
différence dans les méthodes d'estimation. 
Donc. méme si il n'y à pas de variable tndogène décalée en présence d'auta 
corrélation il faut quand même Utiliser une procédurs appropriée. Par contre si il 
n'y a pas d'autu-corrélation, on n'a pas le droit d'utiliser ue méthode d'estimation 
comportant une procédure de Correction. 
Quelle relation y a c-il entre la variance de u €t la variance de v ? 
dé = v, + Dés jp, 

re 2127 2 | 
Eu) =Etv +p? po, ,?} (9 


Eu) = ovil + pt +. +p 


d'où 
: : 
= 0" 
1-p° 


Ce. qui implique que la variance dé u est finie si et seulement si jpi<1 4 


+ Equation avec endogène décatée 
(4) y = aa, +5, +45 bu, 
OÙ nm =pu 17 
En décalant d'une période l'équation 14, il vient : 
(15) Men) = rt À bye-2 TA +, 
ve variable explicative de (14) contient ct n'est pas inlépendante de 4.1 
Or, d'après (1), u, le rime d'erreur de l'équation (1} contient lui aussi u, 
Donc, dans l'équation (14) où a unc variable explicative y qui n'est pas 
indépendant du terme d'erreur u.. 
[sa et LL, ONE CL COMITNIN Ut. 
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L'hypothèse H6 d'indépendance des variahles explicatives et du terme d'erreur 
n'est donc plus satisfaite. 
Dans cé second cas on à à la fois un hiais sur les variances des esttimaleurs mais 
également un biais sur les éstimatèurs ux-mêmes 
Un exemple très simple permet de le vérifier : 
Soit le modèle 

Ye= ei + (8; 

Ur pury+v #1) 

do u=v-Aÿe Ur Vu -aycz (lé) 

vo he PU 

donc ver y. -aÿrie pv - avis) LT 

soit vi ye-fa+ D) vis +4 Piyer (IA) 

ty (at p}yui- (AP) ya + (9 
L'équation (19) traduit l'ensemble des équanuns 114) et (11. 
Dans {gi il v à seulement v, comme terme d'erreur qui satisfait toutes les 
propriétés d'un terme d'erreur. 
—  Estirrons l'équation f15) directement par Les MCO : 
LE Vr- Pr 
Er 


= Estimons {15} ça remplaçant y, prar sa valeur dans (19): 


est l'estimaleur des tratndres carrés ordinaires de l'équation (LS). 


Ÿ 


DIE +03} pe 2 + v, ls. 
DER 
_ (a +P)2 y à PrapE, ar. 1 +Zvire : 
LS 
rer F 
Zhirtri  2vive- 
LV Ne 


ä= 


=(a+p)-aÿ (20 


Donc, on peut déjà constater que à % il. 


Le calcul de fa fimite en probabilité de à, Plim â, permet de déterminer la valeur 
vers laquelle à tend lorsque Le nornbre d'observations, T, augmente indéfiniment. 
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En décalant (15) d'une période 
Vi = r-2 + 4,2] 
on obtient comme valeur de l'estimateur à des MCO 


qui nc diffère que par un {èrme du rapport 
2 Yr-17r-2 

Ze: 
Quand T — + cêtté différence devient négligruble, 
d'où 
Eve H=2 = plimé 

Yr-t 

Quant au dernier terme, son numératwur tend vers 0 puisque vet y,_, sont deux 
variables mdépendantes. 
Par conséquent, la limite en probabilité de à = a+9 -ap Plimä 
donc 
a+p 
1+ ap 


Plimé = 


Donc, même dans les meilleures conditions, l'estimateur des MCQ reste 
asymptomatiquement biaisé, On vérifie que & p = 0, Plim à = à. En l'absence 
d'autn-corrélation, l'estimateur des moindres carrés ordinaires est sans biais. 

—  Calçulons le sens du hiais : 

a+p-a(l+ap) 


Plimä-a = 
l+ap 
d'où 
pisse) 


l+ap 
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Limitons nous à [al < 1: 

- Si p > 0 alors Plim à > à — on surestime à 

_ Si p < 0 alors Plim à < a — on sous-estime a 

Le cas du modèle auta-régressié pur (c'est-à-dire sans vraies variables 
explicatives) correspond au biais maximum {C'est danc le plus mauvais). 

Quand l'équation comporte également de vraies variables explicatives (xt) le biais 
a tendance à être plus faible mais il existe toujours. 

Par conséquent, les résultats d'une telle estimation ne sont absulument pas fiables. 


Mais, facteur aggravant, quand on a ce type de figure (cndogène décalée plus 
auto-corrélalion) le test de DW est lui même biaisé. Si on a un problème 
d'estimation dans {15} alors logiquement on aura aussi un problème dans 
l'estimation de p dans {1}. En cffet Les deux équations sont tout à fail semblables. 


Or, le test de DW tend vers 2-2p cc qui signifie que si À est biaisé alors le 
résuitat du test est lui méme biaisé. 


En elfectuant le caicul de Ê on arive au résultal suivant : 


+ ap 


2 De, 2 
Pimp = ÉPXAP = P= EE 


Pimp = te +p) 


è 
Pimp = PET Ù 


Cas [al < L: 
- Si p<0 alors Plim f > p 


- Si p > 0 alors Plim À < 9 — onsousestime p 


< p dans tous les cas. 
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À est donc hbiaisé vers la valeur ( et donc le test de D'W est biaisé vers sa 
valeur centrale qui est 2. 

Duns tous les cas de figure, de test de D'W à tendance à miniariser l'auto- 
vorréliution {it est opticnistc}. 

Ccla amène à remplacer le D'W par un autre test en présence de variables 
cndogènes décalées. C'est le test h de DURBIX défini par : 


WC | 


RE NE Er 
PT Fvartdi) 


8, étant le cocfficient de régression de la variable endogène décalée. 

La valeur du test de Durbin peut être approchée par : 
DW\T OT 
(+ 2 Nr J'var(à) 

Remarque : pour pouvoir calculer la valeur de ce test, il faut que : 
T variance (411 < 1 
La variable h suit une loi normale centrée réduite 
hNio.l) 
à partir de laquelle on peul fire un test unilatéral 

HO :h=0 

H1:h>0 où h<0 selon le signe de p 
Seuil critique ho = 1,65 ali seuil de 5 % 

la] < 1,65 — H£ÿ donc pas d'auto-corrélation 

(4l> 1,65 — HI donc auto-corrélation 


Si p est élevé (forte auto-corrélation) alors on aura une valeur élevéé pour h. Les 
fortes valeurs de h sont généralement associées au rejet de l'indépendance. 
Exemple : un ajustement destiné à demner une première idée de la fonction de 
consommation. 
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Tunisie (1980 - 1990, données annuelles} 


ce, = 0,939, _1+0,3914, 21,296 {21) 
€0,099 (0.082) (880) 


R? = 0,989 

DW = 1.92 

Les tests sont bons et Je test de D'W semble bon. 

Test du b de DW : 

h= (Er) ue 
V1-TV(;t 1-31+0.095" 

h = 0,26 

h<hü = 1,65 — HO Pas d'ato-corrélation. 

Dans {6 cas où T c> L Durbin propose deux étapes : 

a- Calcul de à, 

b- Régresser ü, sur les variables explicatives et ü; 

Dans notre exemple. ü, = ac, + a2ÿ,+a3 + püt.] (22) 


A 


On teste ensuite +, {T 1-1). Dans l'exemple T - 4 soit 27 degrés de liberté. 
CP 


Remarque : on pourrait sur l'équation {22} faire le test de Breush-Godfrey et 
utiliser Ja table du y? à un degré de liberté. 
C- Les procédures de correction 


Il existe plusieurs méthodes qui ont cormne principe comrmun de passer par une 
élimination préalable du terme d'erreur ut. 


M = 4x, tés +u, (23) 
= Puy yt+v, (l) 


Cette élimination peut se fane à partir de il} : 
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Hg — Jr 4e — 62 

By — Pur] —Vr 
Fr Pin leu = Pr) pre%, 124) 
Je Prin a Por) tal ptér, (25) 
Dans (25) il ne reste que *, qui est un boñ terme d'erreur, présentant les propriétés 
requises. 
e La méthode de D'urbin (2 étapes) 
La première étape consiste à éliminer u, de l'équacion et à estimer (26) de façon à 
né conserver que v, comme terme d'erreur. 
L'opération conduit à : 

Ye = PS + diPtun +e2(—p)+v (26) 
On estime celte équation ce qui fournit des valeurs pour À,4,41p.é. 
Mais les cstimateurs sont 116$ Gnére eux par une contrainte qui est que le 
coefficient de xt, 4 est égal au produit des coc{ficients de ve, et de x,, Or. dans la 
procédure d'estimation rien ne permet de tenir compte de cuite contraint, Les 


éstimateurs obtenus né pourront donc pus être retenus. Is cunslituent un ensemble 
d'estimateurs non cohérents. 


On va alors en retenir un seul : À el connssant À il sera possible de calculer 4) 
Et &: 


Soient les nouvelle variables : 


M =Y- un ON 


rt 


* 


x =x px (28) 


y =@x; +a)l-Bi+v, (29) 
La 26 étape consiste à régresser (29) par les MCO (on à un bon terme d'errour} cc 
qui fournit les valeurs de À] et 42. 
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e Méthode de Cochrane Orcutt 
C'est une méthode itérative. Elle consiste à partir de l'équation (25} en se donnant 
une valeur arbitraire de 9. et à en déduire des valeurs de a; et ax. 


Connaissant ces valeurs de 3, et a, où revignc sur une nouvelle valeur de p,utc. 
Soit PQ la valeur d'amorçage de la procédure : 


Connaïssant 9 4 on forme ver <? ï 


RS (30) 
3 SX Pot (AD 
Puis on régresse 


y = ét x +a(l-p)+v, et on obtient des valeurs de a, et a, soient a? €1 aÿ 


Connaissant 4 et a2 dans l'équation (1} on peut calculer Ÿ, + afx, +a? 


On peut en déduire à, =}, -}, ét on obtient fP, à partir de (|), puisque 
p= Eü à 

san 

Zà 


Connissant p, on récornmence la procédure et on calcule 


=» (32 

xl=a, px (33) 
et dans (29) on retrouve aylet «3! 
Et dans {1} on en déduit fÿ; etc. 


Sargan a montré que cotie procédure cst convergente c'est-à-dire qu'au bout d'un 
certain nombre d'itérations les valeurs de p, a et a, vont se stabiliser. 


Il existe deux 1ypes de programme possibles : 


= On fixe à priori le nombre d'itérations 
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— On fixe un critère de divergence sur D, par exemple 1 %% c'est-à-dire que dès 
que deux valeurs successives de p diffèrent de moins de 1 & on considère 
que la procédure est stabilisée. 

En général la convergence est assez rapide, et, elle sera d'autant plus rapide qu'on 

aura bien choisi la valeur d'amorçuge. Comme le DW est en rapport avec la valeur 


1W 
de p il est commode de prendre Py =1- LT puisque DW—1-1p. 


+ Méthode de Hildreth et Liu | Méthode d'estimation par balayage} 
Elle consiste à estimer une succession d'éguation (25; pour des valcurs 
régulièrement croissant au décroissante de 2. 
Exemple 
p =0,1:0,2:0.3:...:0.9 
On construit pour cela de nouvelles variables sur le modèle de y} et # 
À chacune de ces régressions est associée la somme des carrés des résidus et on 
peut déterminer quelle cat la valeur de p qui minimise À (2. Puis on va repartir 
avce un pas plus fin. 
Ex. : Si la variance minimale est obtenue pour 6 = éou p = 0.7. 
P = 0,61: ; 0.69 
Et ceci jusqu'à abtenis le nombre de décimales souhaité (4 décimales cn général}. 


Les progrannnes de régression utilisent frégaemment la méthode de Cochran- 
Greutt. 


Exemple : 

L'exemple suivant concerne Les importations italiennes en données (rimestrielles 
de 1978,1 à 19924. 

Un premier ajustement a donné : 


LIM, =1,929LPIA,—0,140LP+0,572D90-12.067 (34) 
(4,28) {7.58 (2,313  (-15,68) 


R° =0,97 DW=1,14 Li 20,047 
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LIM, est Le logarihme des importations en volume, L PIB est ic logarithme du 
PIB réel, LP, est le Jogarithme de l'indice du défiatcur des importations de 
prdus munulacturés et D96 est une dummy pour lé 1% trimestre le 1990. 
Las valeurs eriiques du test de Durbin-Watson saut : 

di = 1,48 du = 1.49 
On est donc dans un vas d'aute-cortélation positive. 
L'ostimation par la méthode de Cochrane Oreut( donne 
FIM, »], ISTPIR D Es +0, ue [1,899 (3$} 

MCE 92) {12 #]: ( 16, 14) 

R° - 0,97 D 2,0 &=00S  D-0,417 
Convergence assurte après + ItÉrattons. 
Dans certains cas. il peut arriver que la procédure ne suffise pas à corriger J'aura 
corrélation. Cela peut traduire des problèmes de fond concernant les spécifications 
retenues. À côté d'une reflexion sur la firme de l'équarion, il èst alors utile 
d'explorer d'autres farimes possibles d'auto-corrélacion en particulior d'ardre 
supérieur à ]. Des corrections, génécalisant La méthode de Cochrane Orcurt à des 
processus d'ordre supérieur, éxislent dans Ja plupurt des louiciels et permettent de 
résoudre ce type de difficullé. Le 1esl de Breusch-Godtrey fournit des indications. 


SECTION 4 : L'HÉTÉROSCÉDASTICITÉ ET LA MÉTHODE DES 
MOINDRES CARRES GENERALISES 
À. Définition 
Par défininon, il y a hétéruscédasticité lorsque la variance du terme d'erreur { oc?) 
n'est pas constante sur l'ensemble des observanions (au cours de temps où sur un 
échantillon}. 
S'il n'y a pas en plus auto-corrélation, la matrice des variances-covariances des 
termes d'erreurs E (UU”} reste une matrice diagonale 
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2 
oc? 0 
E(GU = , (D 
û oë, 


Mais Les éléments diagonaux nc sont plus égaux entre eux. 

2 ï: RE. : : ; 
Cette matrice va s'écriré Ojr,7,V # Of ve qui entraîne La remise en cause des 
hypothèses 1 et 4. 


Lorsque l'hétéroscédasticité correspond à la présence d'un élément systématique 
déterminable dans le terme d'erreur, fl est possible de la corriger après l'avoir 
identifié. 

Mais l'hétéroscédasticité peut-être purement aléatoire, ou dûe à des causes 
inconnues et dans ce cas il est pratiquement impossible d'y remédier. 
L'hétéruscédasticilé est un problème assez fréquent qui peut être rencontré dans 
un bon nombre de situarions. 

Exemple : on s'intéresse au taux de rentabilité d'un ensemble d'entreprises de 
tailles différentes. On sac par ailleurs que les grandes entrenrises opt un taux de 
rentabilité plus stable que les petites, mais pas nécessairement plns élevé. 

Uue représentation graphique du taux de profit en fonction de l'actif total de La 
fine fera apparaître La forme caractéristique du nuage de points : 
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Actif 


Supposons maintenant que l'on cherche à évaluer l'influence de la concentration 
du marché et de la taille de La firme sur sa rentabilité, ce qui à constitué un thème 
de récherche important dans les années 60 et 70. 


On cherche à estimer l'équation : 
li =d +4 $ +a3 À +u, (2) 


ox de pro per de poleur de: 
MErChé arrife 


Dans une telle équation on aura de l'hétéroscédasticité. Le résultat sera différenc 
selon qu'on tient compte ou non de la présence d'hétéroscédasticité. De fair, Haki 
et Weiss en 1967 ont montré que les résultats des études sur les effets de la 
concentration étaient subsianticlièment modifiés lorsqu'on tenait compte de 
l'hétérascédasticiné!t}. 
B- Comment faire apparaître la présence d'hétéroscédasticité ? 

Il existe différents tests permettant de faire apparaître l'existence éventuelle de 
l'hétéruscédasticité : 


al et Weiss ‘Firm size and profitabiléts”, The review uf Economics and statistics, Août 
FD 
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— Le test de GOLDFELD & QUANDT 

Le test de GLEJIFR 

+ Le test de Gotdfeld & Quand 

Î purt de l'idée que si il y a hétéroscédasticité ayant un caractère systématique, 
c'est aux deux extrémités du nuage qu'elle sera la plus visible. 


[l faut donc exclusivement s'intéresser aux extrémités du nuage. Pour cela on 
supprime D observations situées au centre dé l'échantillon. 


: , —D - ; 
Si on a T observations : on ne garde que observations aux deux extrémités. 


Sur chacune de ces extrémités on va effectuer la régression et on retiendra »aà 
et Zä, puis on cherche à vérifier si la variance résiduelle est la même à partir de 
ces deux régressians. 

Pour cela on procède à un test d'analyse de la variance qui s'appuic sur le rapport 
des deux variances, en plaçant au numératcur la plus forte variance. 
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Exemple : on suppose que la variance résiduelle est plus forte dans la première 
sous-période. 


Eu, > Lui, 


Yu, 
T=2 
on à 2. > Fisher Sned [TP T-RX) (3) 
Zu, 2 2 
TD, 
2 


?= nombre de variables explicatives 
Soit For le seuil critique de la table de Fisher Snedecor, on leste 
HG :hamoscédasticité 
H1 :hétéroscédasticité 
La règle de décision étant ta suivante : 
*F<Fa - HD 
*F23Fx —HIi 
Si les variances sont vraiment différentes, le rapport sera élevé et donc plus Ja 
valeur de ce rapport augmente, plus on a de chance de conclure H]. 
Remarque : Le lest de Goidfeld et Quandt se généralise aisément à 2 sous- 
périodes de durée quelconque T; et T3. Dans ce cas : 
Eä? 
T, - 
Zi} 
BB! 
e Le test de Glejer 


11 consiste à régresser la valeur absolue ou le carré des résidus sur l'ensemble des 
variables explicatives, 


— FT -{7 1) (4) 


‘14 e Z<s 2 
On utilise généralement le carré du résidu à 


régression 


qui est facile à obtenir après une 


ei] 
on 
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à lou; = A9 + Xi tt 
On s'intéresse au test de Student sur chacun des vocfficients de régression. 
Si le test de Student classique conelli que tous les cuelficients de régression sont 
non sigmuficatifs. aj = QG = 1. [ on conclut à l'homoscédasticité. 

Si ai # O c'est-à-dire qu'il y &« au moins un coefficient significatif on conclura à 
l'hétéroscédasticité er, en même temps on en déduit qu'il y a une des variables 
explicatives qui explique l'ardre de grandeur des termes d'erreur. 


Autrement dit le test de GLEJER donne deux informations : 
— La présence d'hétéroscédasticité 
— La cause de cette hétéroscédasticité 
C- Correction de l’hétéroscédasticité et méthode des Moindres Carrés 
généralisés 
L'hétéroscédasticié né concerne que La maurice des vartances-covariances des 
termes d'erreur ct n'affucte pas les eslimateurs des MCO qui demeurent sans biais, 
Par contre. les variances de ces mêmes estimations et danc les tests de Student 
seront hiaisés. Il st donc nécessaire de corriger l'hélérosvédasticité. 


Pour cela an peut utiliser une méthode proposée par AJTKEN dans les années 26 
connue sous le nom de méthode des Moindres Carrés Généralisés qui permet 
d'obtenir des estimateurs sans biais V'la matrice E(UU'j. 

Supposons donc que Æ(UE")= o°V. V étant une matrice définie positive c'est-à- 
dire que pour tout vecteur Z::7, alors Z'VZ > D (Z'VZ est un scalaire). 


Cette propriété jouc ici le même rlt que la condition qu'une variable soil posttive 
pour adinetlre des racines carrées réelles. 


On peut metire V sous la forme du produit de deux matrices M et M': V = MM 
de forrnat (T,T} M étant une matrice non singulière admetçant une inverse M° L. 
V est une matrice symétrique V = V', par construction. 

De même, si V = MM' alors -! = M°! M: 
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En rultipliant à droite et à gauche on retrouve I: 
M VM' EI 
Il est possihle d'utiliser ccite propriété pour comger l'hétéroscédasticité : 
On pare du modèle: Y = XA +0 (5) 
Soient les nouvelles variables construites en prémultipliant par M°! : 
Y*=MY (6) 
x*= MIX (7) 
U*= MU (8) 
Sur le nouveau terme d'erreur défini on a : 
ÉtU*L'#} = ECM! UU'M"} (9} 
= M EQUU') M"! (10) 
= Mo? vm ! = 9°} 
+ Les nouveaux termes d'erreurs satisfont donc toulcs les propriétés des MCO. 
Vérifions que l'estimateur À des MCO sur les données transformées est sans 
biais : 
À = exe) Xe vx 
ex M Ou x x l'OM ta : MU) QD 
Soit À = A+{X' VX} X'VIU (12) 
donc E (À) = A si H(LN) = 0 
L'estimarcur est sans biais. 


Qusat à la variance de l'estimateur sur l'équation transformée, 
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ViÂ)= rà- raz» | 
x vixi ir vtt euu ve Uxex VIE 3 


s GX VX) NA VOEUX CAN VTIXTT (14) 


HIT Vi QUE y 
=o(X VX) ! 
= ŒUXS XX) car  (X*= MIX) (IS) 
(XVUX=- M IMTIX=XT X) 
- La variance de À est donc donnée pur les MCO sur les variables transformées 
Le problème cst cependant de connaître M. Cela suppose d'avoir des informations 
sur [a matrice V. Mais dans la piupart dos cas ces informations nc sont pas 
disponibles. C’est pourquoi la méthode des Moindres Carrés généralisée n'est 
applicable que dans un nembre limité de cas. 
C'est ce qui fait l'intérêt du test de GLEJER puisque ce test peut fournir 
l'information dont on à besuint. 
Si ce cost indique par exemple que la variance du taux de profit est inversement 
proportionnelle aux actifs, dans l'équation (2j, un peut corriger l'héléroscédasticité 
en utilisant cetre information. 


V= . (16) 
L ‘ L/Ax 
V est donc assez bicn spécifiée. 
On peut mettre V sous la forme Ÿ = MM' d'où : 
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14 A 
} L] 
M “| . 7) 
+ 
| l'VAr 
Va 
M°l= . {18} 
VA & 


Donc, l'application des Moindres Carrés généralisés revient à multiplier toutes les 
- Ds - 

observalions pur : À; 

Et, quand les variables sont transformées à faire La régression sur les variables 

transformées. 

Remarque : La méthode des Moindres Carrés Généralisés (MCG) est aussi 

applicable en ças d'auto-corrélation dés termes d'erreurs puisque Ià aussi on 

contaît la nature de la matrice V. 


1] pes p}1 
£ + L 
val las 
Fc .. il 


Cette matrice admet un inverse V'l 


Me: 1+p° -p 0 0 0 
0 -p 1+p? -p 0 TE 


1 
EU O0 0 0 -pil 


{20} 


Qui peut être mis sous la forme V7] = M") Ml 


L'estimation d'une équation unique 59 


Fe @ 0 0 


1 —f 1 
DVec MU I 0  -pi ol 1) 
F7 
| () 0 -p 1 
L'application des MCG conduit à régresser sur les variables transformées. 
LT ni 
X*=M TX 


Compte tenu de la forme de Ml, appliquer ies Moindres Carrés généralisés 
revient à Lrès peu de choses près à l'applicauen de l'estimation sur l'équation {25) 
de la sectiun 3. 

Dans {a plupart des cas où la matrice des variances-covariances est inconnuc on 
ne peut pas appliquer les Moindres Carrés généralisés. C'est pour cela qu'on 
distingue les Maindres Carrés généralisés &t les Moindres Carrés généralisés 
“aisables". 

Exemple : {on laisse de côté la question du la pertinence du modèle qui séra 
abordée au chapitre 4). 

Une régression simple du logarithme de la masse monélaire (M2) sur le log du 
PIB, le taux de chômage, le taux d'inflation et Le coût d'appornmité sur 1986.1- 
1993.4 donne : 


LM; =] Pre Oc#-- 9,89 nf 0,025C0 0,928 (22) 
201 (2,45) (4H) (ëdi {0,65 


R° =0,92 DW = 0.46 


Après correction de l'autu-Corrélation, le résultat obtenu est aherrant avec en 
particulier une élasticité revenu négative. 


Le test de Glejer sur cette régression donne : 


&? = ©, 0006 LPIB— 8.4 107 ch 0,011inf—1,710 C'O+0,008 
1.94 {—0,48) (1.50) 


R? =0,13 
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Après transformation des variables ea utilisant {eh)2$ on obtient : 


LM: = 1, De 022ch—1{0,74inf+0,022C0- Q, Se (24} 
69) (1.75) 1 4.91) (LA) (01 


R? = 0,99 DW = 0,48 


SECTION 5: PROBLÈMES LIÉS À L'INSTABILITÉ D'UNE 
RELATION AU COURS DU TEMPS OÙ À L'HÉTÉROGÉNÉITÉ 
D'UNE POPULATION 


Il arrive souvent qu'on, travaille sur des échantillons non homogènes. (Exemple : 
population de salariés où il y a des hommes et des femmes}. De même on peut 
imaginér que certaines relations ne sont pas stables au cours du temps ef qu'elles 
se déplacent (Exemple : quand le cadre institutionnel change, certaines relations 
se modifient). 

Dans certains cas il est possible de séparer la période ou l'échantillon en deux 
parties et de travailler séparément. Ce n'est pas toujours possible quand on n'a pas 
assez d'observations Gu quand on s'intéresse justement au déplacement. D'où la 
nécessité d'estimér sur l'ensemble de la période sachant qu'il y a un problème 
d'hétérogénéité ou de stabilité. Dans tous Îes cas, le fait d'ignorer ces phénomènes 
peut conduire à des absurdités. 


Exemple : la courbe de Phillipps sur période longue. 


w taux de croissance des salaires 
t taux de chômage 


Sur longue période, les conditions de l'arbitrage inflation-chômage tendent à se 
dégrader. 


Sur chacun des cycles ainsi décrit on peut ajuster des courbes de Fhillipps qui 
auront une élasticité négative mais la courbe d'ajustement sur l‘ensenible des 
points peut très bien avoir une pente positive. 

Donc l'estimation sur l'ensemble des observations peut aboutir à des absurdités. I] 
faut alors vérifier qu’il y a bien stabilité ou homogénéité. À défaut, il faut adopter 
un modèle un peu différent pour en tenir compte et travailler dans des conditions 
convenabies. 


On peut tester la pertinence de l'opération en utilisant des tests de restriction. 
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Pour cela, on considère que : 
— La stabilité est un cas particulier de l'instabilité ; 
— L'homogénétté esç un cas particulier de l'hétéragénéité. 
À- Comment s'assurer qu'il y a homogénéité ou stabilité ? 
ILexiste pour cela plusieurs méthodes. 


+ Premier test 


Yt 


AT) 
CEE 


1'E sous-période 


2*sous-périodc 


Une prenuère méthode est de chercher si les ajustements effectués sur les deux 
nuages sont diftérenis. Si le nombre de point est suffisant on peut régresser 
séparément sur les deux suus-populations où sous-périodes et à partir de là tester 
l'égalité des coefficients de régréssion. 

Muis ce test est peu performant car cela revient à construire des intervalles de 
confiance à + 2 écarts-types et à regarder S'ils ont une partie commune. 


e Deuxième test 

On peut régresser sur l'une des sous-population et se demander si cette régression 
est aussi bonne pour la deuxième sous-population. On calcule des écarts résiduels 
sur les deux sous-populations et on cherche à savoir si La variance résiduelle de La 
deuxième sous-populatinn est aussi bonne que celle de la sous-popalation 1 (test 
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d'iomogénéité). Sinon il # a hétéroscédasticité : La variance résiduelle est plus 
forte sur la 2° sous-popularion que sur la première. 
e Troisième test 
Le test te plus performant est le test de Chow. Supposons qu'on ait les deux sous- 
périudes d'observation ayant T\ et T2 obsersations respectivement et variables 
explicatives. Ce test demande qu'on fixe à prion la date à partir de laquelle la 
relation change où qu'on définisse ce qu'il y à exactement dans les deux strus- 
populations. 
À partir de Jà le principe du test est de voir dans quelle mesure le fait de régresser 
séparément sur les deux sous-péniodes améliuee Le résultat de la régression. Ce test 
portera sur {es variances résiduclics. 
On commence par régresser sur La population totale soit T observations. 
LETs+T; 
Puis on refait la même régression sur chaconc des 2 sous-popularions séparément 
ec on retient la sormme des carré des résidus de chague régression. 
Soit : 
sur T, — Xäf 
sur 7 > Ÿ 4, 
Eü +34? doit être < X 67 
mais dans des proportions qui dépendent de l'amélioration apportée par la prise en 
compie de l'hétérogénéité des données. 
On va calculer 
A2 2 1 + 4 
Zä (Zi, -E&,) 
Î 

Nombre de degré de liberté du aumérateur = 7 -{—{7 —-127 —-1)=1 
La, +25, 
A +75 -2l 


11) 


que l'on rapporte à 
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247 -(Séf, +24,) 


soil — #1 — =F 
DER + Lé, 
yo +1 — 21 


Ce rapport suit une loi de Fisher Snedecor avec respectivement / degrés de liberté 
pour Je numérateur ec {T - 2 ]} degrés de liherté pour le dénominateur. L'ordre 
numérateur / dénominateur est imposé par le test. 
Interprétation 
Au numérateur on trouve l'amélioration apportée par la prise cn compte de 
l'hétérogénéité et au dénorminaleur la plus faible des variances résiduciles. 
Soit F@ Le seuil critique de la table, et F la valeur calculée par Le test. On pose : 

HG = homogénéité ou stabilité : 

HI = hétérogénéité ou instabilité. 
Si la relation est très instable le gain obtenu est important ct on aura unc valeur 
élevée au numérateur. 
Donc si F > Fu on conclura HN] 

sil < Ft on conclura H@ 

Ce qui fait l'intérêt du test de Chow c'est qu'il est sensible car, au numérateur il 
porte à la marge sur des variations et non sur les variances elles-mêmes. 
Dans cette version, il faut cependant déresminer la partition de ja population 
totale. 
Quand il s'agit de phénomènes chronologiques, ce sont des événements extérieurs 
qui vont donner des indications sur la période charnière. Souvent la période 


chamière sera celle qui maximisera le test de Chow. Certains logicicls permettent 
de calculer une succession de tests de Chow pour différentes périodes, 


Exemple : la fonction d'importation de l’Iralie éiudiée dans la section 3 est-elle 
stable ? 


En observant le graphique des résidus, on à remarqué un changement dans la 
qualicé de l'ajustement avant et après le 4* trimestre 1987. 
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Ta régression des importations sur le PIB, les prix et la constante est refaite pour 
les 2 sous-périodes 1978.1-1987.4 et 1988.1-1992.4, 


Les sommes des carrés des résidus sont respectivement : 
Ea? = 0,0329 


Zi? = 0,0423 


L'ajustement sur la période entière donne +0 =0,1774 (1 
Le test de Chow vaut donc : 
0.1774—{0,0329 + 0,0423} 


3 
(0,012+0,047 #4 
54 
La valeur critique Fa (3,54) = 2,78. 
La fonction est donc instable, 
B-La prise en compte d'un éventuel changement 
+ Premier cas : la foncliun se déplace et garde la même forme. 
Yt 


XI 
Pr 
Ti T2 


tir] s'agit du modèle sans la variable muette D 9, dont la variance résiduelle est nettement 
supérieure à celle de l'équation présentée à la section 3. 
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Exemple 
Soit une fonction d'importation de la farine. 
Login, =a logr, +aslagP,+az (li 
On peut rendre compte de ce déplacement en introduisant une variable muette qui 


servira à caractériser l'une des deux périodes. Tout dépend de ce qu'on choisit 
comme situation de référence. 


Si on considère que la référence est la première période, on caractériscra là 
seconde pénode : 

L=0Dsit<T; y 

Lz0treT, 

L=lt>T ou 

L=lste D 
Ce qui conduit à ajouter une quatrième variable explicative : L dans la régression 
(1) 

Log int = a log, +u: log, +43+a4L 2} 


Log & Log PB 110 


2 : 


RES PRES Et à 
AE 
[l 
[l 
° ll 
Log à 16 
Interprétation des coefficients 
a, mesure le déplacement de La fonclion d'importation qui peut être engendré par 
un événement èn T;, par exemple, un changement institutionnel comme la 
libéralisation du cornmerce extérieur. 
Remarque : 


Si on avail estimé l'ensemble par une relation unique tous les coefficients auraient 
été faux. L'introduction de L rend cempte d'un saut unique sur les importations. 
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+ Second cas : La fonction se déplace et se déforme 
Yt 


Xt 


*+——+ <———+ 


T1 T2 


La pente de l1 fonction est donc différente en T; de sa valeur en T, 
Soit la fonction Log [ms = a4 log y, + aa log Pei+us (1) 


On peut introduire sut Le coeflicient a; un effet propre à la deuxième période ce 
qui va nous conduirt à estimer la régression suivante : 


Eog Lmmç = {a + ag L) log yc + as + log Pau +asL (5) 
ag va mesurcor la variation de l'élasticité revenu des importations qui est ici 
imputable à la seconde période. 
L'équation modifiée va tendre compte de ce relèvement d'élasticité, 
On désigne parfois ce type de fonction sous le nom de modèle “spline" 
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LF 
Log M 1 0 Q. 
° f 
XGr4 = 7 
f 
T 


Log Y7 L } Log Xr 
Pour savoir si ce modèle est performant, il existe deux méthodes : 
— voir si le modèle cost validé c'est-à-dire si a4 est significativement # ( : 


. — Student 


Gég 


Et de même pour as 


— Utiliser un test de restriction à priori. En effet on peut toujours considérer le 
modèle général comme une version appauvric du modèle particulier (avec a4 
el às). Quand on passe du modèle “Spline" au modèle général on diminue là 
qualité de l'ajustement. Si la perte en pouvoir explicatif est importante alors le 
modèle Spline est meitleur que le modèle général c'est-à-dire qu'il y a bien 
deux périndes séparées. Sinon, si le modèle général est aussi bon que 
modèle Spline, la fonction est homogène ou à peu près stable, 


Le modèle sans déplacement est un cas particulier. Cela danne la possibilité de 
tester la validité de l'hypothèse faite sur le déplacement. 


On revient par ce moyen au test de restriction à priori. 
Soit H = matrice des restrictions (q,1) (autant de lignes que de restrictions. 
autant Je colonnes que de coefficients à estimer) 
HA =h (6) 
La contrainte peut être que : 
— Le terme constant des deux régressions soît le même, ce qui implique que la 
valeur du coefficient de régression de la variable muëtte soit nul. 
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Sait l'équation 
ex, ta +aL (7) 
Onimpose 23 =0 


a | 
Â= | H=(001) k=0 


43 


fi] 

(6 01INa; |=0 (8) 
43 

Le modèle dans lequel on impose une restriction est moins bon que [c modèle sans 

restriction, Ce modèle aura donc une Z à? > à celui du modèle sans restriction. 


La différence entre la somme des carrés des résidus avec restriction et sans 
restriction nous renseignera sur le caractère acceptable où non de la restriction. 
Cela servira à effectuer un test sur la validité des restrictions posées : 
HO HA =h c'est-à-dire que les estimaieurs peuvent respecter la 
contrainte. 
HI HAzxh c'est-à-dire que le modèle est non compatible avec la 
contrainte retenue, 


“2 
Soit Zw, la somme des carrés des résidus du modèle avec restriction et Z 2 la 


somine des catrés des résidus du modèle sans restriction. 
On à alors 


È w2 TT Ÿ «2 
— 2 — + Hig,T-1) 
Zé; 
T-i 
Ce rapport mesure dans quelle proportion la restriction dégrade Ie modèle, 
Soi Fix, la valeur critique de la table de Fisher, Si F < Fc alors l'imposition d'une 


contrainte ne dégrade pas beaucoup la qualité de la régression et on conclut denc 
Hg. 
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$SiF > Fa alors la contrainte diminue de façon importante la qualité du résultat et 
on conclura alors H1, 


Remarque : Comme on l'a vu, les tests de restriction à priori sent d'utilisation 
très variée. 


[ls permettent de tester toutes sortés de contraintes de forme linéaire du type 
dia + &3a, = Z telles qué rendements constants, homogénéité de degré 1 d'une 
fonction demande par rapport aux prix etc. 

Néanmoins une utilisation intéressante des Lests de restriction à prioni st de secvir 
à l'analyse de variantes d'une même régression. Par tapport 4 un mudèle de basc 
cela servira à évaluer cc qu'apporte une spécification moins générale, plus précise. 
Cela peut aussi servir à comparer des variantes d'un même modèle explicatif, 
comme par Exemple : 


. Variante avec retatd / variante sans retard ; 

. variamie avec unc struciure de retard / variante avéc une autre structure. 
Tous ces problèmes qui relèvent du choix de l'économètre peuvent être éclairés 
pat des tests de restriction à priori. 


SECTION $ : LES RETARDS D'AJUSTEMENT 
Les variables économiques ne réagissent pas, le plus souvent. instantanément mais 
avec un certain délai d'ajustement. 


11 faut un délai de perception, puis le temps de faire des choix, de les mettre en 
application. Tout cela ne peut pas se faire de façon instantanée. 


En conséquence, l'explication instantanée d'une variable par une ou plusieurs 
autres ne rend compte que d'une partie de la réalité et il faut distinguer les 
réactions de court terme et les cffets de Jong terme. 


Avec un modèle du type 
Ye = Xe + AA, 3 + dat 2 +. +aX , +u, (1) 


. âÿ 
— l'effet de court terme est 3 CT 
X, 
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l'effet de Long terme est Fa. 


En effet, si on applique un acéroissement dé x constant sur toute la pénode de 
référence de i = 0 ài = non retrouve Av = Za,Ax 

Un tel modèle est dit à retards Échelonnés. C'est donç un modèle dans [lequel y 
dépend des valeurs passées de x. On peut aussi avoir : 

Ve = Vote On Von + PiXie bntos +8, (2) 

C'est un modèle auto-régressif, dans lequel lu variable expliquée dépend de ses 
propres valeurs passées et le cas échéant des valeurs présentes ct passées d'autres 
variables explivatives. 

1 est commode d'introduire un opérateur qui est l'opérateur de retard désigné par 
L te] que 


ie = À. 
Fe — 
Lx=x,; 


Cet opérateur possède les propriétés mathématiques habituclles. 

Il peut être utilisé pour réécrire les modèles avec retard. 

Ainsi (1) peut s'écrire : 

ÿe = ACL) x; c'est-à-dire une suite de retards échelonnés sur la variable x. 

La structure générale avec tous Les types de retard pourra s'écnire : 

Ye GYs-le ne +byxt, br +. +bpx-4 +k, (3) 

soit AÛE)y, = BL)X, +u, (bis) 

Un modéle du type (34 ne peut pas dans le cas général être estirué directement, If y 
a, on général, différents problèmes économérriques qui rendent l'opérauion 
impossible. Au lieu d'avair une structure simple avec u, on aura souvent une 
structure de retard cumplexe sur u, et l'estimation par les MCQ sera impossihle, 


Même ta forme {1} fera apparaître une collinéarité des variables explicatives, dans 
un bon nombre de situations. 
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Comme on ne peut pas estimer directement les formes retardées, on recourt à des 
distributions de retards qui vont permettre une estimation à partir d'un modèle 
modifié : 


[existe deux types de structures de retards : 
les distributions à priori; 
-les retards polynomiaux d'Almon. (Polynomial distributed Tags. PDL en 
anglais}. 
À- Les distributions à priori 
Ce sont des distributions que l'on se donne et qui permettent de situer les 
coefficients rerardés les uns par rapporis aux autres. 


Soit le modèle v, = Za;x, ; +u, 
Plusieurs distributions à priori sont utilisables : 


La distribution rectangulaire 


On suppose que l'influence de la variable décalée est la même pendant k périodes 
puis s'estompc. 


8, = 8 pour i = 0, …., k 
a, = 0 pour i >k 


v=alx…. xx] + u (4) et on n'a donc qu'un seul coefficient à éstimer. 
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La distribution en V renversé 
Elle permet de traduire une influence croissante jusqu'à un maximum puis 
décroissante, de la variable explicative. 


Supposons que celte influence augmente pendant les deux premières périodes de 
retard, passe par un Maximum pour i = 2 el redescende ensuite. 


La spécification de cé processus est : 
HT aolx +2x +3x,-2 + 2x, 3 + x=4] +4 (5) 
Il y a ici également un seul coefficient à estimer. 


Cette formulation simpliste est assez bien adaptée au problème de 
l'investissement : l'influence des variables décalées {taux d'intérêr. profit) est plus 
importante que celle de la période contemporaine. 


Remarque : Le modèle pourrait aussi être cunsidéré comme une version contrainte 
d'un modèle plus général ct on pourra en (ester lu validité. 
La distribution à poids géométriquement décroissanis 
a = À avec O« À <1Les coefficients retardés sont fonction des puissances 
successives d'un coefficient À <1 
Exenple : ÿ,= ax, + ares +aÂx,2.. aAxs, (6) 
Un tel modèle peut être transformé en modèle auto-régressif par l'intermédiaire de 
la transformation de KOYCK. 
En décalant d'une période l'équation (6) et en La multipliant par À, il vient : 
kyey=a ke +aMxip+.satix ss (7) 
En soustrayant membre à membre, i] vient 
Ye AY = ax al) 18) 
Soil encore 
Ye ax + À yes - a Atlxp ni (9). Le dernier terme tend vers ( quand n 
augmente. 
d'où 
Yes a+ Ayei (10) 
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La transformation permet donc de passer d'un modèle à retard échelonné à un 
modèle auta-régressié. 
Lau transformation de KOYCK pose un problème économétrique puisque si dans 
l'équation 4 les termes d'erreurs on! Les propriétés requises {3-4} i]s ne peuvent 
plus les avoirs dans (10). Dans (10) on a cffectivement un terme d'erreur composé 
&u- Aus} 
Il faut distinguer les effets de CF et de LT dans {6} 

A 
ACT ea 

Àx, 
* A LT: supposons que l'on applique une variation de x égale à Ax sur toute la 
période considérée. 
Xt devient x{;+Ax Vi 
En conséquence y, augmente de Ay, que l'on peut calculer : 

Ay,+r, aix, — Ax)+aÂ(x._] +Ax)+...aÂ"{x,_, + Ax) 


-|ar + aÂx,_ +...+ax, a F aix + up Ax+...+aÂ"Ax 


Le terme situé entre crochets n'est autre que y, lui-même. 


D'où il vient: Ay, = a[Ax+ 4x... 4" ac]soit en supposant que n est suffisament 
grand. 


1] y a donc bien 2 effets distincts 


Av 
etlet de court terme Lac 4 
% 
äy, mn 
effet de long terme —L =—— 
; Ax  1-À 


4 étant le coefficient de l'endogène décalée dans (10). 
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On voit que l'introduction de l'endogène décalée y, dans une fonction peut être 
interprétée comme la traduction d'en mécanisnx de retard avec des poids 
géométriquement décraissants. L'endogènc décalée est ie signe d'un retard 
d'ajustement, 


Soit le modèle y, = aix, +325, +uav_] (11) 

Pour faire apparaître ses propriétés de long terme on applique un Ax, permanent 
d'où : 

Yi + Av = a Üxis + Ax) + asxx + à2 Cyr. + y) puisque en régime permanent y; 
augmente de Ay 

Yr+Ày=axir + 42a2t + ay +A]AxI+ a;Ay (12} 

d'où y =a, Àx] +a,4y 


4) 


Soit =; A (13} 


Un mécanisme voisin : l'ajustement partiel 
#,= valeur courante 


#4 = valeur de long terme (dans certain cas c'est la valeur d'éguitibee ou la valeur 
désirée). 


Entre les deux valeurs il ÿ a un mécanisme d'ajustement partiel qui peut s'écrirc : 
es = AY *—73, 1) avec D<%A<l (14) 
La valeur de long terme peut se définir par Une certaine fonction par exemple : 
VY= a +aix, 
Po Vrni © Aa +uxe)— y 
soit 

M = +aax, +(1-24)5 à 


Cette équation traduit l'ajustement du niveau effectif de la vatiable À ses 
déterminants, ii xt. 
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Si on cherche la réponse de long terme de y, à x, on trouve ici 42 : C'est épal au 
coefficient de régression x, sur un moins le coefficient de régression de 51.1. 
jette 

1-{1-À4) 


— On voit que deux choses assez différentes conduisent finalement à {a même 
spécification : 
. un phénomène de retard dans la réponse d'une variable à une autre 
{équation 6) ; 
. un phénomène d'ajustement partiel {équation 14). 


Le mécanisme d'ajustement parüel s'écrit souvent sous une forme à élasticité 
constante ; 


y +Ÿ 
[2] = [22] À = coefficient d'élasticité 
ei Jr-1 


On peut l'appliquer, par exemple, à l'emploi : dans un premier temps, en inversant 
une fonction de production, on détermine l'emploi nécessaire, 


Q:f{K.L) 
d'où L*=f!(Q, K}. 


Mais on sait que l'emploi effectif ne s’ajuste pas immédiatement à l'emploi désiré 
d'où l'appel à un mécanisme d'ajustement, Qn fera ainsi apparaîlre des vitesse 
d'ajustement selon Les branches, selon Jes pays... 


Remarque : si à chaque période on ajuste 20 %, le délai d'ajustement sera 


22 soit 4, dans ce cas. 


Toutes ces modalités de retard sont des retards à priori. On préférera généralement 
utiliser la méthode des retards polynomiaux d'Almen, 
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B- Les retards polynomiaux d’Almon 


Soit le modèle ; 
M = GX, ++. tax, +, (15) 
ou 


# 
sax. itu, 
EL 


Comme on ñ6 peut pas directement estimer (15) on utilise un attifice consistant à 
estimer une fonction polynomiale qui sera suffisamment proche des valeurs vraies 
des coefficients. Soit F cette fonction polynomtiale. On remplacera les vrais 
valeurs des coefficients par leur valeur dans a fonction. Îl faut donc que cette 
fonction soit une bonne approximation. 

Après avoir estimé la fonction on pourra déterminer ainsi n'importe quel 
cosfficient 4. Cette opération est une application du Théorème de Weirstrass qui 
établit que n'importe quelle fonction continue et bomée sur un intervalle peut être 
approchée par un polynôme de degré convenable. 

On suppose donc qu'il existe une vraie fonction inconnue qu'on peut approcher 
par une fonction continue. 


Une fonction polynomiale est unt fonction de type : 
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(16) a; = ag +@iit di. i" — polynôme dé degré r 

Pour connaître les pararnètre de F il faut connaître les ot; afin de pouvoir calculer 

n'importe quelle valeur de a. 

Cette opération n'a de sens et d'intérêt que si le nombre de coefficient & est 

inférieur au nombre de coefficient a car : 

— Ce serait illogique d'estimer plus de coefficients que dans la forme de départ. 

— Sachant que l'on à (n+1) coefficient a et (r+1} coefficient &, si r > n le 
polynôme n'est plus défini de façon unique, c'est-à-dire qu'on pourra faire 
passer plusieurs polynôme par l'ensemble des points qui ont été identifiés, 

La méthode consiste à estimer les coefficients & de la relation (16) pour ensuite 

calculer les a. Pour cela il existe plusieurs méthodes possihie. 

La méthode simplifiée d'Almon consiste à réécrire l'équation de départ en fonction 


des coefficients œ el à les estimer par les moindres carrés. 
Supposons que l'on ait : 
n=4 
Ê =2 
a=œ+@i+@i {17} 
Ye = dort 41Xe] + A2Xe2 + 43X a + A4Xt4 (18) 
ag = Yo 
a = Co + + 
a3 = y +20 +4@; 
a3 = Gp +34, +90 
da = Cp +40, +16 


Remplaçons les ai par leur valeur en fonction des &: 


Yu = Cox + CO, + Os + Ds) Xe + Lo + 20 + data) Xe + (Ge + Roy + Otis) Xe 2 + 
(es + 40 + 1603) xn4 (I9) 


L'estimation d'une équation unique 79 


En mettant en facteurs les & dans (19) : 
Yi = Mgr + Xe-i + Xe. 2 +2 + Xr 4) 
Hors +240 +32, 3 + 4x, 4) 
+d>Üxs 1 + 4x2 +9x,.a +164, 4) (201 
Dans chacune des parenthèses on retrouve de nouvelles variables soit Zo. Zi, Z2. 
On peut donc écrire : 
= doZote +22 (21) 


Lä méthode simplifiée consiste à régresser y sur Zp, Zi, Z2. pour ohtenir des 


estimaleurs @t,@i ef. 

Ces estimateurs servent pour calculer les valeurs de a. 

(Remarque : if existe une autre méthode de résolution, fatsant appel aux 
polynomes d'interpolation de Lagrange.) 

Exemple : on s'intéresse aux déterminants de l’etnploi en données trimestrielles 
dans Le secteur de l'énergie. Les premiers résultats font apparaître une relation 
négative entre l'emploi et I capital productif ainsi qu’une influence négative du 
taux de salaire. Considérant que la réponse de l'emploi aux variations du taux de 
salaire n'est pas instantanée, on introduit un polynome de retard sur cette variable. 
La durée totale des retards est de 8 trimestres et le degré retenu est deux. 

n=8 

r=2 

Sur la période 1906. L - 1994. IV le résultat obtenu est : 


8 


EN=-1330LK+ Ÿ a, LW_; +25,864 
(CAL 1 (L4173 
R? = 0,998 


DW=185 


=0,179 
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ag =0,07857 a =0,05812 4: = 0,040 
(4,40) {5.40} €7,$0? 
ax =-0,02601 as =-0,01435 a; = —0,00562 
7,92} 3.00} (9,901) 
=90,0x18 =6,00305 = 0,00298 
ds (9,026) Pi (0,532) "4 (0,847) 


Za; = -6,21706 
(8,83} 


Les variables sont en logarithme. Les chiffres entre parenthèses sont les t de 
Student. 


Dans cette procédure on peut contrêler Le choix du degré du polynôme. Il ÿ a deux 
paramètres importants : net r : 

* nest souvent dicté par la théorie économique ; 

* le degré du pelynôme : plus il est élevé, plus {es calculs sont lourds, plus on 
augmente le risque de collinéarité dans ta méthode ci-dessus puisque chaque 
degré correspond une variable explicative. 

{en passant de r = 2 à r = 3 on rajoute une variable explicative}. 

D'un autre côté plus le degré du potynôme est élevé et plus on a de chance de se 
rapprocher de la vraie distribution, E y a donc un arbitrage à effectuer entre la 
lourdeur des calculs et Ja précision. 

Uac bonne politique consiste à minimiser le degré du polynôme pour un degré de 
précision acceptable. 
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« On peul utiliser les informations à pren disponibles sur la classe de fonction 
ostimée. 


— Ex. : Si on estime une fonction d'investissement : la répartition des 
coefficients retardés à tendance à monter, puis diminuer avant de remonter. 


Dans ce ças une fonction du troisième degré autorise deux extrermums locaux. 
Pour l'investissement : r = 3 est un hon point de départ. 


Pour une funétion d'exportation : ka répartition des coefficiènts à la forme d'un 
V renversé et un pulynome du second degré sera généralement suffisant. 


« On peut également suivre le pouvoir explicatif de l'équation (21) en fonction du 
degré choisi pour Le polynôme. 


Plus le degré du polynome est élevé, plus le pouvoir explicatif augmente, Mais 
cette augmentation peut être très rapide où au contraire (ente. 

Dans le premier cas, 11 semble que l'on ait un bon résulrat avec un degré faible. 
Dans le second cas, 11 peut être utile d'explorer des degrés plus élevés. 


R° 


Dans la pratique on dépasse rarement le troisième degré en essayant de maximiser 
l'écart entre r et n. 
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Sur le plan pratique pour obtenir une estimation par Lés retards polynormaux 
d'Almon, les bons logiciels se contentent d'une seule instruction (instruction 
PDL)}, qui demande que soient fixés le nombre de retacds et le degré du polynome. 
En comparant l'équation de départ ei l'équation estimée : 
Me ox +@x tt tax, (Si 
Y = Agir te FE, Lr (22) 
L'équaion {22} peut être considérée comme une version de l'équation (15} avec 
des restrictions. On peur donc essayer de valider globalement 14 procédure des 
relards d'Almon en utilisant La procédure générale de rest de restriction. 
L'équation (22) est nécessairement moins bonne en terme d'ajustement que 
l'équation (15), Dans quelle mesure est-clle moins bonne ? C'est à cette question 
que répond le test de restriction. 
Pour cela on calcule 
RSS(22}- RSS(15) 

_ nr 

Fa-r, 7-1 = RSS(15) 
T-1{ 
Soit Ex le seuil critique, on teste 
HS : procédure valide 
H : procédure non valide 
Së F calculé» Æ, -> Hi 
St F cakulé<#, >» HG 

RSS (22) dépend du degré du polynôme. 11 est donc possible de tester différentes 
variantes de dègré. 
(n - r} es! Le nombre de restrictions qui joue le rôle de q dans lé test précédent. 
“La procédure des retards polynomiaux d'Almon est d'emploi général : 
pratiquement toule fonctien qui présente des retards peut-être traitée par des 
cetards de type polynonuaux, C'est une procédure par défaut qui peut être utilisée 
quand il n’y a pas d'idées préconçues sur la véritable nature des mécanismes. 


Chapitre IL: L'estimation des modèles à 


éq uations simultanées 


Avec Les modèles à équations simulianées, de nuuvelles catégories de problèmes 
se posent, à côté de ceux qui or’ fait l'objet du chapitre L Ils conduisent à 
lutitisation de méthodes d'estimauon spécifiques dunt les mérites sont peut-être 
plus théoriques que réels. 

On commencera par présenter Les propriélés internes des modèles avant d'aborder 
la question de leur estimation, 


On peut définir un modèle comme un ensemble de variables organisées selon un 
système cohérent de relations mathématiques. 


C'est un ensemble formahisé, 
Cest un ensemble cohérent (intéressant Cn cas d'interactions ou de cu- 
détérnnnarions), 
Les modèles sont tion adaptés pour rendre compte des intérdépendances multiples 
entre les variables économiques. t'est une des raisons de Jeur utilisation. [s 
peruéttent de prendre en compte simultanément un ensemble de phénomènes qui 
doivent être étudiés conjointement. 
ILexisce : 

- dés modèles théoriques : 

-des modèles empiriques qui sont ceux qui nous intéressent principalement. 
Parmi les modèles empiriques on s'intéressera à ceux dont les paramètres sont 
stiniés économétriquernent. 
Le champs couvert par les modèles &sl très large. identique à celui de l'analyse 
économique elle-même. 1 existe donc des modèles macro-Éconamiques et des 
modèles micro-Économiques. 


L'avantage de cohérence peut être illustré par l'exemple suivant. Soit le modèle : 
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Cr = Loges) 
Hg.) (2) 
M=Gthta (y 


En travaillant séparément sur {1} et (2), men n’assure que les valeurs de c, et I, 
sont compatibles entre elles. 


Dans lé cadre d'un modèle (1} (2} (3) cette compatibilité est assurée 
immédiatement par la relation (3). 


- Historiquement, on a commencé par travailler sur des fonctions séparées ce 
qui imposait de vérifier la cohérence des résultats. 


— Ce n'est qu'ensuite que l'on a utilisé de grands modèles macro-économiques 
qui assuräient le résultat souhaité, tout en posant d'autres problèmes. 


Dans un modèle il y à donc plusieurs choses qui peuvent être classées en : 


. une partie formelle : des variables organisées en équations ; 
. une logique sous jacente : logique inteme que traduit lé système. 


IL faut done évaluer différemment el sclon dés critères différents : 
. la qualité des équations ; 
. Je caractère plus ou moins adapté à la situation de la logique interne. 


Une même situation économique peut être correcternent décrite par des modèles 
d'inspirations différentes. 


L'économétrie n'est pas assez sélective pour faire le tri entre les modèles, d'où 
l'analyse à deux niveaux. Le choix des modèles ne se fera donc qu'avec une 
analyse théorique. 

Ainsi par exemple la situation de l'économie française pendant les années 70-80 
pouvait être décrite par des modèles d'inspirations très différentes(1) : 


1} Sans compter un grand nombre d'autres modèles muins connus... dont la liste sorait à 
peu près impossible à établir : autres modèles construits par l'adtinistration, modèles des 
instituts de recherche publics ct privés, modèles universitaires, modèles des grandes 
institutions intemationales OCDE, FMI, modèles construits par les inglituts dé recherçhe 
étrangers etc. 
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Modèle de David (Banque de France) : Monétariste ; 
Modèle STAR (DP} : Néo Cambridgien ; 
Modèle DMS (INSEE - DP) : Néo Keynesien. 


Les trois modèlés impliquaient des politiques différentes, parfois même opposéts 
en terme de dépense publique ou de politique suluriale. 

I n'y a pas de relation assez étroite entre Un état de la réalité et un type de 
modélisation. 

1 peut arriver qu'un modèlc dont les équations sont irès bonnes en termes 
éconvmmétriques se révèle inférieur à un modèle dont individucllement les 
équations sont nettement moins bunnes, parce que la logique théonque du second 
est plus conforme à la situarion gui est décrite. Cetic évaluation n'est donc pas 
simpte à faire et ne peut se faire qu'en recourant à des critères variés concernant à 
la fois la qualité économétrique des équations, mais ausri La logique interne du 
modèle. Sur le premier point, certaines équations mal spécifiées peuvent tirec vers 
le bas la qualité d'ensemble d'un modèle. Cela pose le problème du traitement qui 
peut être fait pour certaines variables difficiles à modéliser, que l'on sait expliquer 
mais que l'on préfèrera considérer comme exogènes ou traiter hors modèles. 


Ce sont surtout les relations les moins bicn établies qui finalement détérminent la 
performance globale, 


SECTION 1 : LE MODÈLE COMME SYSTÈME D'ÉQUATION 


La partie visible d'un modèle, ce sont des variables organisées en &quations, On 
s'intéressera d'abord aux variables puis aux équations. 
A- Les variables et les équations 

Les variables qui figurent dans un modèle sont réparties en deux catégories : 

les variables endogènes ; 

. les variables exogènes. 
» Les endogènes sont celles dont la valeur est déterminée par le modèle lui-mêrne. 
C'est la résulntion du système d'éguation qui fixe les valeurs des variables 
cndogènes. 
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Certuines endogènes peuvent intervenir avec retard : ce sont les endogènes 
décalées. 


Quand on résoud le modèle, on te résoud en 1, par conséquent les valeurs des 
variables endogènes décalées sont déjà fixées et donc les endogènes décalées ne 
peuvent pus jouer Le rôle de vraies variables endugènes. Sur Le plan de la logique 
du système elles sont analogues à la catégorie des variables éxvgènes. Sur te plan 
économétrique elles posent des problèmes particuliers en compliquant La structure 
des térmes d'erreur. 


+ Les variables exvgènes sont des variables dont Ja valeur est déterminét on 
dchors du modèle. Pour certaines d'entre elles, elles sont exogènes par nalurc car 
ciles nv relèvent pas d'une explication économique (Exemple : Ja climatique, Ja 
démographie, ….). 

Mais les plus nombreuses sont des variables qu'on choisit de traiter comme des 
variables exugènes, c'est-à-dire qu'on ne cherche pas à les expliquer. De ce fait, la 
frontière entre variables endogènes et variables exogènes tradait le partage entre 
ce qui est expliqué et ce qui ne l'est pas, C'est un choix que Fait le constructeur du 
modèle. 

+ La prentière question à laquelle devra répondre l'écanomètre est celle du choix 
des variables et de Ja catégorie à Haquelles elles appartiennent, 


Parmi tes exogènes on distingue parfois : 


. les variables d'envirannement: 
. lés variables de commande, 


— Les variables d'environnement décrivent par exemple la situation de 
l'économie mondiale (prix de l'énergie, taux de change du $, etc.) 

— Les variables de commande sont celles sur lesquelles les autorités peuvent 
agir et qui décrivent des choix de politique économique. 
Cette distinction est utilisée quand on fait des exercices de prévision qu'en 
appellé des scénarios : on s'efforce dans un premier temps de caraceériser 
différents états possibles de l'environnèment : l'intérêt est de les caractériser de 
façon aussi nette que possible. Sur chacun de ces états on peut alors faire des 
simulations «vec des politiques économique qui sont elles aussi caractérisées. 
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Puis on évalue Le type de politique économique qui est adapté à un type 
d'environnement. 
C'est un exercice qui st pratique à l'occasion des principales décisions de politique 
économique : 


loi de finance: 
. préparation des plans à moyen terme, 


Parmi tes équations on peut distinguer deux catégories : 

J- les relations de Compurtement ; 

2- les conditions d'équilibre et les identités. 

Les relations de comportement 

Elles traduisent la réaction d’une variable aux variations d'une où plusieurs autres. 
Exemple : les équations (11 ct (2) 

Dans un modèle macra-£convmique, fa fonction de production, de consommation, 
la demande de monnaie, etc. sont des relations de comportement. 


La valeur des coefficients des relations de comportement résulic de la procédure 
d'estimation et c'est ce qui les différencie de la catégorie suivante. 
Les identités et conditions d'équilibre 
Les coefficients des identités et conditions d'équilibre sont, en effet, connus 
d'avance. 
+ Soit que cela résulte du la définition des variahles {cas des identités}. 
Exemple 

K,=K,1+1Inv, 
“Soir qu'il s'agisse d'une condilion gue l'on impose (cas des conditions 
d'équilibre}. Exemple : équation {3} du modèle précédent. 
La différence est que dans cenains &as ün peut introduire une variable d'écart sur 
ua condition d'équilibre et assacicr des résultats à cette vanable d'écart. Cela a un 
sens pour unc condition d'équilibre, mais n'a aucun sens pour une idenlité. 
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Un modèle complet va donc se présenter sous la forme de : 


n vatiables endogènes : 

m variables exogènes ou prédétérminées 
organisées dans h relations de Cotnportement : 
et k identités ou relations d'équilibre. 


Pour que le modèle soit dit fermé vu complet 4l faut qu'on puisse calculer de façon 
unique la valeur des variables endogènes. Les endogènes vont jouer ke rûle des 
inconnues dans le système d'équations. 

» Si toutes les équations sont de farme linéaire et indépendantes on sait que Ja 
condition pour qu'il y ait détermination d'une solution unique est que : 


Cette propriété doit être rapprochée du classement dés variables : touté 
augmentation du nombré d'endogènés entraîne automatiquement uné 
augmentation du nombre d'équations c'est-à-dire de la taille du modèle. 


Donc pour pouvoir Himiter la taille d'un modéle, il Faut limiter te numbre 
d'endogènes. Li y aura done un certain nombre de variables qui scront traitées 
comme exoènes. 


» Si certaines des relations ne sont pas linéaires on ne peut plus établir strictement 
Ja propnété précédénte. 11 faut raisonnér par élimination et dire que ss n£#h+kil 
est sûr qu'on ne peut pas avoir une solution unique alors que si Là propriété est 
vénfiée elie n'implique pas nécessairement une solution unique (des équilibres 
multiples sont possibles, cé qui traduit des propriétés particulières du système). 
Un modèle dont toutes los équations sont linéaires se résoud facilement par 
invorsion de matrice. 


Un modèle comportant des non-linéarités se résoudra en général par itérations : on 
se donne des valeurs d'amorçage pour certaines inconnues ct on regarde si ces 
valeurs peuvent être solutions, ce qui nécessite des moyens de calculs puissants et 
ne garantit pas contre l'existenec d'autres solutions locales associées à d'autres 
valeurs d'amorçage. 
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Quand la propriété n'est pas réalisée. c'est qu'en général F y a un problème de 
logique inierne dans le modèle : si par cx : n < bh + k alors en général il n'y aura 
pas de solution car on aura 1MpOosé trop de contraintes sur l'ensemble des 
endogènes ct sauf exception an aura vraisemblablement des mécanismes qui ne 
seront pas cohérents entre eux. Qn redonne alors des degrés de liberté au modèle 
en transformant des exagènes en endogènes : on endogénéise le nombre 
nécessaire de variables. 


Où peut illustrer ce problème par un exemple très simple : 


= 4 (4) 
= b:$ + br (5} 
M=c+é (6) 


On suppase deux endogènes : ÿ, ct ce, 
kin<h+k.lin'y a pas de solution. 
Une variable y, ici est déterminée de façon indépendante par deux mécanismes : 

. l'offre dans l'équation (4}; 

. Là demande par (55 et 6}. 
On peut alors redonner un degré dé liberlé au modèle en considérant l'une des 
deux variables LL ou 7 comme endagène, 
Si où choisit L4, y, sera dérerminé par la démmande dans {5) et (6) ; et (4) sert 
uniquement] à déterminer l'emploi. 
On peut aussi rencontrer la situation inverse : quand il y à trop d'endogènes par 
rapport au nombre d'équations. Dans çt cas il y a une indétermination : les 
solutions ne sont plus déterminées de lagon unique. il existe un infinité de 
solutions 11 faut donc réduire 1& nombre d'endogènes en canstdérant certines 
variables comme exogènes, 
Un modèle s'écrit équation par équation et donc il n'est pas évident quen=h +k. 
On nc peut que le vérifior à la fin. 


De plus, la variable de gauche n'esl pas forcément une endogène {ellé lé sera 
souvent cependant). 
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Pour faire apparaître les propriétés internes d'un mudèle on s'intéresse à deux 
choses : 

- L'organisation de la causalité dans le modèle. 

- Les propriétés dynamiques qui sont incorporées dans le système d'équation. 


C'est la traduction de la logique contenue <lans le modèle. 


B- La causalité 
La causalité porte sur les relations instantanées entre variables, 
1l existe différentes approches du phénomenc de causalité, Cela est vrai aussi en 
langage quantitatif, 
En éconamétrie on retient surtout la définition de Granger qui cepose sur l'idéc 
que si x et # sont deux processus stationnaires. x est considéré comme cause de y 
si Ja prédiction de y qui est faite à partir des valeurs passées de v elle-même et des 
valcurs passécs de x est meilleure que la prédiction de y faite uniquement à partir 
des valeurs passées de y. 
Cette notion a donné narssance aux tests de causalité usuels. 
Cela s'applique bien à des modèles à une seule équation. 
Exemple : relation entre les prix et la masse monétaire. 
Dans le cas des modèles à équation simulranée on retient une autre définition qui 
cest restreinte à l'ordre de délermunation des variables endogènes. Ce qui implique 
que l'analyse de la causalité sera restreinte aux soules endogènes. On retiendra 
comme définition que y; est cause de y; si y; est détérminée avant et sert à 
déterminer y. 
Pour faïre apparaître les relations de câusalité 1} faut donc faire apparaître l'ordre 
dans lequel seront déterminées successivement les différentes variables 
endogènes, ce qui revient à faire apparaître l'ordre dé résolution du motèle. On 
raisonne toujours dans le ças d'équations supposées linéaires et on utilise les 
propriétés des systèmes et sous-systèmes d'équätions linéaires. 
On peut donner une représentation simple de la méthode. Pour cela utilisons un 
exemple. 
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Soit le modèle suivant : 


(7) C= as, + à 

(8) 1 = buim,+biyt + b3 
{9} Im =cix-âr 

{10) tlmec+l+x 
{11) = dé 


ct = Consommation 
Yu = ICVETUS 
Im, = inporlation 
X4 = CXpUrATIONS 
l, = investissement 
= emploi 
Ar, = variation de réserves cxtéricurcs 
Il y à donc cinq équations. IL faut cinq variables ondagènes : 
CELA LRU 
Le modèle est donc fermé c'est-à-dire qu'il admet un ensemble de valeurs uniques 
pour és vanahles cndogènes. 
Les éxogènes sont : |, ye1, x, An 
Comment faite apparaître un vrdre dans les causalités ? 
Définition 
Un sous-système d'équations linéaires est dit rimimum sil contient un nombre de 
variables êt d'équations suffisant pour le déterminer. 
Un grand système d'éguation peut ou non contenir des sous-système plus petits. 


Déterminer là causalité revient à hiérarchiser [es sous-systèmes par taille 
croissante. 


Pour cela on peut s'aider d'un tableau dans lequel figurent en ligne les équations et 
en colonnes les variables. 


92 Économétrie appliquée 


Ô= non {8) 0 (+) l | 0 
| û 0 û 0 

‘ {0j 1 1 1 | 9 

{11} | ÿ il LL 0 0 | l 


On peut donc déterminer les sous-systèmes : 

(9) ne contient qu'Img — c'est donc le plus petit sous-syslème car il n'a besoin 
d'aucune autre information pour fournir son résultat. 

C'est un sous-système d'ordre 0. 

On peut réécrire le tableau en permutant les Tignes er los colonnes de façon à 
placer (9) et (8} aux deux premières lignes et Im, et 1, aux deux premières 
colonnes. 


Im,ihiy C, Î, 
@) 110:0 © 6 
{8 11:11:10 0 0 
A) Tilri 1 0 
7) 8:01 1 0 
{D 701011 6 1 


Connaissant la solution de (9) an en déduit la solution de (8). 

On a alors un autre sous-système qui est d'ordre | 

(10) et {7) constituent alors un sou$-système d'ordre 2. 

(11) constitue un système d'ordre 3 

L'ordre dans lequel vont être déterminées les 5 variables est donc : 
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lei, on à en fait une structure de pré-ordre. (La matrice n'est pas triangulaire). La 
méthuce s'apparente à une diagonalisation de la matrice des (0,1}. Quand on x une 
structure de pré urdre c'est-à-dire que des variables sont déterminées ensemble on 
ne peut pas diagonaliser la matrice. Si on peut la diagonaliser on obtient un 
système récursif qui à comme propriété d'être entièrement hiérarchisé sur la 
toralité des variables. 


Un ensemble d'équaliun peut en fait contenir des logiques différentes. Conserrons 
ke même système mais considérons les variables différemment. Considérons que 
ls cing cndugènes sont ct, ÿt. Ir, 1mt, &re on prénd 1 comme exogènce : on 
impose un niveau d'emploi}. Le modèle va donc incorporer une logique assez 
différente. 
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a —_……. 


Faisons apparaître un sous-système d'ordre D qui est {Lt}, on a: 


sous-système d'ordre 0 


sous-système d'ordre 1 


sous-système d'ordre 2 


+ Une même structure formelle peut donc incorporer des logiques différentes, On 
pourrait encoré trouver d'autres logiques en considérant xt comme cndogène par 
exemple. 


À partir de 1à on peut faire apparaître une typologie assez sommaire des systèmes 

d'équations : 

— [y a d'abord les systèmes décomposables qui sont des systèmes dans 
lesquels une partie des variables est déterminée dans une partic des équations 
et le reste des variables dans les autres équations. Un modèle décomposable sc 
présente en fait comme la juxtaposition de deux sous-modèles indépendants, 


L'estimation des modèles à équations simultanées 95 


Le modèle peut donc se décomposer en deux sous-parties indépendantes, 
Historiquement cela Cotrespond à un modèle intéressant, lc modèle original de 
Harrod Domar. 


Le dépassement de ce type de modèle se fait en établissant des relations entre Les 
deux parties. Dans 1e modèle Harrod Domar en endogrénéisant le taux d'épargne 
par EXCmpIC. 

— Les modèles ou systèmes de causalité : ce sont des modèles dans lesquels on 
peut au moins partiellement faire apparaître l'enchaînement des causes et des 
effers. Notre exemple correspond à un système de causalité. 

— Les Systèmes d'interdépentance : ils constituent dans leur totalité Le plus petit 
sous-$système admettant une solution unique. Cela revient à dire qu'on ne peut 
pas hiérarchiser l'ordre de résolution des variahles : l'ensemble des endogènes 
est déterminé simultanément au moment de la résolution du modèle, Sur le 
plan de 1a logique, l'interdépendance peut s'interpréter comme l'existence de 
relation de retour (leed back) dans 1e schéma de causalité du type : 


C 
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a ——_—— 2 —— — - .—— 


Yi 


Y] houcle ou relation de feed back sur l'arbre de causalité 


À 


[] faut malheureusement remarquer que cette technique d'analyse est peu 
perfonmante cat la plupart dés modèles ou système de taille importante sont des 
systèmes d'intérdépendance. Plus La taille du modèle augmente plus on à de 
chance d'avoir un système d'interdépendance. 
Une autre faiblesse est que la méthode se bloque rapidement quand l& nombre 
d'équations et de variables augmente. Ce qui est facile à faire apparaître avec un 
nombre restreint d'équations devient beaucoup plus complexe avec un grand 
notnbre d'équations et de variables. 
+ T1 pout alors être utile d'élaborer des maquettes du modèle de taille réduite pour 
analyser de façon plus commode la causalité ou Les propriétés dynamiques. 
Les propriétés dynamiques des systèmes d'équation 
Par définition on peut dire qu'un modèle est dynamique si dans un environnement 
donné qui est caractérisé par la stabilité des variables exogènes, ur état de 
l'éconunie en t entraîne un autre état en t + 1 fun état de l'économie étant 
caractérisé par un ensemble de valeurs des variables cndogènes). 
Sicnti+ Iles valeurs sont les mêmes qu'en t on à un Étal stationnaire. 

Soit Y = l'ensemble des variables endogènes. 

Yi pour X constant (X = ensemble des exvgènes). 
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Sur de plan de l'écriture (plan formel) un modèle est dynamique s'il contient dés 
variables ondogènes décalées. 
Tout modèle qui inclut des vanablus endoyènes décalées possède des propriétés 
dynamiques qui vont sc traduire par des trajectoires suivies par Les endogènes au 
cours du temps. 
C'est pourquoi les propriétés dynantiques peuvent être analysées et conslituent un 
critère d'évaluation du système d'équation. 
Parce qu'on sait que dans la plupart des cas les phénomènes économiques 
évoluent assez lentement, on demandera le plus souvent au système des propriétés 
minimum de stabilité. 
{Par cantre si on veut modéliser un système très instahle on privilégiera une 
structure qui a cette prapriété}. 
Qu peut donner différentes définitions de la stabilité qui sont plus ou moins 
pertinentes selon le type de modèle utilisé. 
Sachant qu'on à : 

Y l'ensemble des variables endoygènes 

X l'ensemble des variables cxogènes 
On appelle Yet X! les valeurs de références des variables endogènes et des 
variables éxogènés (valeurs du compte cençral). 
La stabilité statique 
Si on applique un écart pelit aux valeurs des exogènes X-X! < © alors on modifte 
neu les valeurs des endogènes : Y-YT <fÿ 
Le système est amorti, il a une forte incrtic. 
La stahilité dynamique ou asymptotique 
On dira qu'un modèle présente cette propriété si après un écart X-X! porit : X- 
XT<u, le vecteur des endogènes revient vers le compte central quand Fe temps 
passe : 


Y— Y'quandT ?+6 
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Le multiplicateur Keÿnesien peut servir à illustrer cetle propriété : si on applique 
pendant un notnbre limité de période une variation du niveau des dépenses 
autonomes, lesqueiles sont par définition exugènes, alors 1e revenu, etc. 
s'éloignent du compte central puis + reviennent progtessivement. 

La stabilité en trajectoire 

Un système est stable en trajectoire si après un écart X-X'< @ les valeurs des 
variables endngènes rejoignent une trajectoire parallèle à celle du compte central. 


Le même multiplicateur Keynesien quand on applique un changement permanent 
de la dépense autonome : revenu et consommation vont suivre la même évolution 
mais en dessous ou au dessus du cœmpic central. 

Exemple d'analyse des propriétés dynamiques d'un modèle : J'oscillateur de 
Samuclson. 


Comme on sait, if peut être représenté par : 


= ar + &Yel {12) 
L=b1 +b3 (y - Yr2) (13) 
= + +G {14) 


En remplaçant (F2) et (13) dans (14) où obtient : 


Fes ai + 2ÿr) + bi + b2ÿi1- eo) + Gt 
Ye= (ao + ho} y - Do vis + a + bi + G 
Yu (ao + ba) yes + ba vs = a + bi +G, {15} 


C'est donc une équation de récurrence du secotul ordre avec second membre. 

La solution est donnée par y{ = À! et pour calculer l'équation caractéristique on 
remplace à gauche dans l'équation sans second membre et on divise par AT: 

À = (as +) BAT 2 20 (16) 

À-(as+b;)Ath=0 (17) 

La solution générale est donnée par : 


3 = AA + 4,46 où A, et À, sont solutions de l'équation caractéristique {17}. 
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À et A3 sont des valeurs calculées en utilisant la partic droitc de l'équation 
(Second membre). 


Les propriétés dynamiques de la solution générale dépendeni des valeurs de À. 
l- Si jet À: sont réels, le système engendre des évolutions dynamiques 
IMONHUNES 

la- [A1] et (À 2|<1— évolutions monotones amorties (— stable } 

LE [2.1 ou [à 2 >] évolutions monotones divergentes ou explosives 

2- Si Ajet À complexes, alors le système engendre des évolutions cycliques. 
Za- modules de À, et À, <1— oscillations amorties 

2b- moduiss de À, et 2, >1l— oscillations explosives 

Sait sur l'exemple de l'oscillateur : 

À = discriminant = (a + b2)2- 4 ba 

Étudions Le signe de A : 

A>0 ai :(@7 +b3)° > 4h; 

Soit a +b3 > 24h 

a >24{b 2 (8) 


Calculons quelques points de La fonction (18) : le maximum ct atteint par : 


da 1,172 
M9, 12 109 
db: 2 
Soit —=——1=0 d'où b)=1 
vb 
{b;=l &5 =i 
si ba =0 a; =0 
ba = à a3 =0 


+ Graphiquement 


Au-dessus de la courbe îes solutions sont réelles donc lés évolutions sont 
monotones et en dessous ellés sont complexes, donc tes évolutions sont cycliques. 
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Il 4 


* Pour voir si les évolutions ont un caractère amorti où explosif, on commencera 
par raisonner dans fe cadre de racines complexes. On sait qu’un nombre complexe 
peut se mettre sous la forme : 


picosô+i sin6) 

p= module  8= argument 

Quand on élève un nombre complexe à une puissance t on obtient : 

p'(cos r8+i sin 48 ) 

Deux complexes conjugués ont ta même partie réelle mais deux parties 
imaginaires opposées suit : picos 8 +isin Pet pécos(-8}+rsint-#)} 

Le produit de deux complexes conjugués donne : 

p(cos? 8— 57 sin? 8}= p°? 

Dans je cas où les racines À, &t À: sont complexes, ce sont des complexes 


conjugués et le produit de ces deux camplexes vaut alors p° c'est-à-dire le carré 
de Jeur module commun. 


Mais on sait aussi que À,1, =b: dans l'équation (17} 
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Donc si b3 > 1 alors p?> let dans cc cas Ics oscillations sont divergentes on 
cxplosives, 
Sib3 < 1 alurs p< 1 = cas n° 2a usoillations amortics. 


On peut représenter graphiquement les différents cas possibles : 


4 1b 


b 


Si les solutions sont des nombres réels, on peut utiliser la propriété qu'un nombre 
récl est un nombre complexe dont la partie imaginaire esl nuffc, ce qui permet dc 
retrouver la mêtme propriété et la même conclusion. 


L'analyse des propriétés dynamiques peut se faire par une méthode plus générale : 
+ Lu méthode générule 

Elie utilise les propriétés des équations de récurrence matricielles. 

Soit LA un vecteur de variables endogènes et Y-1 ce même vecteur décalé d'unc 
période. 


L'équation Y= AFP, (19) constitue nne équation de récurrence matricielle du 
premier ordre, À étant une matrice régulière de format (n x n}. 
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La solution d'une telle équation est donnée par: Y =r y: (20), V étant un 
vecteur de même format que Ÿ el r étant un scalaire, 
L'équation caractéastique de cette éqtiation de récurrent s'écrit : 
tV=Arlv (21) 
soitrV=AV (22) 
ou(A-rV=Q (22his 
Cette équation caractéristique exprime que r est une valeur propre de la matrice A 
et V est le vecteur propre qui lui est assucié. 
Si là matrice À est une matrice régulière alass 1] # aura autant de valeurs propres 
que de lignes ou de colonnes dans À soit n valeurs propres. 
[A-ri|=0= n valeurs propres 11... 1. 
La solution générale est alors : 
Y = A qu") #+...+ An r LE (23 
—  Lesri… 1, étant les valeurs propres de À 
— Les A, … A, étant des coefficients de pondération qui indiquent l'importance 
respective des différents modes dynamiques indiqués par r, à rs 
Ces coefficients sont calculés par rapport aux conditions de dépait : en faisant 
t=0onobtient Yo=A;Yi+...+ AV, 
y ici n équations et n inconnues à calculer, ce qui permet de calculer 
les valeurs des cocificients Al … À, 
Chaque valeur propre carattérise un mode d'évolution dynamique. 
Qr, il y a autant de valeurs propres qu'il y a d'équations c'est-à-dire de lignes, et 
denc autant de valeurs propres qu'il ; a de colonnes c'est-à-dire de variables. 
La dernière étape consiste alors à associer des propriétés dynamiques c'est-à-dire 
des valeurs propres à des équations c'est-à-dire des propriétés structurelles du 
modèle. 
+ Cela permettra de savoir quel mécanisme est à l'origine de telle ou telle propriété 
dynamique et ça permet aussi de modifier les propriétés dynamiques du modèle 
qui ne conviennent pas en agissant sur la Cause. 
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Appliquons cette méthode à l'exemple précédent : 


{12} Ce à] + à Foi 
(6 E)) L= by +b2 Gr - 12) 
(14) = € ++ Gt 

C? 


Le vecteur Y = | 1, 
Ÿr 
La difficulté cst la présence de Y..2 qui impose de faire un changement de variable. 
On considère que (14) est satisfute en toute période d'où on peut éctire ! 
Yet eee tt + Ge (24) 
On peut définir une nouvelle variable z, telle que : 
= #hei +81 (25) 
et Z1=G-2+h.2+G 253 2 (26) 
Le 
Posons Y =| {, 
L 
Ce vecteur constitue le nouveau vecteur des vanahles entogènes. 
En remplaçant y,2 par sa valeur tirée de (26), l'équation (13) devient : 
L= br + bo ce + bol + boG1-h2z) (27) 
Donc, en remplaçant v+.1 et y, le système devient : 
= ac +sl 1 +a5G,_, (28) 
L=h+bsc +bl,1 +b2G, 1 Bou à (20 
a =Githt+G_, (25) 


que l'on peut mettre sous la forme Ÿ = AY 
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Cr DE l 

Y=|£, Y-1=| 1: 
Z êe. 1] 
dj 3 Q 

Ab; b2 -b 

(3,3) ; 1 0 

Donc 


F A F] {32) 120 
dau de 


autre LÉ membre 


e Calvul des valeurs propres 
det{tA—r/}=0 


UT  & 0 
lA-ri|< b, Pr —b; {29} 
1 1 —+ 


detéA.-r/}=-r(as -r)(b; - #) wub thoag rt ashsr 135) 
=-rfas-r)(D; -r)< bor+anbor 
= -rr? (a; +b)r+b 50) 
r=N0 
detiA-#i)=0 4, 
r—(as +bsdr+b; =0 (31) 
L'équation (31) est identique à Féquation (17) étudiée précédemment. 


Comme le montre cet exemple, dans la plupart des cas la mise sous forme 
dynamique Ÿ = AŸ., ne se fait pas directement, Célia impose des transformations 
sur le modèle de sortt que l'analyse des propriétés dynamiques se fait surtout sur 
des versions dérivées c'est-à-dire retravalllées. 
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Deux exemples permettront d'illustrer l'apport de l'analyse des propriétés 
dynarmigues d'un modèle, 
Exemple 1 : modèle de Klein pour l'éconemie US : {1} 
e =16,78+0,0211+0,231L, + 0,80(wi +w2) (32) 
1=17,79+0,2301+0,5S01: -G1SX ; (33) 
wi =1,60+0,42(Y +T-Wwa)+0,16(Y+T-w3) 1 +0,13(6—-1931) (34) 
F+T=C+1+6G quasi condition d'équilibre, (35) 
F=lltmtw (QD 
AK =1 (38) 
IT = profit des entreprises 
Wi = salaires versés 
K = stock de capital 
T = impôts 
t= indice du temps 
[= investissement 
c + consommation 
Après (ransformation ce modèle peut être mis sous la forme : 
[T 0,863 -0,063 (ltd 0 î à I 
Y 1,489 0,174 -0,283 0 Or; 
K 0,746 0,015 6,816 O0 
c|"|a,743 0,189 -008 0 
H] 0,626 0,237 0,119 © 
1 | 0,746 0,015 184 0 


oÏK, 
As +Z (39) 


(e] : 4, 
ol 1 
On constate que sous cetle forme, le modèle présente la particularité d'avoir trois 


colonnes nulles, ce qui montre que les trais dernières variables n'interviennent 
dans aucun processus dynamique. 


9CES CC 


(1) Cet exemple est tiré de Bridge "Applied Econometrics" 
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— Les éléments dynamiques sont donc situés dans La partie haute de la matrice et 
Fanalyse se réduit à cette parlie. 


0,8R63—r  —-0,063 0,164 
Calcul du déterminant | 1,489 0,174-—-7 0.283 1=0 (4) 
0.746 0,015  (0816-7r 


ñ = 0,321 {mude monotone amorti) 
ines :45 =0,77- i : ; 

ERA à : ee "x Ses imodes ostillants amurtis} 
On met en évidence les propriétés du modèle mais on pourrait aussi se demander 
les raisons des différents modes. C'est ce que permet d'aborder l'exemple suivant. 
Exemple 2 : modèle DECA{D? 

K=C +1,+$,—1lm, (condition d'équilibre) {41) 

Im, =, +ail+as hs —bil.1 (42) 

€, =R, +ax +2 + CR {43} 

l,=udY, (44) 

Si=e[r-fr] (45) 

R= gl +82 tesks (46) 


St= stock CL = consommation 
Im, = importation I, = Investissement 
R, = revenu Y,= PIB 


Certaines de ces équations conticnnent des propriétés dynamiques. C'est Le cas des 
équations (42) (44) (45) et (46). 


Après chiffrement des paramètres on obtient sur le modèle de départ trois valeurs 
propres distinctes : 


(À ff s'agit d'une version simplifiée de DECA. L'exemple est tiré de l'ouvrage de Deleau et 


Malgrange "L'Analyse des modèles macrn-éconnmiques quantitatifs", Economica 1978 page 
100 et suivantes. 
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ñ = 0,26 

r =0,31 

r =1.22 
On a donc 


. un mode altermatif (+ - + -): 

. un mode monotone amorti; 

. un mode divergent. 
On peut se poser la question de savoir quelles sont Les équations qui engendrent 
ces modes £t comment on peut les modifier. 
Si on modifie l'équation (46). Cette équation contient des vanables décalées qui 
lui donnent des propriétés dynamiques. On abtient alors : 


ñn = -0,26 
r = 0.36 
r = 0 


On peut aussi modifier l'équation (43), On obtient : 


ñ =-0,27 
rm =-0,02 
r =1,23 
On peul aussi modifier (45). On obtient : 
n = 0.06 
r = 0,35 
ñ =1,2} 
Enfin on peur modifier (42) et on obtient : 
n =-0,53 
ñ = 0,29 


7 — 1,23 
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À chaque transformation on peut associer uné valeut propre qui à été fortement 
affectée. Le mode d'évolution monotone divergent caractérisé par la valeur propre 
1 est associé à l'équation (46). 


Le mode d'évolution monotone convergent caractérisé par r> est associé à (43) 
c'est-à-dire à la fonction de consommation. 


Le mode alternatif caractérisé par r, ést associé à l'équation {45) c'est-à-dire au 
mécanisme de formation des stocks. 

L'équation (42) ne change pas là nature des mudes d'évolution et donc cette 
équation n'est associée à aucune dynamique en particulier. Ce n'est pas surprenant 
car on avait quatre équations pour 3 racmes. [l y avait donc une équation de trop 
qui semble être ici (42). 

Qu'il s'agisse de l'analyse des relations de causalité ou de l'analyse des propriétés 
dynamiques on se heurte vite à des difficultés pratiques importantes dès Jors que 
la taille des modèles augmente ce qui fait qu’on ne peut pas analyser ies modèles 
cn grandeur réelle mais seulement des versions réduites et surtout transformées 
pour Jes besoins de la mise en forme {type Ÿ = A Y ). 

Les deux exemplès qui ont été présentés sont des modèles anciens, mais les 
versions transformées sOn {ÈS mauvaises. 

Des études ont été faites également sur lé modèle STAR pour la France. 

Quand il y à des non-linéarités (rapport de variable par exemple...} les méthodes 
qu'on a vues ne sont plus directement applicables. Au total, l'analyse des propriété 
dynamiques, comme du reste celle dé la causalité, reste un domaine difficile et 
finalement assez peu satisfaisant, en dehors de la typologie un peu simpliste des 
cas de figure possibles. 


SECTION 2 : L'ESTIMATION DES MODÈLES À ÉQUATIONS 
SIMULTANÉES 
Comme on l'a vu un modèle se présenté cernme un ensemble de variables : 


. ni variables endogènes : 1... Yn: 
. M variables éxopènes : x) … Xn: 
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organisées dans un ensemble de k + k = n équations. 
Le modèle peut s'écrire en isolant Les variables endogènes à gauche : 
Pur = 2) 2 te ta Var pts + bip Ÿ M 


ne = pl} te nn -(Yn-)s + bike +. + Drm Xait +ünr 


C'est la forme qui vient généralement en premuer lors de l'écriture du modèle et 
qui est censée traduire comment les variables réagissent les unes par rapport aux 
autres. On l'appelle forme structurelle, On peut aussi isoler à draite les termes 
d'erreur en respectant le même ordre pour [es variahles endogènes dans chaque 
équation. 


La forme structurelle s'écrit alors : 


Me M2 2 aa Xe bidon SU 
Ant ar À dore Ga Vas — hote am lon © M2 


Tri is À Vae Paie om X ot = 
On peut réécrire le modèle en notations matricielles. 
Soit Ÿ= l'ensemble des variables endogènes : 

1 
Kn.) = 

Ys 
Soit X = l'ensemble des variables exogènes : 

A 

Xm,1) = 


Xm 
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4] 
Soit Len, 1) = 

CA 
Soit À : matrice des coefficients des variables endogènes dans le système 
d'équations : 


LG. —-&x 

—47 | “dy 
Ajan) _ i j 

nl 1 
De même 

UT ‘ bin 
Enr) un 

ba — bn 


La forme structurelle peut alors se réécrire 


AY +BX=UC . 


Remarque : Dans le chapitre 1, Yet X désignaient respectivement Les valeurs 
prises par li variable expliquée et les variables cplicanes au cours des T 
périodes utilisées pour l'estimation. 


{et Y et X représentent l'ensembles des variables endogènes et l'ensemble des 
variables exogènes du modèle. 


A- Le biais des moindres carrés ordinaires 
Un premier résultat important est que l'estimation d'un modèle à équations 
simultanées par les moindres carrés ordinaires conduit à des estimateurs hiaisés et 
non convergents. 
On peut le montrer à l'aide d'un exemple très simple en s'inspirant d'unc 
démonstration proposée par Haavelmo. 
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Soit le modèle suivant : 
ray, +h+u, (]) 
É =C+z (2) 
En ee (2) dans (E) où abuent : 
Gr e=— Ts* D me 
Si on reporte dans (2} il vient : 


1 
& = 
Ra en Lee 


D'après (4), la variable explicative y, n'est pas indépendante du terme u dans (1) 
et donc l'hypothèse H6 n'est plus vérifiée. 
b lis 
E(y,u,)= dc nr | ea {5} 
L'estimateur des MCD dans (1) 
ë= 2 EX ET) {6} est biaisé et non convergent. 
2 - 5) 
On peut pour cela remplacer c et y dans : 
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és fs D (7) 
1-4 1 


Le dénominateur vaut 


L = M F 
El —— (2-7) 2 


l-a 1-4 


1 =, î 
Er - x£ + 
de eq 


ac? +02 +(a+licovzu 
Soit Ê EE à à 
O; +03 +2cov zu, 


è 2 ï = 
7 Yu + 2 2}, 


L'estimateut à est donc biaisé. De plus il n'est pas convergent c'est-à-dire que le 
biais demeure même si le nombre d'observalion augmente indéfiniment. 


La limite en probabilité de â est : 
. . a0?+6? | 
Plim = —$—ÿ" donc ne converge pas vers a même pour T grand. 
gi + 


La cause du problème est bien la non-vérification de H6. 


En effet, si y;ne contenait pas u, il n'y auçait pas au numératcur et au 


dénominateur de termes en oc? 


+ Quel est Le sens du biais ? 
ac? +o2-a02 02 {1-a)o? 
TE rt. 


Piimäâ-a= 
o? + 6? 6? + a? 
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Si donc : a>l L-a<0 Plim à < à — on sous-Cstinc 
a<l l-a>0 Plim à > à — on surestime 
Dans l'exemple de Haavelma : 
On surestime la propension à consommer. 
Il n'est donc pas correct de traiter un système à équations simultanées par les 
MCO directs équation par équation. 
Il Éuut donc recourir à des méthodes d'estimation spécifiques. [l en existe plusieurs 
dont 2 présentent un intérêt particulier. 
B- Les moindres carrés indirects 
Soit : AY + BX=UC 
Siilexiste Al alors il est possible d'isoler Y : 
Y=-ABX + AU 


On pose 
C=-A'!B 
et 
V=alu 


Alors il vient 


Y=CX+V 


Dans cette nouvelle forme les variabies endagènes c'est-à-dire Les y n€ dépendent 
plus que des variablés exogènes, les x et donc il n'y à plus de biais de 
simultanéité. 

La malrice € est de format (n x m) et V est un nouveau vecteur dé variables 
aléatoires de format (n x 1}. 


Si on développe cette forme réduite, elle s'écrit : 
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ir = EnAe-Cimdnur ŸVir 
É 

, 

l 


Var = CE ee Con tnt À Va 


dy; 

dans le système : c, = RAA 
l dx; 

F 


Les c,; sont alors les multiplicateurs d'impact des variables caogènes sur les 
variables endogènes. 
La méthode des moindres carrés indirects (MCT) consiste à : 
- Dans une I étape : estimer les coefficients de la matrice € par les MCO 
— Dans une 2F étape : remonter de € vers Ta forme structurelle c'est-à-dire vers 
les coefficients de A et de B. 
C'est l'apération d'identification du modèle sous su forme structurelle. 


Comme le montre la relation C = -A°TB, lcs cocfficients de € sont des 
combinaisons des coefficients de À et de B. 


| : a, by bi ] 
. [l 


a Ain ni Braun) = 


b, [ | LP ! 


nl il 
El faut donc dans À calculer n {n - 1} coclficients inconnus. 
Dans B il y a {n x m) Cueilicients inconnus. 
Cela Éait donc (n - 1} 0 + nm cocfficients à calculer à partir des éléments de € qui 
sont au nombre de (n x 111). 
Dans le vas général ce n'est pas possible cat il y à trop d'inconnues par rapport au 
nombre de relulions. Le modèle ne peut pris être identifié. 
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L'identilication du modèle n'est possible que si il existe un certain nombre de 
restrictions à priori : n fn - 1) restrictions à priori qui permettent de réduire ke 
noënbre de coefficients à estimer au nombre de relations c'est-à-dire à am. 
Les restrictions à priori sont de nature différentes de celles qui ont été présentées 
au chapitre T. ici ce sont toutes les informations qui permettent de réduire le 
nombre de cocfficients à calculer. Ces restrictions à prion pcuvent être de 
plusieurs types. 
1- Le coefficient est connu d'avance (cas des conditions d'équilibre et des 
identités} donc il n’y à pas besuin de le calculer. 
2- Une variable ne figure pas dans l'équation. Le cocfficient correspondant cst 
donc nul. 
3- Les relations possibles entre coefficients : exemple : un coefficient est la 
somme de deux autres. 
De plus, l'identification ne se fait pas sur l'ensemble du modèle mais équation par 
équation. Le nombre de restrictions doit donc être équitablæment réparti sur les 
Squations soit {n - 1} restrictions par équation, 
Dans Le vas où les restrictions à prion prennent la forme de variables ahsentes, 
l'identification n'est possible qu'à condition d'avoir (m+1) variables dans chaque 
équation. En effet, soit l'équation k, Cette équation à au maximum n + m 
variables, On lui impose n - L restriction à priori sous la forme de variables 
absentes. 
Il y a done in + m)-ût- 13 m + 1 variables présentes dans l'équation. 
On comprend aisément que cctle contrainte ést inacceptable car beaucoup trop 
stricle, en macro-écunomis par EXxcmple. 
La plupart des modèles seront donc soit sous-idencifiés, soit sur-identifiés. 
-—  Sur-identiftés quand le nombre de restnclions est > {n - 1) par équation. 
—  Sous-identifiés quand le nombre de restrictions est < {n - 1} par équation. 
L'exemple très simple des équations {1} et (2) permet d'illustrer ce problème. 


En (1) 
Y=G+tz (2) 
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Sani AS 


Ce modèle comprend : 


[ef 
+ deux endogènes donc Y = | : | É 
Yi 


* deux exogènes donc X = | . : 
L 


|" 
“etU= . 
(4) 
En isolant les termes d'erreur à droite, il vient : 


€, —aÿ-bru, 
—G +} =0 


D'où a=[} ri et 8-[° Ÿ 
Le modèle peut donc se réécrire : 
NGC TZ LM sav esxeu 
US cp l4 Flo oit AY + BX = 
Ceci est une autre façon d'écrire la forme structurelle. 


Pour la première étape des MCL il faut calculer La forme réduite : 


Pa #4 


La Hu 
à He 0 Aa ie 
: | 


PE AR 
Pen au 


AT! 


1 b 

1 -a 1-2 

C'est une matrice dont les éléments sont des combinaisons des cocfficients de A et 
B comme indiqué précédemment. 
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: Aa #-u 4, | ven 
: H=é K-a F f., 


La forme réduite est : 


Ike Wa EL D) 
| PAL 


Qn vérific bien que la partie droite ne contient que des variables exogènes. 


Peut-on idenûütier ta furme structurelle, c'est-à-dire connaissant ta valeur des 
coetficients de €. peut-on remonter vers les coefficients de À ct de B ? 


Il y à quatre relations constituées par les 4 coefficients de C, pour calculer : 
2 cocfficionts dans À 
4 coctficicnts dans B 
Le nombre de restrictions à priori sur À et B: 
dans À : L restriction (le coclficicnt de &, est égal à +1 dans 
l'équation (2} ; 


dans R : 3 restrictions (] variable absente dans chaque équation et Le 
cosfficient de z, dans (2}}. 


Nombre de coefficients restant à calculer : 6 - 4=2 
Nombre de relations + 4 
Le modèle est donc sur-identifié globalement, 


+ Cela signifie que plusieurs coetficients peuvent être calculés par plusicurs 
mécanismés indépendänts 


n es 
Appetlons c,; les éléments de C: |" 91 
C4 CA 


Chacun des coefficients a et b peut être calculé de plusieurs façons diftérentes. 
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La 
€ 
Exemple : a =“ 
ES 
c21 
nl x 
Les SE = 65 
nv € 
Pat 
de même b=2 
#“ 
C21 
cha 
et b= << 
€21 


D y a donc diftérentes façons indépendantes de calculer à et b. C'est ce que veut 
dire la sur-identification : trop d'équations par rapport au ncnbre d'inconnucs. 
Cette question de l'identification permet d'expliquer un problème important qui 
concerne l'économétrie des marchés où économétrie de l'offre ct de la demande. 
Supposons qu'on chetche à estimer une fonction de demande et que pour cela on 
dispose d'un certain nombre d'obsérvations sur le prix êt la quantité échangée au 
cours de plusieurs périvdes. 


P 
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] 
ia d'où ae 


{21 
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Un premier cas de figure serait celui dans lequel les fonctions d'offre et de 
demande sont stables. À chaque période la quantité échangéc sera proche de 
l'équilibre. 

4, = quantité vflerte=a+h? (9) 

Gp =c+dP (10) 


On obtient un ensemble de points groupés autour de l'intersection des fonctions 
d'offre et de demande et il est Impossible d'esrimer ces fonctions. 


Un deuxième cas est celui où l'unc des fonctions se déplace : 


P 


O1 


O2 


D1 


Q 


La fonction d'offre se déplace ct l'intersection décrit la fonction fixe, c'est-à-dire la 
fonction de demande. 

Dans ce cas la régression de la quantilé échangée sur le prix d'équilibre donne 
bien unc fonction de demande, 
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Lc cas symétrique est celui où l'affre est fixe ct où Ta demande se déplace : 


P 


O1 


D1 
DÜ 


Q 


Le point d'intersection décrit la courbe d'offre et donc ia régression de Q sur P 
douneta la fonction d'offre. 


Le cas général est celui dans lequel iles 2 fonctions se déplacent. 


P 


O1 
O2 


D1 
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Le prix d'équilibre passe de 9,7 Di à O3 ND. Si on régresse Q/P on obtient 
AA qui n'est ni une fonction d'offre, hi une fonction de demande. 
Dans le cas général. quanct Les fonctions d'alfre et de demande se déplacent, une 
régression de Q sur p donne une combinaison ou un Mix de Qu et Qg qu'an peut 
interpréter comme une forme réduite. 
Ceci amène a poser le problème du recourt au modèle structurel à partir du modèle 
réduit. 
Le rectour possible sur Le modèle structurel affre-demande est donc dans une 
certaine mesure analogue au problème de l'identification. 
Cas simple 

Q =23+bp+0R pour autoriser les déplacements de la fonction 

d'offre (11) 

Qu=c+dp+fR f12}pour la fonction de demande. 
Le prix abservé est le prix d'équilibre et à l'équilibre 

Qo = Qu = Q, = quantité d'équilibre 

Î0,-a-br] 


€ 


Daas (11) on calcule R = (là 


Et, ea remplaçant R par sa valeur dans (LA on retrouve une relation encre la 
quantité ei le prix : 


Q, = c+dP+ TT, — a bP] 
e 


a[i-£]-e-#+pla-#] (14) 
€ ê L C2 


= ec-af | pld-t] (LS) 
ef e-f 


Soit 
Q,.=a+BP (16) 


On voit hicn que (1$} apparaît comme une comhinaison des coefficients de (1 Li et 
(12) qui peut être assimilée à une forme réduite. 
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Connaissant & ct B peut-on rttrouver ab, c, d. e, F7? 


En général. il n'est pas possible de remonter à la forme srructurelle, comme on 
peut le voir en comparant { 15) et i L6. 


Par contre, à La foncnon de densande est stable alors f = ( et dans ce cas (15h se 
réduit à Qe = c + dp ec qui correspond bien à La fonction de demande (10). 


Symétriquement si l'offre est stable on a e = À et on retrouve bien l4 fonctiot 
d'oftre Qu =a+bp. {9) 


Le problème de l'identification permet de mettre en évidence un certain nombre de 
ditficultés concértant l'estimation des fancuons d'affre ec de demande : 


-  Onne peut pas estimer directement des fonctions d'offre ou des fonctions de 
demande à partir d'observations sur les prix et sur les quantités échangées. 


— Pour que l'on puisse remonter vers des fonctions d'ottre où de demande 1] faut 
que <es fonctions smient identifiables. ce qui n'est pas le cas du modèle (1H), 
(12) ici. 

— Dans certains cas très particuliers il n'existe pas de comportement d'offre ou 
de comportement de demande. On peut alors légitimement estimer celle des 
loncüons qui reste. 

Par exemple la demande de dépôts à vue ou réglementés (caisse d'épargne), 1 
n'y a pas de comportements d'offre c'esç-à-dire que tes conditions faites par 
l'institution ne changent pas quelle que soit la quantité déposée. Dans çe cas 
on peut estimer des fonctions de demande. Mais on n'aurait pas le droit 
d'estimer dans fes mêmes conditions une fonction de demande de dépôts à 
terme non réglementées puisque là il ÿ à un comportement d'offre de la 
banque qui lixera Les rémunérations en fonction des dénôcs qui lui sont 
proposés. 

Autre exemple, lorsqu'on à une offre rationnée, cette fonction pourra être 
estimée directement. 


C- Autre méthode d'estimation : les doubles moindres carrés (DMC) où 
TSES {Two Stages Least Squares) 


C'est une méthode simple el performante. 
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Cette méthode procède en deux étapes : 

- Dans un premier temps elle consiste à régresser (es variables endogènes 
figurant en variables explicatives dans la forme structurelle sur toutes les 
variables exagènes. Ceci permet d'obtenir des estimateurs de ces mêmes 

ï h# n 
variables endogènes };...v,. 


— Deuxième étape : on remplace les variables endogènes qui figurent en 
variables explicatives dens là tonme struclurclle par leurs valeurs caiculées 
dans l'estimation du modèle équation par équation. 


Cette procédure fournit des estimateurs convergents. 
Exemple 
Soit le modèle suivant : 
ee = bo + by, +hoc,_ +ux (17) 
= ap tv + av, tt, IS) 
=c+l+sg (9 

n à | ‘à variables endogènes «,:1:, 
Ce modèle contient » 
4 variables exogènes Leé,:%, 138 

Ii y a bien un biais de simulranéité car dans 17} et (18) y, Üigure en variable 
explicative. 
L'application des doubles moindres carrés consisit, dans la prétnière étape. à 
estimer : 
HE Foto rar tf3& (201 
ce qui fournit des estimateurs Ÿ5,Ÿ ?2 er Ÿ3 permettant de calculer à chaque 


période Ÿ, d'après (20). 


Puis dans la seconde étape, cette valeur calculée LA remplace ÿ, ce qui conduit à 
estimer : 


D 


124 Écomométrie appliquée 


= hp + $,+bgc, 1 ras (21 


l = 40 +4%-1 +42 Ye tu (22) 


Cette dernière estimation se fait par les moindres carrés. 


En remplaçant yt par ÿ on introduit une certaine dose d'imprécision sur la 
variable explicative mais l'effet qui en résulte sur tes estinateurs n'a pas un 
caractère Systématique. 
Remarque : La méthode des DMC est une application particulière d'une méthode 
plus générale qui est la méthode des variables instrumentales, méthode dont le 
domaine d'application déburde largement le problème de simultanéité. 
La Méthode des variables instrumentales (M.V.I} 
Quand dans une équation l'une au moins des variables explicatives n'est pas 
indépendame du terme d'erreur, l'hypothèse 6 n'est pas remplie et donc les 
Moôiadres Carrés ordinaires ne doivent plus être utilisés. 
SotY=XA+U (23 
X = matrice [X...X,]. La matrice des variables est constituée de ! colonnes 
correspondant aux différentes variables cxplicatives. 
Ei, il existe au moins X; tel que E{X'; U}# 0 
Soit W de format {T, k} un ensemble de K variables observées sur T périodes 
telles que : 
— les variables W soient indépendantes du terme d'erreur U; 
— les variables W sont corrélées avec les variables X ; 
— les variables W sont en nombre au moins égal à celui des variables X, 
E(W'U)=0 
E(W'X)=2Z K>1 
Cette méthode procède également en deux temps : 
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« Première étape 
Elle consiste à régresser les variables X sur W. Ceci permet d'obrenir des 
estiräteurs par les moindres Carrés : 


X;, = WiB+v, 


mu (Ted 


Li] 


et donc B, =(W'W) LU" x; 
ts) 


Pour chacune de ces variables X i on obtient l'estimateur B;carrespondant, 
Soit au total : 


Bo =[8.8t] (24) 


h 
À l'aide de cette matrice on peut calculer les valeurs de X. 


X= Meizu 8 ir (253 


Cela permet d'obtenir une nouvelle matrice de variables explicatives non plus 
observées mais calculées. 


D'après (24) et (25) on a: 
X=WB=W{WW) !WX=MX (26) 
en posant M = WOW Wy lw 

+ Deuxième étape 

On régresse Ÿ sur Æ (et non plus sur X£. 
Y=XC+U (27) 


L'estimateur € des moindres carrés &st donné par : 
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CR(XX)IXY (28) 
soit d'après (26) : 


C=(XMMX)IX MY 
Par conseuction la matrice M st idempotente et donc M'M = M 
En exprimant M en fonction de W' : 


c=ix wow WW XX WOW WW Y. 09.) 


C'est l'estimateur des variables instrumentales dans le cas général. 
Remarque : 
SiK<£:ily à moins d'instruments que de variables X. 


L) 
Donc dans ie calcul de X° il faudra calculer { vecteurs à partir d'un nombre de 
vecteurs indépendants inférieur à {. Les vecteurs obtenus ne sont pas linéairement 


indépendants et donc «x x } n'est plus inversible. 

La condition K > ? renvoit donc à une question de format. 

Dans le cas ou K = 4 ({ étant le nombre de variables de départ) alurs X' W est de 
format (x TE CT x #) soit. (E #} et donc l'inversion peut se faire par blocs à partir 
de {29) : 

C=tW x)! EWWEX MY UIX WOW WI LW'F (30) 

Soit 


|C=(WXTIWY Gi) 


Estimateur simple des variabies instrumentales. 


Dans un programme d'ordinateur il est généralement prévu une procédure 
d'estimation par les variables instrumentales, 1.V, à partir de Jaquélle l'ordinateur 
demande la liste des instruments, c'est-à-dire les variables de W sachant que les 
variables X qui sont cxagènes sont Leur propre instrument. 
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L'estimateur des V.[ est sans biais car en remplaçant Ÿ par XA r L dans 
l'équation (29}, il vient : 


La -] 7 
+ E(C}= ex Wow wy ! wx| X'WCWOWY IW'XA+U) | (32) 
= À +| XWwWwyiWw x] X'WOW WI EL) (34 
E(U}=02 ÆA(C)= À ce qui établis que l'estimatcur çst sans biais 


» 62, - (XX) où 


(e} 

à à RSS 

Étant STUMÉ Par -— 
T-k 

On peut donc calculer La variance des estimateurs par Ja méthode habituelle sur la 
dernière régression. 
Un domaine d'application important est constitué par l'économétrie des 
andicipations rationnelles. En effet lorsqu'on utilise une variable anlicipée cumme 
variäble explicative, celle-ci résulle de l'application de la théarié éconcmique 
pertinente appliquée aux données dispanibles, Les variables résuliant 
d'anticipations rationnelles sont donc, pour cette raison, assimmlables à des 
variahles endogènes. 


Par exemple, on sait que Le taux de change intègre le taux d'inflation anticipé. 
ch=:t(.., ff,..)4u, (44} 
Mais le taux d'inflation anticipé, dans l'optique des anticipations rationneiles, est 


lui-même une fonction du taux de croissance de la masse monétaire, du solde des 
finances publiques, de la pression de la demande, etc. 


ds T-G 


LT] x 


fé = gt. 4 (39 


Dans l'équation (34) la variahle explicative "ne peut plus être indépendante du 
terme d'erreur puisque celui-ci contiens v, d'après (35). 
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Plus généralement, on peut aussi se poser la question de savoir quelles sont, dans 
une équation. les variables explicalives qui sont véritablement exogènes et quelles 
sont celles qui ne IC sont pus. 

Cette question cernmie on f'a vu, est le choix de l'économètre dans un modèle à 
équation simultanée mais peut se discuter dans le cas d'une équation omique. Pour 
cula on peut utiliser Ev rest d'exogénéité de HAUSMAN - WU. 

Soit le modèle standard Y = XA + TC 


Se 


On peut avoir des doutes sur l'exogénéié de ceriaines variables explicatives : 
supposons qu'il y ait m variables véntablemment cxogènes &t donc ({- m} variables 
pour lesguclles 1l y & un deutc : 
soit 

22 : ensemble des (2- m) vatiables douteuses 

Z2 : ensermhle des m autres variahles 


{2 2 


Le test consiste À tester : 

Hg : E(Z'y U=û 

HI:E(Z') Dixû 
On peut réécrire le modèle en faisant apparaître Z1 et Za 
Y F, At ZA y (36) 


IT Mn: 
La régression (36) permer d'ohtenir Zé = #1 


On peut ensuite par la méthode des variables instrumentales calculer des 


# 


estimateurs des variables douteuses, les 71. Ces estimations vont être rajoutées 
dans l'équation {36}, ce qui conduit à estimer 


Y=ZA +2A1+Z BV (37) 
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è 
On obtient Z 5: = $, à partir de certe régression avec 83 < So. 
So : 2 eprésente lu variance résiduelle d s à loi 
Le rapport 7 LÙ Oo, représente la variance résiduelle dans So. suit une loi 
u 


de x2-, (à m-} degrés de liberté). 


Principe du test 
Si les variables douteuses (c'est-à-dire les Z1) sont vraiment exogènes, leur 


La] 


explication par les V.I n'est pas bonne et donc l'introduction de Z1 dans (37) n'a 
pas beaucoup d'effet sur la somme des carrés des résidus. 


Par conséquent S, sera peu inférieur à So. 


Si par contre les Z; sont vraiment endogènes, alors l'infroduction des Z; 
diminuera sensiblement la sornme des carrés des résidus d'uù la règle de décision 
suivante : 


Sy —S 
Lex onaccepte HA 
6 
MS. > X2 = on rejette H@ 
Exemple ; 


Soit une fonction de consommation élémentaire : 

Ce 29 + A1 Fit az Col + Ur (38) 
On peut douter de l'exogénéité de y car le revenu dépend fortement de la 
cansommation, Dans ce cas, les variables explicatives peuvent être réparues 
entre : 

Z1=h 

22 = Cr 
En régressant (38) et on obtient So. 


4 


On estime ensuite Zi c'est-à-dire ÿ, dans ce cas. 


Pour cela on peut utiliser comme instruments 
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. les dépenses gouvemementales ; 
. Les exporiations. 


Soit y=bg+ by gels x+v, (39) 


Ce qui permet de calculer $ ,; Qui est rajoutée dans (38) 


La nouvelle équation est c, — ag ax, +a2cs +83, (49) dont la somme des 
carrés des résidus fournit Sy. 


lei le y? aura un degré de liberté. 


Il existe encore d'autre méthodes qui permettent d'estimer les modèles à équation 
simultanées. Par excmple on peut estimer l'er.sembie des coefficients d'un modèle 
en une seule étape par La méthode du maximum de vraisemblance. 
La méthode est une généralisation de éelle qui s'applique dans le vas d'une 
Équation unique : 
Dans Le modèle à une variable explixative 
Y=ax,rb+u, (41) 
u=y—ax,—b (42) 
d'après l'hypothèse H7 
u, — N(0,62) 
2 
I 4 — FE ty ) 
2e « 
I e H 
oY2ll 
Plur..ur)= Plus)... Plur) puisque ces événements sont 
indépendants d'après #14, En tenant compte de H2 et H3: 


et Pl }= 


y © -trÈt 
e 


Piur..us)= FR (43) 
are 
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NB. : Maximiser la vraisémblance d'avoir observé Ü,.. ü, revient bien à 
minimiser Lu?, d'après l'équation (43}. Ce qui vérifie l'équivalence entre la 
méthode des moindres carrés et celle du maximum de vraisemblance, dans ce cas. 


Cette démarche peut être étendue à un modèle à plusieurs équations. 


Vin 2 Ve an ist bikes = 14, (44) 


ls. | - NES D. 

nt Tee Gaine): BaiXt—e booXmer © By 
Sur ce modèle on peut calculer une fonction de vraisemblance qui traduira la 
probabilité d'observer lu valeur dés résidus des T pérniudes simultanément. 


FA = PC Ts dnqe Bar) 


un T 
L=ITIIRA, surlesT pénudes et les n résidus 


e-lt-l 
D'où 
f T FA #2 
] CFA j —k Es, + E | 
L-{ a) ——|e -* 51 (45) 
V2 # 
5, 


Les csumateurs du maxintim de vraisemblance sont les valeurs des cocfficients a 
et b du modèle (44) qui maxumisent La fonction (4$}. 


Qu les obtient en dérivant la fonction de vraisemblance par rapport aux 
coefficients a et h. 


gL 
da 2 
“dE 
= 0 
am 


L est dune possible d'utiliser la méthode du maximum de vraisemblance pour 
estimer simultanément tous les coefficicnts d'un modèle. 


132 Économétrie appliquée 


# 


La méthode du maximum de vraisemblance à information Jimitée ou 
méthode de la Cowles comission 


C'est une méthode qui procède équation par équation. Il faut bien la distinguer de 
méthodes voisines. 

En particulier, le maximum de vraisemblance à information limitée est différent 
du maximum de vraisemblance appliqué équation par équätion qui n'est autre que 
la méthode des MCO, comme cela vient d'être rappelé, 


Il existe d'autre méthodes qui n'ont pas beaucoup d'intérêt pour les applications 
simples de l'économétnis. 


En théorie, la méthode des MCO étant incorrecte il faut par conséquent recourir à 
l'une des méthodes appropriées, mais 1a comparaison des résultats obtenus par 
différentes méthodes dont les MCO ne confirme pas totalement ce résultat 
théorique. 

Quand on compare les résultats obtenus par ditférentes méthodes on devrait à 
priori trouver les résultats des méthodes correctes d'un cêté {c'est-à-dire assez bien 
groupées) et d'un autre côté le résultat des MCO qui serait alors manifestement 
faux, comme en (a). 


Mco 
{a} 
méthode 
correcte 
Mco 
{b) 
méthodes 
correctes 


Mais on constate en réalité que dans la plupart des cas lé résultat des MCOQ se 
situe au centre tandis que les résultats des méthodes correctes sont assez 
dispersées, comme indiqué en (b}. 
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De toute façon une telle comparaison n'a guère de valeur probante car il est 
impossible de connaître la vraic valeur du coefficient estimé. 


Pour comparer de façon plus pertinente les propriétés des différentes méthodes, on 
utilise la méthode dite de Monte Carlo dans laquelle on se donne d'avance les 
valeurs des coefficients recherchés. On crée ensuite des sérics de variables 
expliquées en donnant des valeurs aux variables explicatives. Enfin on ajoute à 
ces séries un terme d'erreur dont on a spécifié la distribution et dont les valeurs 
sont obtenues par tirage. 


Exempie 
CZ 40 + Yi 
- On sc donne 49 ct ay On üre une série de y 
Ceci permet de calculer une chronique certaine de c,= ay + 45: 
. On spécifie la distribution de u, 
Exemple 
u, — N{0:a2) 
On tire au hasard T valeurs de u, à w. 
On calcule c,= 35 +21yt +ut 
. On régresse 
Ct= âg+AIYt+ Ut 
Ce qui permet d'obtenir les cstimateurs 49 et à que l'on compare avec les 
vraies valeurs de ap et a. 
Ceite comparaison permet d'apprécier les résultats d'une méthode 
d'estimation. 


+ En général quand on l'ait des simulations par cette méthode, un multiplie les 
séries, par exemple en modifiant la distribution du terme d'erreut ou les valeurs 
des variables explicatives et on s'intéresse aux propriétés en moyenne des 
ctimatcurs. 


De même on peut modifier T : la longueur des séries pour analyser leurs 
propriétés asymptotiques : $e rapprochent-ils de leurs vraies valeurs quand on 
augmente le nombre d'observations ? 
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Des études par la méthode de Monte Carlo ont été faites et permeLtent d'avoir un 
jugement au moins partiel sur l'ensemble des méthodes d'estimation. 


Globalement ces études sont assez canvergentes. L'une des premières a été 
réalisée par Basmaun : 


Étude de BASMANN (1958){1} 
Le modèle utilisé est un modèle à trois équations dont l'une est : 
Vi - 2y2 + LS ya + 42), - Oées + lU=u; 


Dans cette équation on sait d'avance que les valeurs vraies sont : 


4 =—2 
a = LS 
a =3 
LOF = -{,6 


y = variables cndogènes 

z = variables exogènes. 

Rasmann à réalisé 200 expériences et a calculé Les valeurs moyennes obtenues 
pour ces cocfficients, Soient ai les estimateurs ohtenus pour le coefficient a. 


tr 
Xi = ta ik 
Méthodes LS 2 «3 4 
d'estimation 
MCO - U,89 1,13 233 - 0.59 
&} DMC -1,84 147 2,95 - 0,65 
MVIL -2.69 1,66 3,52 - 0.60 


On véribe que les estimatcurs des MCG sont biaisés et que leurs valeurs sont 
assez différentes des valcurs vraies des cueflicients de régression. Des trois 


(A) R.L. Basmann "An Experimental Investigation of Some Small Sample properties of 
Generalized Classical Linear Estimations of Structural Equations, Some Prelimipary 
Kesults". GE. Company ; Hanford Laboratories, 1988. 
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groupes d'estimateurs, es MCO sont en moyenne les plus éloignés des valeurs 

vraies des pararnètres. 

Ensuite Basmann s'est intéressé à la dispersion des coefficients «tj autour de leur 
= : 1 = 

moyenne &;.Et pour cela il a calculé 390 21% -&; ÿ. 


Variance des eslimateurs @; : 


LEA {x &X4 4 
MCO 1,91 0,05 1,08 0,46 
DMC 2,33 O,14 2.10 0,83 
MVIL 17,39 1,08 14,97 1,74 


Le tableau montre que les éstimateurs des MCO sont de loin les micux groupés, 
alors que veux des MVIL ont une variance considérible, 


Dans un troisième temps on peut calculer la dispersion autour des valeurs vraies : 


là 2 | dite 2 
700 24 —a;)" = 300 —U;} +(œi-a;) 


& ü3 ü3 ü4 
MCO 2,24 0,19 É,53 0,46 
DMC 2,43 0,14 2,10 0,63 
MYIL 17,89 1.1] 15,25 174 


Les estimateurs dés MCO ne sont pas tellement plus mauvais que les autres et ceci 
cn raison de leur moindre dispergon, (dans trois Cas sut quatre). 

Parmi les deux autres méthodes dites exacles la meilleure de très loin est la 
méthode des DMC. La méthode du MVIL parait effectivement peu able avec 
une énorme dispersion. Les méthodes correctes sont asstz bicn centrées sur les 
valeurs vräies mais avec une dispersion qui peut être importante. 

Les estitnatcurs des MCO sont biaisés mais pas forcément de façon très 
importante ot ils sont bien groupés. Passer des MCO à une méchode appropriée 
c'est passer d'un hiais faible mais certain à un biais important mais pas cerlain. 


Qn peut en donner un illustrahon graphique à partir du problème du cir à l'arc. 
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Si on veut atteindre le centre de la cible. Les méthodes appropriées sont centrées 
et perinctient d'atteindre n'impurte quel point de la cible. 

La méthode des MCO ne permet d'atteindre que le petit cercle. [1 n'est pas centré, 
mais le résultat n’est jamais très éloigné du centre. 


Cr, l'estimation d'une équation ou d'un modèle pour un ensemble de données, se 
fait une fois et non 200 fois. Par conséquent l'estimation s'apparente au tir d'une 
flèche unique vers Ja cible et dans ces conditions les propnétés de dispersion ont 
autant d'importance que les propriétés des valeurs moyennes. 


Les résultats de Basmann sont relatifs à un modèle correctement spécifié. Mais la 
plupart du temps il y aussi des problèmes de colinéanité, d'autocorrélation, … Que 
peut-on dire quand il y a l'un où l'autre de ces problèmes ? 


Les résultats sont un peu différents et pas toujours faciles à interpréter. 


En présence de forte collinéanité, les estimateurs des MCO sont meilleurs que 
ceux des DMC, Alors qu'en présence de faible collinéarité c'est l'inverse. 


À partir d'un enserable d'études, Koutsoyiannis a effectué un classement pour les 
différents cas suivants(L) : 


(D) À. Koutsoyiannis "Theory of Econometrics“ 2° édition, Mac Milan 1981, pp 506-511. 
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Omission d'une Malticolhnéarité 
Cnitère variable Mulicoilinéarité + mission 
DMC DMC DMC 
ds MVIL MVIL MY 
rs MCO MY MCO 
MY MCO MVIL 
MCD MCO DMC 
Variance autour MY MY MCO 
de Ja moyenne DMC DMC MV 
MVIiL HVIL MVIL 
DMC DMC DMC 
Dispersion autour MIT. MY MCO 
de Ja valeur vrac MCO MCO MY 
MY MVIL MVIL 
CONCLUSION 


—  Iln'est finalément pas complètement absurde d'estimer un modèle à équatians 
multiples équation par équation par les MCC. 

En effet, certains, voire la plupart des modèles continuent à être esiimés par 
les MCO bien qu'on sache que c'est mcorrect. 

— Parmi les méthodes théoriquement correctes la seule qui ait un intérêt pracique 
est la méthode des DMC. Quand on traväille sur un modèle à équations 
simultanées il est prudent de faire une estimalion par les DMC et de comparer 
les résultats avec ceux des MCO, au moins pour certaines équations 
importantes. 


Chapitre IIT : L'économétrie sans modèle, les 
séries temporelles 


Dans les deux premiers chapitres l'estimation se l'ail à partir de relations 
structurelles censées traduire de’ mécanismes £conomiques. La modélisation 
consisté à traduire en éguälions la théorie économique. Dans ve chapitre ces 
éléments structurels disparaissent. De ce fait, on ne cherche plus à expliquer mais 
simplement à décrire ou à prévoir. Pour sela on uulise des séries &'esl-à-dire des 
observations d'unc variable au cours du temps. 

L'utilisation de ces séries se généralisé pour tout un ensemble de problèmes pour 
lesquels La modébsation traifhionnelle apparaît peu satistitsunte. 

C’est Lo cas pour des phénomènes complexes das lesguols il y 1 de norubreuxes 
actions Gt ééactions simulranées pour lesquelles il est difficile de faire apparaître 
clairement un enchaînement de causes ei d'effets. (Exemple : cours d'une valeur 
mobilière au jour le jour). Cela peut aussi être le cas de variables extrémement 
volatiles, Les moilèlcs explicatifs sont en général incapables de prévair 
correctement Les points de relournement Ces modèles sant performants À moyen 
rme mas peu pérformants à cout ternx, lc taux de change, qui subit des 
fluctuations as$6Z Courtes et assez amples, fournit dés illustrations de ce type de 
difficultés. 

[ne des causes de mauvaise performance des modèles explicatits vient du Fait 
qu'ils sont beaucoup plus performants pour analyser des tendances que des 
fluctuations. Une grande partie de [eur pouvoir explicarif cs tiré du parallélisme 
des évolutions. Cet argument laisse entendre que certaines cérrélations nt 
correspondent pas véritablement à des relations entre variables, mais peuvent tout 
simplement résulter d'une évolution semblahle dés variables sous des influcnècs 
n'ayant rien à voir avec le problème étudié, Cest La question des corrélations 
fortuités sur laquelle on reviendra dans le chapitre suivant. Un test intéressant #st 
de comparer lé pouvoir explicatif de la même relation macra-économique pour 
des pénodicités annuelle, trimestrielle où mensuelle. 
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Dans tout ces cas, quandt on a hesoin de faire une prévision à caurt terme et quand 
Les enjeux de cette prévision sont économiquement importants on préfère utiliser 
des modèles non explicatifs mais assez performants en prévision qui sont donc des 
modèles de série temporelle, 


€c sont des structures légères qu'on peut facilement récstimer, fe modèle peut 
donc être mis à jour tous les jours. Je cafit de la prévision est relativement faible. 
Un économiètre biun formé avec nn bon logiciel peut assurer au guotidien la 
prévision au Jeur fe JOUr, 

Un des domaines d'application est le domainc des variables financières (taux 
d'intérêts, faux de change, taux de rendement). 


+ La méthode consiste à rechercher dans l'histoire de (a variable des régularités 
susceptibles d'aider à prévoir ses valeurs futures. 


Remarque : On peut auvsi associer les séries temporelles à une équarton 
structurelle : c'est ce qu'on appelle pafois des fonctions de tranxfert où des 
ARMAX et cela permet d'améliorer les qualités en prévision d'un modèle 
structurel, 


Dans ce chapicre, on commencera par présenter une classe particulière de séries 
qui présentent des propriétés intéressantes. On analysera ensuite les prapriétés des 
différentes catégories de modèles de time séries avant de présenter quelques 
éléments perméttant leur estimation. 


SECTION 1 : LES PROPRIÉTÉS DES PROCESSUS 
STATIONNAIRES 


L'utilisation des sérivs temporelles conduit à rechercher des régularités dans les 
valeurs passées de la série, Pour que cette démarche ail un sens pour la prévision, 
1 faut que le processus présente une Certaint stabilité ou un certain degré 
d'invariance au cours du temps. 

+ Si le processus était complètement critiqué on pourrait peut Être tirer quelques 
enseignements du passé mais ces Éléments n'auraicnt aucun intérét prédictif.. 

C'est cctte idée de stabilité ou d'invariance au cours du temps qui est traduite par 
la notion stalisixque de stationnarité. 
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On peut définir la slationnarité d'un processus de série temporeile par Îles 
3 conditions suivantes : 


= LS condition : un processus stochastique de série temporelle est invariant s'il 
a une moyenne constante au cours du temps : 


EGpy=m vtr 
— 2 condition : ce processus a une variance constantt CH = ft = m}| 
PRE D 
soit GC, = Yo Vt 


— _4$ condition : les covariances entre deux observations de y pour deux périodes 
différentes ne dépendent que du nombre de décalages. 


El -mWy, 4 me EG mhtr ee -m)]= 4 
COV (Pipe) CONCT Tea RD TR 


On peur retenir l'image d'un processus régulier qui s'inscrit dans une bande 
horizontale, centrée sur sa moyenne. 


À partir de ces propriélés on peul calculer empiriquement et théoriquement des 
cocfficients d'auto-covarniance Y4 ot des cocfficients d'auto-corrélation qu'on 
désignera par P4 coefficient d'auto-corrélation entre ÿ, et},_y OUY,45 . 


On peut définir des auto covariances qui seront désignées par Y4 
Ye = CON Pr ne) = COVCyeYe rad. 


Si le processus est stationnaire cette covariance ne dépend que du nombre de 
décalages. 


Covariance particulière : yp.=cov(y,r,)= oi 

De même on peut définir des coefficients d'auto-corrélation qu'on désignera 
par Pr = Pr) 

Par application de la formule statistique, 


L CV Yea À 
Or OYi+x 


Si le processus est stationnaire la variance de y. est constante au cours du temps. 


Pk 
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Alors 


COM VAL} Ye 

Ya To 
Tous les cocfficients d'auto-corrélation s'ebtuiennent facilement à partir des 
coufficionts d'auto-covariance. I suffic de les rapporter à%a. 


OH OYs8 = OX = y d'où Py = 


Si un processus si Slationnare. son évolution subit une sorte d'amortissement au 
cours du temps, ce qui Se traduif statistiquement par une décavissance rapide des 
couflicients d'auto-corrélation. 

Pout une Série quelconque, (par exemple le taux d'intérêt du marché monétaire) 
an peut loujours calculer des cocfficients d'auto-corrélation empiriques. 


Ces coefficients peuvent être représentés graphiquerrent en fonction de k : 


Si ie processus est starionnaite ces coefficients d'auto-corrélation diminueront 
assez rapidement. 


Un processus non stationnaire peul avoir Une résonance de la période t à une 
période plus éluignée, de telle sorte que des coelficients d'auto-corrélution élevés 
puissent exister pour des valeurs élevées de k. 


Seuls les processus stationnaires sont intéressants our ils sont stables. Il est 
possible de tester la scationnarité d'un processus c'est-à-dire en fait la décroissance 
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rapide de ses coefficients d'auto-corrélation. C’est la première opération à faire 
lorsqu'on étudie une série, 

Pour cela il existe différents tests. 

Définition 

On sppelle “bruit blanc” (white noise) un proccssus stationnaire de moyenne nulle 
et doni les autocovariances sont nulles à exception de 79. 

Exemple : ÿ,= u, sous ts hypothèses faites sur u, est un bruit blanc. 

Si le processus étudié est un bruit blanc. Bartleti a montré que les coefficients 
d'auto-corrélalion suivent une loi normale de moyennc nulle et d'écart-type 


Ke P-"0 fr 


+ Le test de Bartlert consiste donc à vérifier si les cocfficients d'auto-corrélation 
successifs sont significativernent différents de Q. 


Exemple 


Pour T = 10, tout coellicient inférieur à - 0,2 ou supérieur à + 0,2 sera donc 
réputé significatif à 95 %. 
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+0,2 


F 0,2 


On utilise aussi deux autres tests qui portent sur un énsembie de coefficients 
d'auto-corrélation et non sur les cocfficients pris mdividucllement. Ce sont les 
tests de Rox et Pierce et de Ljung ct Rox. 


Test de Box et Pierce 
Il consiste à calculer la somme des carrés des k premiers cocfficients 
d'aulocorrélation, puis à multiplier cette somme par le nombre d'observations. 


& 
Q= TÈ6E Ce produit suit une loi de y°. 


0 — Fa à K- nombre de paramètres estimés {soil k - p - q dans le cas général). 
T = nombre de périodes 
Le test se formule de la façon suivante : 
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Si e caléulé < ge on conclut Hp, c'est-à-direque les coefficients 
d'autotorrélation ne sont pas sigmifcativement différents de zéro. 


Avantage du test 

1 permet de ne faire qu'un nombre Jimité de tests et resté valable quelle que soit la 
moyenne, à la différence du test de Rartlent. 

Mais ce test est peu performant sur de petits échantillons. On le remplace alurs par 
le test de Ljung et Box. 


Test de Ljung et Box 
: £ 
Q°=T{T+ DAT-KIT* PÈ 
L'interprétation est Ja même que pour le test de Box et Pierce 
L er - 
Q 4 k—p-9 
Exemple : on s'intéresse à la série des taux d'inflation trimestriels mesurés par le 
taux de variation des prix à la consommation, gn France de 19851 à 1953.II. 


Sur les 3(] premiers coefficients d'auto-corrélarion, donc pour k = 30, on calcule 
les tests de stationnatité. 


Box et Pierce = 35,75 

Ljung et Box = 53,39 

Aucun coefficient n'ayant été estimé, ces tests suivent une loi du 7 à 30 degrés 
de liberté. La valeur critique de la table est y? = 43,77. 


On constate que les conclusions divergent puisque pour Box ct Pierce, la 
stationnarité est acceptée, alars qu'elle ne l'est pas pour Ljung ot Bux. Compte 
tenu du petit nombre d'observations, c'est cette dermière conclusion qui doi être 
retenue. 
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Mais, un certain nombre de processus ne sont pas stationnaires. Dans le cas d'un 
processus non stationnaire il faut essayer de le rendre stationnaire. Pour cela, la 
méthode la plus courante consiste à appliquer l'opérateur différence” à La série de 


départ. 


Soit y: la série de dépan, 

Alors w=Yi-ÿe1 = ÂY 

En ütilisant l'opérateur retard L ceci est encore égal à (1-L) yr 

Sur w, (la nouvelle série) on peut à nouveau effectuer les tests de 
stationnarité : 

. W, est stationnaire, il est inutile de continuer. 

. W, m'est pas stationnaire. Il faut continuer la procédure en calculant une 
nouvelle série : 

= AW =Ww,-W, == Liw, =(1-Lÿv, 

De même, on leste la stationnanté de 2, 

Si nécessaire on pourra appliquer une truisième [fois l'opérateur 
“différence”. Thécriquement, un processus qui esl stationnaire après 
l'application de l'opérateur différence reste stationnaire après toute 
nouvelle application de l'opérateur c'est-à-dire que si w, est stationnaire, zx 
ost stationnaire. 


Définition 
Un processus qui est stationnaire aprés une différenciation est dit intégré 
d'ordre |: 


— Si west stationnaire. "y, est intégré d'ordre |" 


— Si w, non stationnaire, mais, 7, stationnaire "y, est intégré d'ordre 2" 


Le degré d'intégration d'une série est donc le nombre de fois qu'il faut 
appliquer l'opérateur "différence" pour rendre le processus stationnaire. 


Remarque ; L'opérateur retard “L' s'utilise comme un opérateur mathématique 
ordinaire : 
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a =Q-LP y, 
=(-2L+12?}y, 
=y,-2Ly +Ly, 
ÆYy = 2 Ver # Yi-2 
Tous les processus temporels ne peuvent pas être rendus stationnaires aussi 
facilement. 


Les processus peuvent être représentés par différents modèles qui ont des 
propriétés particulières. 


SECTION 2 : LES PROPRIÉTÉS DES MODÈLES DE SÉRIES 
TEMPORELLES (TIME SÉRIES) 
Il existe différents madèles dc séries temporelles : 
AR 


ARIMA 
(AR) = modèles auto régressif (Auto Régressive) 
(MA) = moyenne mobile (Moving Average) 
(ARMS4) = combinaison d'un AR et d'un MA 
(ARIMA) = processus non stationnaire aulo-régressif intégré. 
Sous sa forme générale, un praccssus de série temporelle fait dépendre une 
variable de ses propres valeurs passées ainsi que des valeurs passées du terme 
d'erreur. 
Soit en utilisant l'opérateur retard : 
@{L)y,; =a{E)u, +8 (1) 


Ê = terme constant, nul si la série est un bruit blanc. 


À 
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Pour pouvoir estimer économétrigquernent un tel processus il faut d'abord 
caractériser sa structure. En particulier déterminer le nombre de termes qui 
figurent dans les deux polynômes de retard. 


En l'absence de cette précision la structure n€ peut pas être estimée. 

N'y à danc une opération préalable à l'estimation qui consiste à identifier ou fixer 
la structure à estimer. Et, pour fixer cette structure on peut s'aider des propriétés 
des processus élémentaires (MA, AR, ARMA). Ces propriétés servént à 
reconnaître un type de modèle. 


Le modèle MA (q) 
C'est un pracessus qui à de la mémoire sur le termc d'erreur. 
Yr= M + Ut: A1 We + - à q y {2} 


ou de façon équivalente y, = a(L)u, + m 
Sous les hypothèses habituelles concernant les termes d'erreur, ce processus à 
comme moyenne m :E (yr) = m 
Et la variance du terme d'erreur est constante E (u2) = oÈ 
Les termes d'erreut sont eux-mêmes indépendants au cours du temps 
Etuuk)=0sk#U0 
Propriétés permettant de l'identifier et de connaître q : 


Ce sont les propriétés du coefficient d'auto-corrélation du MA (q). Pour les faire 
apparaître, il faut passer par les auto-covariances puisque 


Yo 
Calcul dé Yo : 
E{r, -m} = Eu, — au, —.…. Age 4 y 
= Etu? +au} +. +aQu +,..—24j44 | +... 
+Laqua us -gue-g) (3) 


Gr, les termes d'erreurs satisfont aux hyporhèses habitueiles donc 
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E Qu, u,,4) = 0 par hypothèse (si k # 0}. Ce qui entraîne que tous les doubles 
produits sont nuls. 
Sous l'hypothèse d'homoscédasticité Alu? }= Etui )=... 07 
D'où Yo = 02 (l+af+..+aÿ) (4) 
Pour que cètle variance soit Linie {conditen de la slationnanté} il faut que 
Ltuf+...+as <oe. 
Calcul de y] 
F1 = COV(yYe 1) 
n= Ets, mm) t M] {5} 
= ElCu, — ques 7. 1 Gt 2... ait, 1] (6) 


Puisque les termes d'erreurs sont indépendants au cours du temps il ne reste que 
les termes relatifs au même décalage. 


= c?[-a +87 +agag | (2) 


Calcul de 

pet 

+ Soit le proccssus simple MA {1} 
M=mtne-anm.| (9) 

* E(y -m} = variance de y, 


= Eee) ]<ott+a)= 7 (0 
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eh= EG -"Xv1-m)] (D 
D =) és {12} 
d'où Yi = ElGe 8-1) m8 2)] (13) 


soit Y1=-402. (14) 
== (15) 
1 


+ Y2= EG -mXy 2 -"m)] (16) 
Via =Mitu om; nr (17) 
Ya = EG au Mes qu, 3)]=0 (18) 
or, p=2 donc p3=0 (19) 
Yo 


Donc, pour un processus MA (1), les coefficients d'autc-corrélation #4 sont nuls 
pour tout k > k, 


De façon plus générale, pour un processus MA (q), les coefficients d'auto- 
corrélation sont nuls pour tout décalage k > q. 


Soit un MA (q} avec k > q 

Ya = E[(u, — ht] re Ggté-g Xi -4 “A4; Aytr-s-e)| (20) 

Pour k > q il n'y a aucun élément commun entre les deux parenthèses et donc 
Yx =0 c'est-à-dire 9 k= 0. 


Înversement, quand ta fonction d'auto-corrélation chute brutalement au bout de g 
décalages c'est l'indice que le processus pourrait bien être un MA (q). 
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MA{2) 


Un processus MA pose un problème d'estimation car il contient un terme d'erreur 
composé qui ne peut pas satisfaire les hypothèses habituelles de l'écunométrie, 
même si individuellement les u, satisfont ces propriélés. 

Supposons maintenant qu'on connaisse une fonction inverse de a (L) soit [atL}[!, 
alors y, (a(L)7!'= u, est un processus qui présente des conditions d'estimation plus 
favorable car il n'y a plus que u, comme terme d'erreur. C'est l'opération 
d'inversion d'un processus MA qui peut être nécessaire pour cstimer cc P'OCCssus. 
Exemple 

Soit le processus MA (1) supposé être un bruit blanc. 

Ye=u au. (21) 

soit y, ={l-aLju, {22) 


Pere 
ou Ta u, (23) 


On fait ainsi apparaître 


l —4L 


* Si a<l 


1 peut s'interpréter comme la somme des termes d'une 
—&jl 


progression géométrique : 


=t+aL+afL2+4..+afL" (L étant l'opérateur décalage) 


1-aL 
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d'où 
1 


l-a 


: =, ({+aL+..+ai if} (24) 


Avec Ja définition de l'opérateur retard, cela peut encore s'écrire : 
VAS TIURE +@f Jin = 4, (25) 
Le processus (25) est alors l'inverse du processus (21). 


On voit ainsi que l'inverse d'un processus MA est un processus AR de durée 
infinie, 

Ce calcul fait aussi apparaître la condition d'inversibilité d'un MA : cette cunditiun 
a été exprimée ici sous la forme a, <f. 

Si on reprend le polynôme de retard 1- 4, L cette condition a été exprimée sur a : 


à < 1 pour que le processus soit inversible. Mais on peut aussi l'exprimer par 
L>1. 


[-a .=0— ra 
a 


L est ia racine ici unique, du polynome de retard 1 - a, LE. 
La condition d'inversibilité peut aussi s'exprimer par L > 1. 


Dans un MA (1) la condition d'inversibilité peut s'exprimer alternativement sous 
la condition a] <1 ou bien racinc en L du polynôme de retard > }, 

Cette condition est facilement généralisable et fournit La condition d'inversibilité 
dans le cas général. 


Un processus MA (a) est inversihle sj Jes racines du polynôme de retard en L sont 
situées en dehors du cercle unité. Cette condition rend compte des cas de racines 
réelles et dé racines complexes. 
Exemple 

MAG2)  yf=u,-ai Ur - azur 

soit y, u(l-81 L-a L?) Le polynome de retardest 1 - a1L - a2L? 
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En posant l—-«a£-as F2 20, on fait apparaître 2 racines L' et L” qui 
peuvent être réelles où complexes. 
Pour que le MA (2} soit inversible il faut que ces racines sortent de 
module > 1. 
Remarque : Cette question à été abordée dans le chapitre 11 avec les propriétés 
dynamiques des modèles. La comparaison enire les méthodes suivies aux chapitre 
I et HT permer de mieux situer les deux démarches sur un exertple. 


Chapitre LL Palynôme de r 


Ye Tr 55 


Ve Ye 3 — A2Yr.2 


ÿ= rt polynôme él a La) 
Pure 0 caractéristique PT 
A = a? +4a) A = a? +4a» 
._a+Ya p-ntva 
2 —24t 
PL Ÿs racines À VA 
2 -2@; 


On vérifie aisément que r'L" = r"L'e 1, Ilest donc logique de retrouver sur les 
racines en 1 des conditions de stationnarié qui soient inverses des Conditions 
sur r. 
Ou voit que la stationnarué peut s'aborder de deux façons : 

. de façon empirique: 

- de Façon mathématique. 
Il se peut aussi que le polynôme de retard admette des solutions où L= 1. 
On dira alors que “lc polynôme en L admet une ou plusieurs racines unitaires". Le 
prucessus ne peut alors plus être stationnaire ét céla donne une troisième façon 
d'aborder 1e problème de lu stationnarité : en utilisant des tests de racine unitaire. 
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B- Propriété du modèle auto régressif AR(p}. 

Comme précédemment, ces propriétés serviront à reconnaître de tels processus. 

Soit le mudèle AR{1}) : 

Me +Ë+w (26) 

Ea valeur moyenne de ce processus est donnée quand y, = #1 V,, Dans ce cas : 

m=am+ê (27) 

d'où 

ê 


mm = 
1 


Si on pose & = 0 alors m = Ô (Hypothèse de bruit blanc}. 
ro = EO)= AGayi+u ]e Efoty, 2+u? +2ax,.ju,] (28) 
Er, 2) = 70 (cl.hypothèsedestationnarité) 
Ye 1 né peut pas contenir u, d'après(26) et donc E(y,_çu, }= 0 
d'où 


à 
= x 
Yo = ? 
ns EG 1) 


= Efay,-; +u, )x1] 
= Eleu?. +1] 
476 

= Yo 


LA 
A=—=a, (29) 
( y 1 


ü 


L'économétrie sans modèle, les séries temporelles 
V2 = EC:-2) 
= Ef(air,+4)32] (30) 
= ir 
Or, = & To 
donc 72 = 270 
P2= ai 
a = Ely7s-3) 
= Ef(ay1+4)]»-3 (GD 
Soit Ya =@1}2 


Or, Y3 = a? Yo 
Donc }3=ai}9 


ar 
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Pour un processus AR (1) les coefficicors d'auto-corrélation décruissent de façon 


gévrmétriqué à partir du prérmier. 
Ce résultat se généralise pour un ÀR {p). 


+ Dans un AR (p} les coefficients d'auto-corrélation sant géométriquement 


décroissants au-delà de k = p. 


Symétriquement un processus qui fait apparaître une décroissance géométrique 
des cocfficients d'auto-corrélalion au-delà de k = p a de honnes chances d'être un 


AR (p). 
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— Le processus représenté pourrait être un AR {5}. 
Un AR{?) s'écrit : 

Yi = A} +22 +4 (32) 

(En supposant qu'il s'agit toujours d'un bruit blanc) 
ô 


Ê 8 #0 alors m= 
TETE 


— d'où l'intérêt de poser 8 = 0 


E(y)= Yo = Foy + ay, 0 +43) M +724 Etyu) UD 
or, E(yu)= E[(ay +@y, 2+uKu)]= 05 (34) 


Donc 
Yo=Mi+0Y»+08 (35) 
= EC») 
= ElCaye 1 +623. +4)3,1] (36) 
= Yo + 2h 
D'où 


Yr= Yo +71 (47) 


P = L- a+e:P (38) 
Yo 
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LCA Î 
= 39 
Pi I- a (39) | 
Va = Allons + @2v 2 + y-2)]e Gin < Yo 
soit 
fr Mr Ya 


Pr = CP + (40) 
2 
4 (4 
2 


ou P2=- 


Remarque : Dans un proccssus ÀR (2) il v a deux coefficients à calculer 
(&1 ét 2). Les équations {39} et (41} constituent un système à deux 
inconnues (œ, el œ3). Connaissant les valeurs numériques de 4, ct p3 ces 
équations pcrmettent Lhéoriquement de calculer & ot &. 

Ces équations sont ditcs équations de Yule-Walkor permettant en effet de 
calculer les coefficients du modèle à partir des coctficients d'auto- 
corrélation. 


Toutefois pour utiliser ces équations il faut savoir que lé processus est un 
processus AR (2). Sinon il n'y a aucune raison de s'intéresser au système des 
équations (39} et (41). 


F3 = El(aixs-1 +@2Y-2 tb 52] =&Y1+07; (42) 
P3 = Gifs + GP; 
Et les cocffcients successifs seront obtenus en multipliant par Œ, et C2. 
Fa = EG + av tu 4]= &73 +02 
Pa = }sta 
Pr = Mir + A2Pa-2 (43) | 


Four AR {2} au delà de k = 2 on retombe donc dans Ja situation de décroissance 
géométrique des coefficients d'auto-corrélation. 


Ce résultat se généralise pour tout AR(p) pour k > p 
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C- Le proctssus ARMA (P, Q) 
C'est la combinaison d'un AR (p} et d'un MA (q}. 
&{L)y, = 8+a(Llu, (44) 
Il contient un polynôme de retard sur y et un polynôme de retard sur u avec pet q 
éléments : y, =@ÿ,-1 ++ die pt Ô+u, — LE Poe LP EE (45) 
Comme pour les modèies élémentaires, certaines propriétés permettent de 
reconnaître un processus ARMA. 
Soit un ARMA (1.1). I] s'écrit : 
M = Gi +4, ain, 1 (Eypothèse de bruit blanc) (46) 
Condition de stationnarité du procéssus ARMA. 


De façon générale un processus ARMA est stationnaire si sa partic AR st cile- 
même stationnaire. La partie MA d'un processus ARMA n'est pas susceptible, 
normalement, d'entraîner des évolutions non stationnaire. 


Dans le cas d'un ARMA {lion a: 
1- &,£ = 0 et la condition de stationnarité est donc : 
L>1 
_ VITE 2 
Yo = Ely; )=Ff(a y, +441) (47) 
L Elair +4 +afuf +20 = 2j y, — Zee] 
Sous les hypothèses habituelles 
Yo= &?Yo+02+a/o2-2aa0 (48) 
En effet 
Feet = Mide-c +4. 17 44,;-2 € donc 
; 2 
E(ys48.1}5 où 
d'où ÿ9 = a?(1+ a? - 24,4 )+ {Yo 
Soit 
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__@/{l+u? -2aa) 


1-@? 


Yo (49) 


A =E(ayru-au 1); 1] (56 
n=ar-ac (SD 
Ya = EGen ve + ex, |] 
donc Y1=<&@71 (52) 
h3= EK ŒiVe-i +4 — GW; be] 
Y3=0@Y2 (53) 
ou 
7" af té 
À partir des équations {50} à 153) on calcule aisément les coefficients d'auto- 
corrélation : 


: 2 
PP EE (LUE EE 7") 
To. L+ai 2 

tpm am (55) | 

L LL RON EUR 


paleain Go | 
Ya | 


On constate que dans un processus ARMA {p, q}, la partie MA disparaî lorsque 
k> q. À partir de k= q, le processus ARMA aura ics mêmes propriétés qu'un 
processus AR simple et donc on retombe sur la décroissance géométrique des 
coefficients d'auto-corrélation. 1 n'est pas toujours facile pour cette raison de 
distinguer un processus ARMA d'un processus AR de même ordre sur p. 

C'est pourquoi il est intéressant d'utiliser aussi nn autre indivateut pour identifier 
les processus ARMA 


160 Économétrie appliquée 


Exemple : si en k = 511 y à décroissance géométrique cela peut être AR (5) mais 
aussi un ARMA (5,3). (Le décrocherment sur k = 3 n'est pas toujours évident}. 
D- Le modèle ARIMA (p, d, q) 
C'est un modèle non stationnaire mais qu'on peut rendre stationnaire en passant 
aux différences d fois, C'est par conséquent un processus intégré d'ordre d, 
Exemple : ARIMA (4, 1, 3) 
On le rendra stationnaire en le différenciant une fois. 
Si yt= ARIMA (4, 1,3) 
{1- Lir=w=Ay = ARMA (43) 
Par conséquent ARIMA (p, 0.9} = ARMA (p,q) 


Jusqu'ici on s'est appuyé sur La fonction d'auto-corrélation, appelée aussi fonction 
d'auta-corrélation totade, pour mettre en évidence certaines propriétés des 
différents modèles étudiés. Mais ce critère peut être insulfisant ct on ulilise aussi 
unc autre fonction qui est Ja fonction d'autae-corrélation partielle, 


E- Fonction d'auto-corrélation partielle 
On s'intéresse uniquement au processus AR (p). 
Suit un processus AR (p} qu'on va chercher à identifier . Cela revient à détérminer 
b- 
«Sip= 1 alors le processus va s'écrire : 
= Mer té 
donc Y1 = A3 1 0 Xy:-1)]= Yo 
d'où p, =, comme on l'a vu dans F équation (29) 
Si le processus est un processus AR (1) il n'y a qu'un seul cocfficient à calculer, 
c'est Œ et il vaut p. 
+ Si le processus est un AR (2) : 


Yi Mes + Mrÿs-2 tu, Et 
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P2 = af + (40} 


Si le processus est un AR(2) Les équations (39) et (40} permettent de catcuter ct, et 
œ (Équations de YULE WALKER). 


+ Si le processus est un AR (3) 

Ve = Aiÿsni + MaYs-s + Maÿi_z +4, 

On obtient alors trois équations permettant de calculer &!, æ2, 3. 
Plus généralement pour un AR(P} qui s'écrit : 

Ye = Yet pi D + Ur 

Yo = Ely, (CAPES, #,]= af + AY AYpt ci (57) 
hd A (Ye Cpdip + ue 1] = fa + Ye dp7o1 (58) 


6 a PET CAE +u,)] FM +@2Y2.. 0 Ype2 (59) 


Ya = ER CŒ it pop +4 )] Sie + os. CpY,-e (60) 
Soient les coefficients d'auta-corrélations 9,...9, obtenus en divisant les 
covatiances par 79 : 

Me +0pir. ap, 1 (61} 

Pr ii + ÆPat.. UP (62) 


Pr = Pr 1 + Copa te LD pt (63) 


On appelle coefficient d'autocorrélation partiels les valeurs de &, 2... @y que 
l'an trouve en posant successivement p = 1,2, 4, ...k dans le système des équations 
(61) à (63) qui sont les équations de Yule-Walker, 
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Cette fonction d'auto-corrélation parlielle pout aussi être calculée pour les 
processus ARMA et elle présente certaines propriétés qui peuvent servir à 
l'identification d'un processus de série tempurcile. 

C'est un critère qui viendra compléter le critère dé l'auto-cotrélation totale, Dans 
la réalité, le choix est facilité par Le l'ait qu'on n'utilise pas de processus trop long, 
pet q < 5, sauf en cas de saisunnalité où des décalages de 4 où 12 périodes 
peuvent être névessaires. 

Le tableau ci-dessous, tiré de Cuthbertson, Hall et Taylor) } donne quelques 
éléments permettant d'identiGicr kes processus simples 

Critères d'identification à partir des propriétés des fonctions d'auto- 

corrélation totale et partielle 


à Fonction d'auto-corrélation 


| partielle (ag) 
MA (1) »  Décraissance exponentielle 


ARE) ‘+ Picpourk=1 


Fonction d'auto-corrélation 
totale (9, ) 


*  Unpic pour k=] 


* Une décroissance 


exponentielle 
*  Unpicpourk= I! suivi par ARMA + Un pic pour k = 1 puis 
une décroissance 


! 
{L1} | décroissance expoñéntielle 
exponentielle 
»s Deux pis :k=]letk=2 
suivi par décroissance 
exponentielle 


ARMA +  Unpicpourk= 1! puis 
{1,2} décroissance exponentielle 


ÀARMA ,+  Deuxpics pourk=1letk=2 
{2,1) | puis décroissance 


+. Unpic pour k= 1 puis 
décroissance cxponentiel{e 


| exponenticile 
* Deux pitsk=letk=2suivi| ARMA ‘* Deuxpicspourk=lek=2 
d'une décroissance (3,2) suivis par décroissance 
cxponenticlle exponentielle 


L'appréciation de ces critères suppose une certainc expérience. 


{D x. Cuthbenson, S. Hall, M. Taylor. “Applied Econometric Techniques". Philip Allan 
1992. 
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Au-delà de k = 2 Fidentitication procède à partir des mêmes principes. 

Les critères ne Sont pas très précis el il existé en général plusieurs modèles 
possibles avec différentes valeurs de p et de g pour une série chronologique 
donnée. 

Il faut estimer les différents processus possibles et vérifier à posternion si tes 
processus convicanent, Adtreméent dit il ÿ à un contrôlé à posterion possible sur 
Les choix qui ont été faits. 


En effet si le pracessus utilisé est correct, les résidus de la séne c'est-à-dire les à, 
ne contiennent plus aucune autocorrélation significative car lout a été absorbé par 
le modèle. 


[1 faut donc tester les coefficients d'aulo-corrélation avec les tests habituels 
(Bartictt, Bux et Picrec...) appliqués à La série des résidus. 


Si dans la série des résidus on observe des auto-corrélations significatives il faut 
alors rejeter le modèle. Parmi les modèles possibles, certains ne passeront pas ce 
test et 1] restera un nombre assez limité de modèles acccptables. 

Toutefois il est exceptionne] qu'il ne reste qu'un seul bon modèle, Pour départager 


les bons modèles on peut aussi examiner ce que ces modèles donnent sur La 
période récente. 


Box et Jenkins ont proposé la méthode suivante qui résume les différentes étapes 
de la démarche : 
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Détermination de p et q 


Estimation des processus (en généra] un ARMA ({p, q)) 


Contrôle sur les résidus 
ACCEPTATION , REJET : le modèle ne 
(Le modèle passe le test convient pas et donc 
sur l'autocorrélation des la procédure 
résidus} recommence 


SECTION 3 : L'ESTIMATION DES PROCESSUS ARMA 
Un processus ARMA est en général une double distribution 
&(L}y, = ô + a(Lja, 
+ Quand q = 0, üä se réduit à un AR (p} et ne pose pas de problème d'estimation 
particulier. 
En effet, sous les hypothèses habituelles un AR (p} à un bon terme d'erreur et peut 
être cstimé directement. 
+Siq#0:c'est un ARMA (p, q}. 
Si le terrne en u, possède les propriétés habituelles aiors le polynôme de retard en 
u, ne peut pas Les avoir. 
L'estimation d'un ARMA (p, à) suppose qu'on puisse se ramener à une équalion 
dans taquelle on isole à droite u.. Pour cela on suppose qu'il existe [a(L}]'! le 
polynôme inverse du polynôme à (L). 
L'application de cette transformation sur ia partie gauche de l'équation ramène à 
uné formulation du type : 
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[atLiT 'o(Liy, =u 

C'est donc sous cette dernière forme qu'on peut stimer un modèle ARMA. 
Comme on peut le constater, l'équation ainsi obtenue fera apparaître de 
nombreuses non-linéarités et ces non-linéarités rendront impossible l'utilisation 
d’une méthode d'estimation simpie comme par exemple, les MCO. 

On doit danc recourir à des méthodes itératives donc des méthodes dans lesquelles 


on se donne des valeurs d'amorçage qui néuvent étre tirées dés conditions de Yule 
Walker. 


Dans certains cas, ces méthodes ditération sont cunvergentes. Il peut arriver 
qu'elles ne Le soient pris. 

Îl peut amiver aussi que la convergence se fasse sur des solutions locales et qu'en 
fait le modèle admette d'autres solutitms pour des valeurs d'amorçage différentes. 
La piupart des logiciels d'éconnmétrie ont dés méthodes d'estimation des 
processus de séries Lletnparciles, 


Chapitre IV : Introduction à la coïntégration 


La plupart des propriétés statistiques des méthodes d'estimation s'appliguent à des 
variables stationnaires, C'est-à-dire qu'elles ne Sont pas valables pour n'importe 
quel type de données. 

Oc on applique indifféremment ces méthodes d'estimation à des variables 
stationnaires et à des variables nun stationnairés. D'où Ja possibilité qué ces 
propriétés statistiques ne soient plus valables pour des variables non stationnaires. 


Cet aspect de la méthode économétrique a été ignoré pendunt longtemps avant de 
ressurgir au Cours des dèrtüères années. Li pose des problèmes redoutables dans la 
mesure où la plupart des séries représentant des variables économiques sont 
affectées par une évolution de long terme. C'est par exemple le cas de 
pratiquement toutes les variables macro-écunomiques, y compris des variables 
comme ie chômage ou Je taux d'inflation sur 20 à 30 ans. C'est aussi vrai pour les 
donnécs de beaucoup d'entreprises ou d'unités individuelles. En fait, {2s séries 
dent on peut logiquement penser à prion qu'elles devraicnt remplir les conditions 
dé la stationnanité sont peu nombreuses : les taux d'intérêt, certains taux de change 
(Fdollar mais pas F/Dm par exemple). fes prix de certaines matières premières 
com l'or ou l'argent, 


On s'est intéressé au cas patticulièr des variables avec racine unitaire. En cifet, le 
modèle avec racine unitaire constitue un ças limite puisqu'il n'est plus stationnaire 
au sens du chapitre TT. 

(y = Ye + 
est un modèle à racine unitaire (L= 1} 

G-Liy=u 

— ]-L=0 
Remarque : Un modèle de ce ivpe s'appelle une “marche tu hasard” (“Random 
Walk") 


Ce n'est pas un processus stationnaire : 


Ye =Yr2 tuer (2) 
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Yi=Yt2tüc+uer (3) 


La variance ne peut donc pas être constante, comme on le constate en poursuivant 
le développement de l'équation (3). Ce n'est donc pas un processus stationnaire. 


GRANGER ET NEWBOLD (1974) ont mancré à l'aide de simulations empiriques 
que la distribution du t de Student dans un modèle simple du type y, = ax, +b+u, 
en présence de racine unitaire ne peut plus être interprétée selon la loi de Student 
Fisher. 

C'est-à-dire que si x, et y, peuvent être mises dans la forme d’une AR{1} avec 
racine vüitaire, comme dans l'équation (1), la distribution n'est plus standard. 


De plus, entre des séries admettant des racines unitaires il existe des corrélations 
fortuites ("Spurious Cuorrelation") qui n'ont aucunc signification véritable et qui 
n'ont pas de véritable fondement. 

L'éconarètre peut être amené à conclure à l'existence de relations qui en réalité 
n'existent pas. 

Granger et Newbold ont montré que cela est particulièrement probable lorsque la 
régression comporte un R° élevé et un DW faible. 

En 1986, PCR Phillips a confirmé de lagon mathématique ces résultats et a 
également mantré qu'en présence de racine unitaire fc D'W était aussi biaisé vers 
la valeur zéro. De ce fait la question de savoir s'il élait légitime de continuer à 
travailler sur des séries non stationnaires s'est trouvée posée, En toute hypothèse 
la première chose à faire consiste à vérifier si cette prapridté est remplie. C'est à 
cela que servent les tests de raçine unitaire. 


SECTION 1 : LES TESTS DE RACINE UNITAIRE 


Pour vérifier la stationnarité des séries, il faut pratiquer des tests de stationnarité 
ou des tests de racine unitaire. Dans l'équation {4} 


(4) yt= ay +u 


on va tester : 
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& 


H&:[a]=1 
| H1:|al <1 
Pour des raisons de commodité on transformne l'équation (4) de façon à pouvoir 
tester par rapport à la valeur zéro. 
Yr-ytista-Dyit+u ($} 
Ayr = fa-l) ve +u {6} 
On teste alors : 
Hé:a -1=0 
Hl:a-i<0 
Pour ce test on utilise le test de Studeur mais sachant que la dismihution du test de 
Student sous l'hypothèse HS n'est plus standard, 
Sous H@ cette distribution a Été tabulée par Dickey et Fuller d'où ie nom de test 
de Dickey-Fuller (1976) qui a été donné au plus courani des tests de racine 
anitaire. 
Il'existe une version dans laguelle on introduit dans (#} les accroissements décalés 
de la variable expliquée 
dy =(e- y + Aÿ,_1+...+Ap,, +8 (6) 
connue sous Je nom de Augmented Dickey Fuller (ADE}, qui est une version plus 
générale du test. 
Les valeurs critiques du test se tisent dans la tahle de Dickev-Fulicr. 


Supposons que l'on ait une équation avec 2 variables explicatives et 
25 observations, 

T=2s 

î = 2 = nombre de variables explicatives 


Les valeurs critiques de la table de Dickey-Fuller diffèrent sensiblement de celle 
de Student. 


Seuil unilatéral de Student : - 1,71 


t>- 1,71 on conclur HO c'est-à-dire a = 1 
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————— 


t<-1,71 — HI c'est-à-dire qu'il n'y a pas de racine uaitaire 
Dans la table de D & F le seuil correspondant vaut - 3, 
1>- 3 — H@ c'est-à-dire racine unitaire 
t<-3 — Hi c'est-à-dire pas de racine unitaire 
J] existe donc une différence substantielle entre les deux tests. 
Si la valeur calculée du test est par exemple -2, il existe deux conclusions 
opposées sclon la table utilisée. d'où l'importance de faire k: bon choix. 
Exémplé : on s'intéresse à la variable log PIB pour la France en données 
trimestrielles sur la périvde 1980.1-1992.1V, 
Le test pratiqué est le test de Dickey-Fuller : 
DF =-0.416 
Valeur critique à 5 % = -2,9{4 
La variable log PIB n'est donc pas stationnære. 
En prenant Les différences premières de lug PIB, Le fést danne Le résultat suivant : 
DF = -5,964 
Valeur critique à 5 % : -2,915 


Les différences premières sont stationnaires. La variable log PIB est donc imégrée 
d'ordre I. 


Le test de racine unitaire (test standard} s'intérprète comme un tést de 
staüonnarité. En effet, si une série possède au moins Wie racine unitaite elle ne 
peut pas être stationnaire tandis qu'une série qui n'a pas de racine unitaire peut être 
stationnaire. 

On peut donc utiliser le test de Dickey-Fuller pour tester la stationnarité d'un 
processus temporel au sèns du chapitre précédent, (ARMA où ARIMA) au livu du 
test de BOX et PIERCE. 


Lés deux tests posent en fait le même problème. 


Il existe une autre façon de tester une racine unitaire 4 partir du test de Durbin 
Watson pratiqué sur y, + yes + 0, C'est le test CRDW (Coïntégration Régression 
DW), méthode proposée par Sargan & Bhargava (83). C'est un test moins usuel] 
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que le test de Dickey-Fuller dont les valeurs critiques ont été obtenucs par une 
rclabulation de la table de DW. 


Le test de Dickey-Fuler se pratique aussi sur les résidus dus équations. Quand il 
est appliqué à des résidus, il porte le nom de test de ENGLE et GRANGER 
{ISRT). 


SECTION 2 : LA COÏNTÉGRATION 
Face au problème qui a été rappelé au début de ce chapilre, la théorie de la 
coïntégration permet de préciser les conditions dans Jesquelles il est légitime de 
travailler sur des séries non stationnaires. 

Définition 
Soit deux séries yr et xx toutes les deux intégrées d'ordre d: 


É = 1} 

Va 

S'il existe un vecteur {@, À) Æ Ü tel que z, = x, + Br, 

soit intégré d'ordre d - bh avec o<b<d, alors on dira que x, et y, sont coïntégrées 


d'ordre (d, b) et (x, Bi s'appelle vecteur cuïntégrant. 
Par exemple, avec deux sériés qui sont intégrées d'ordre 1 si une combinaison 
linéaire de ces deux séries est par exemple, intégrée d'ordre O on dira que ces deux 
variäblés sont cv intégrées d'ordre (1, 1). 
Concrètement si on prend © = à y + Cu 
si une combinaison linéaire des variahles ©, ct y, ust stationnaire elles sont 
coïntéprées. 
Dans le cas le plus usuel où les variables sont intégrées d'ordre 1 avec deux 
variables il existe deux relations de coïntégration qui sont réciproques l'une de 
l'autre. 

Y,=ax,+u (|) 

x =ay;+v (2) 
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Une combinaison particulière de x, et y, par rapport à La relation (1} peut être 
a, =}, -ax,. Le vecteur coïniégrant &st alors 1. - a), On peut tester lé fait que y, 
- ax, - b cst intégrée d'ordre O en pratiquant un test de racine unitaire sur la série 
des résidus. C'est le test de Engle et Granger. 
Si par exemple on travaille sur des variables intégrées d'ordre |, après régression 
le rest de racine unitaire pratiqué sur les résidus môntrerä si ceux-ci sont 
intégrés d'ordre | : il v a présence de racine unitaire el non coïntégration. 
intégrés d'ordre 9 : il n'y a pas de racine unitaire — Jes deux variables 
sont coïntégrécs. 
Ce résultat ou cette définition de (a coïntégration pour deux variables se généralise 
à un nombre quelconque de sérics, mais le secteur coïntégrant n'est plus urrique. 
Exemple : En données trimestrielles, les variables suivantes sont intégrées d'ordre 
1, en France sur la période 1980.1-19S3.IV : 
- log PIB 
- log M; 
- taux de chômage 
- taux d'inflation 
- coût d'apportunité mesuré par la différence de rémunération entre les bons du 
trésor à 5 ans et le taux des livrets À, 
Un test de Engle et Granger pratiqué sur une relation entre ces variables donne le 
résultat suivant : 
. valeur calculée : -2,567 
. valeur cnüiqué à : -4,673 
Par conséquent Les variables né sont pas coïntégrées, en particulier le log de la 
masse monétaire M; et celui du PIB nc suivent pas des évolutions parallèles sur Hi 
période, 


SECTION 3 : PROPRIÉTÉS DES RELATIONS DE COÏNTÉGRATION 


Quand des variables sont coïntégrées on retronve certaines propriétés habituelles 
de la méthode d'estimation par les MCO. En fait les estimateurs sont qualifiés de 


introduction à la cointégration 173 


super canverpents c'esr-à-dire qu'ils convecgent plus vite que les MCO sur 
données stationnaires, vers les valeurs vraies des paramètres. 


Cela semble impliquer que dans un cas de coïntégration on puisse fégitimernent 
utiliser le test de Student normal. 


Ea fan, l'imerprétation est un peu plus délicate. 
Dans le modèle 

(D se =ao+aix, +4, 

avec x,=x, 1+v, (2) 


si les termes d'erreur u, et v, sont normaux, de moyenne nulle et indépendants, la 
distribution des tests de Student et de l'isher est bien la distribution standard. 


Lorsqu'il existe des corrélations non nulles entre u, et v,, ces distributions ne sont 
plus standards, mais il est possible de les corriger ca spécifiant la relation existant 
Cntre U et V. 


Pour cela, on suppose que u, dépend de *, pour des décalages de -k à +k 


k 
(3 U — de ét 


R=—k 
et que les comélations sont nulles pour des décalages supérieurs à k, c'est-à-dire si 
l> 4. 


Eu reportant ectte expression dans (1}: 


k 
Yi = 69 + X, + )HA "Tee 
—# 
soil puisque d'après (2) 
Àx,=v, 


d') y,=4 +4 x, +Y rar +e, 
#“ 
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Le nouveau terme d'erreur, e, de la régression (1°} a été purgé de toute corrélation 
avec v. En effec si u, dans {1)n'est pas indépendant de v,, e, est hien indépendant 
de r, d'après (3}. 


La correction proposée par Saikhonen! 1 (91) consiste à multiplier respectivement 
+ 


les F de Fisher et Les t de Student ohtenus sur La régression {1°} par le rapport Eu 
À 
: Te 
pour Fou us pour ? 
4 
: x : l à 
o?étant obtenu à partir de (1°} par l’estimateur habituel] ——56 où lest le 


T-1 
nombre de paramètres éstimés dans {1”}. 


À est obtenu à partit d'un modèle AR appliqué à &,. 


(4) ê, = Œê, LD pËge p +£ 


ñ 
avec où = 2, Ée d’après (4}. 
1-p 

On vérifie que si les termes d'erreur ui et v, sont bien indépendants, y, =0 Van 
et donc 0, se confond avec e, qui à sont tour ost égal à €,. Dans cs vas le facteur de 
correction est égal à 1, ce qui impliqué bien la possibilité d'utiliser Jes tests 
standards de Student ou de Fishér, 

Dans le cas (rès courant où les variables sont intégrées d'ordre 1, le fail qu'elles 
soient coïntégrées implique que leurs évolutions au cours du temps soient 
parallèles c'est-à-dire que l'écart entre les deux séries soit stationnaire, C'est ce que 
traduit ]3 stationnarité des résidus. 


(6) Saikhonen, Peruti : Ascymplotically Efficient Esthnation of Cointegration Regression"” 
Econometric Theory 1991, p 1-21. 
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Qn peut montrer qu'il y a alors équivalence entre équations de coïntégration ét une 
catégorie particulière de modèle structurel appelé modèic à correction d'erreur 
{MCE}. C'est le théorème de représentation de GRANGER (1983). 


Les modèles à correction d'erreur représentent une catégorie de relations 
particulièrement intéressantes quand on travaille sur des variables intégrées 
d'ardre 1, Nous reviendrons sur cette question dans les chapitres V et VL 


Les résultats qui viennent d'être exposés permettent de préciser Jes conditions de 
validité de l'économétrie. 


Les conditions de validité 

L'économétrie traditionnelle sur des variables non stationnaires reste légitime à 

deux conditions principales : 

— Ne pas mélanger dans une équation des variables ayant des ordres 

d'intégration différents. 
On rejoint une question souvent abordée de façon intuitive : celle de 
l'homogénéité d'une équation quand dans une même équation il y à des 
variables qui sont affectées d'une tendance et des variables stationnaires, il est 
en général impossible d'obtenir une estimation satisfaisante. 

C'est un problème connu depuis longtemps pour des fonctions comme : 

1, =@}, + ar; +a, pour lesquelles il est généralement impossible d'obtenir des 
cstimateurs significatifs de az. 

— . Dans une équation de nature dynamique quand on introduit des décalâges sur 
les variables 1l faut respecter le méme numbre de décalages sur toutes les 
variables où au moins sur les principales variables. 

Par conséquent, la théorie de La coïntégration valide plutôt la démarche 

£conométrique sous des conditions pas trop restrictives, 

Dans un travail d'économétrie, il est donc recommandé de commencer par 

Chercher l'ordre d'intégration des variables (test de racine unitaire puis différences, 

puis test de racine unitaire etc...) 

Ensuite, il est bon de s'assurer que les relations éconumétriques qui ont été 

obtenues appartiennent bien à l'espace de cuïntégration des variables (non unicité 
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en général de la relation) en déterminant l'ordre d'intégration de là série des 
résidus el en montrant que l'ordre d'intégration des résidus est inférieur à l'ordre 
d'intégration des variables. 

Le plus souvent on obtient différentes relations possibles pour estimer et on 
retiendra de préférence celles qui ont aussi le caractère de relation de 
coïntégration c'est-à-dire celles dans lesquelles Les variables sont coïntégrées. 
Moyennant ces précautions on peut considérer que la plupart des propriétés 
statistiques restent vajables. 


Chapitre V : Application à la consommation 


La consommation est un domaine qui se prête particulièrement bien aux 
applications de l'économétrie. Les variables qui s'y rattachent, qu'il s'agisse de 
la consommation elle-même ou du revenu qui constitue la principale variable 
explicative, font preuve d'une assez grande stabilité et d'une forte incrlie, Elles 
peuvent être traduñes par des spéctficauons relativement simples. 


La consommation peut être analysée à 1rois niveaux qui correspondent aux 
trois sections de ce chapitre. 


SECTION I ; LA CONSOMMATION GLOBALE 


Le modèle de base. déduit de ta thévrie générale possède une propension 
moyenne à consommer décroissante ct une propension marginale stable, mais 
inférieure à l'unité. La traduction la plus simpie ost donnée par : 

(1) &=av, 42 +u, 

avoc U<a <1 et a: >Ù 

où €. représente la consommation par tête à prix constant et y, le revenu par 
tête à prix constant. 

Les très nombreuses études empiriques faites à partir de ce modèle ont bien 
confirmé la stabilité de la propension marginale à court terme, mais ont 
infirmé la décroissance de la propension moyenne. En particulier, une étude de 
Kuzneïs sur la consommation aux États-Unis à étahli La quasi-stabilité dc la 
propension moyenne à consommer depuis 1865 ‘!}. 

On sait que parmi les très nombreuses relermulations, la plus inéressanté est 
celle proposée par M. Friedman autour du concept de revenu permanent {2. 


US. Kusneis “National Product since 1869" NBHR 1946. 


{21 On trouvera un panorama de ces hypothèses et tentatives d'explication dans l'ouvrage 
de 1€, Eicher "Consonunation et épargne" Ed, Sirey 1961. 
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Toutefois cette approche pose le problème de la mesure du revenu permanent &t 
de la consommation permanente. M. Friedman fait l'hypothèse que les 
consommations transitoires sont nulles en moyenne pour une population 
donnée et que les révenus passés interviennent dans la composition du revenu 
permanent avec des poids géométriquement décroissants, {l en résulte une 
spécification à retards échelunnés : 


(2) e =kU- Ar + ap An, | avec 0<A<1 


L'application de la transformation de Koyck permet de passer à une 
formulation auto régressive, ef supposant que np tend vers l'infini. 


Ac KG A5 À" 5] d'où 

(ce =K(l-A)y, +de, 

dans laqueile la propension à consommer ssl : 

à court tete : K{]- À) 

à long torme : K 

On peut remarquer que dans cette spécification, l'effet d'une variaion 
transitoire de revenu est parfaitement déterminé, ce qui est en opposition avec 
la théorie de départ. En effet, celle-ci postule qu'une variation de revenu 
transitoire affecte la consommation dans une proportion variable, Ce n'est à 
qu'un exemple du caractère réducteur de la formalisation par laquelle un 
ensemble de propriétés ou d'hypathèses est traduit par une rclation 
mathématique. 

La théorie du revenu permanent conduit à estimer une équation du type : 

(4) ec, =ay +acs_1 + as +4, 

avec 4 =K(l-4) 

d; = À 

Cette spécification, qui s'ajuste généralement bien sur données 
macroéconormiques brutes, présente de nombreuses faiblesses. 
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a - L'évolution du revenu et de la consommation suivant des trajectoires à peu 
ptès identiques et les comporiements faisant preuve d'inerue, les variables 
explicatives ti et y. seront très souvent Colinéaires. 

b - Si la fonction de départ, équation (2), comporte des termes d'erreurs 
indépendanis au cours du temps, l'application dé la transformation de Kayck 
rend probable la présence d'autocorrétarion dans l'équation (4}. De plus ie 
modèle comporte une variable cndogène décalée, le test de Durhin Watson n'est 
pas fiable ut les coefficients de régression sont biuisés. 

&- La formule de l'éguation (4) ne contient que les déterminants principaux. I 
convient de la compléter pour lui permettre de prendre en compte l'influence 
éventuelle de la répattition des revenus {!!, du chômage qui peut exercer un effet 
de frein (21, des liquidités ou du taux d'intérêér dont l'effet s'exerce 
principalement sur la demande de biens durables. 

d - Les variables y, et c; sont généralement intégrées d'ordre 1, parfois 
d'ordre 2. Il faut s'assurer que les autres variables ont un ordre d'intégration 
équivalent et vérifier que Ja relation obtenue appartient hien à l'espace de 
coïntégration des variables. Les travaux anciens ignorent cotte précaution. 
Les modèles précédents peuvent être ufilisés pour unc première approche, mais 
restent peu performants. 

Exemple : Estimation d'ünc fonciion dé consommation sur le Niger en 
données annuelles 1962-80, Avec le modèle le plus simple : 


c, = Ô,88 PIRRU — 6,98 
146,7} 1,64) 


R?=0,997 DW=1,72 p=0,475 


1} Soit en distinguant les revenus du travail et les revenus du capital dans une approche à 
la Kaldor qui peut être pertinente pour des économies caractérisées par de fortes 
inégalités de revenus, soit en séparant les prestations sociales qui sant presque entièrement 
Cunsouuriées, dés autres révenus dans le cas des économies d'Eurupe occidentale. 


(2) Ce qui montre que la traduction du revenu permanent utilisét est bien imnparfaite, car 
elle n'inlègre pas oct élément d'anticipation commune ele le devrait, 
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PIBRC = PIB moins les revenus de l'uranium. Ceux-ci sont en majorité versés à 
l'État ou transférés à l'étranger et ne servent pas de support à la consommation 
des ménages. L'utilisation directe du PIB ou du revenu réel en variable 
explicative donne des valeurs de la prapension marginale à consommer trop 
faibles, généralement inférieures à 0,7. 
Le modèle tiré de la spécification la plus simple de la théorie du revenu 
permanent donne : 
c, = 0,853 PIBRL'+0,113c,_ — 9,621 

196.43 (10,45) [19,583 
R?=0.998 DW=1,69 
Estimation obtenue en normant l'équation par PIB RU. 
L'analyse graphique des résidus fait apparaître un écart important pour les 
années 1973 el 1974 correspondant à 2 années de grande sécheresse au cours 
desquelles la consommation à moins baissé que le revenu. L'introduction d'une 
variable muetté CL qui vaut ! pour 1973 et 1974 et à ailleurs permet d'effacer 
ces écarts : 


, =0, CA CERS ED Oil 1 +6, 993CL-11,035 
{76.4} 12,80} TÉ) 


RÈ=0,9988 DW=1,57 


La propension à consommer à lang terme est alors de 0,97, ce qui correspond à 
un taux d'épargne moyen de l'ordre de 3 %, compatible avec les chiffres 
Crégistrés pour cette variable pendant Ia période. 


Remarque : dans cet exemple, on peut aussi reprendre le problème à l'envers, 
en supposant que l'on connaît mieux la consommation que le PIR ou le revenu. 
Le comptable national utilisera les estimations des chiffres de la 
consommation pour faire son évaluation du revenu ou du PIB en se donnant à 
priari les coefficients de la fonction. on touche 1à une Hmite évidente de la 
pracédure d'estimation dans çe cas particulier. 

Cette classe de fonctions s'ajuste bien dans la plupart dés cas, en données 
annuelles ou trimestrielles. Toutefois. leur qualité est sans doute surestimée du 
fait de la stabilité des évolutions. On constate, en effet, qu'elles sont 
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incapables de décrire les reltournements ou Ics inflexions importantes. Ce 
sont, pur conséquent tout ou plus, des spécifications exploratoires. 

La nécessité de mieux suivre les variations de la consommation a conduit à 
diverses tentatives de refarmulalion : 

- En ceconsiruisant la variable revenu permanent dircetement sur Ja base de 
revenus passés et anticipés. C'est par exemple ce qui est l'ait dans le modèle 
ICARE de l'IPECODE. Le revenu total est scindé en 5 composantes done la 
fraction permanente est cstimée à partir des valcurs des 4 trimestres précédents 
ei des trimestres courant {1} 


- En spécifiant unc fonction dont la variable expliquée est la propension 
permanente à consominer. C'est ce qui 4 été fait dans le modèle PITI de la 
direction de la prévision {2} On peut dire aujourd'hui que cetie solution a 
l'avantage d'utiliser des variables stationnaires. 

Mais la tentative de reformulation à la fois la plus ambitieuse et la plus 
intéressante est celle qui a été menée principalement en Angleterre à partir de 
la fin des années 70, par Hendry et Davidson. En travaillant sur la 
consotnination, ils ont construit un nouveau type de spécification. déduit d'une 
catégorie d'équations désignées dans la littérature sous le nom de modèle à 
correction d'erreur (MCE) qui présentent des propriétés économétriques 
particulièrement intéressantes. L'analyse de la consommation 4 été l'occasion 
d'un intérêt renouvelé pour le type de modèle. 


La démarche de Hendry et Davidson se présente comme une généralisation de la 
démarche de Friedman avec la (héorie du revenu permanent. Conime ce dermier, 
ils Se proposent de rendre compte simultanément de la stabilité de la 
propension à consommer à lung terme et de son instabilité dans Le court terme. 


(CARE mouète conjoncturel de l'éconcnie française" Revue de l'IPECODE n°1. 
Mars 1983 et "ICARE L98G le modèle ICART eu noavelle base de comptabilité nationale” 
Revue de l'IPECODE n° 22. Déc, 1948, 


(2Y <pITI : éléments pour un modèle à prix implicites trmeslriel intéuré”, Statistiques et 
études financières n°41, 1980). 
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Si la propension à long terme csi stable, lé taux de croissance de ja 
consommation est cn longue période égal à celui du revenu : 


LS ee 

CE EU 

Cette condition est par exempte remplie par 
C, = K(glr, soit JogC, = log* +logK(e) 
en posant : c, = log C, et y, = log, 


(5) 


Î vient 

(6) Cor +k 

En passant aux différences 

[T} G-cs =}; 8 puisque ci = log = togcleÊEL } 
1-1 1=1 


AC, 


dont l’approximation est =g 


3-1 
Le modèle {6} se généralise par l'introduction d'un polynôme de retard sur €, et 
d'un autre sur y, : 
(8) ali =B(L)y,+k 
Dans la pratique on retictit des polynômes du premier dégré dont une forme 
particulière est : 
D ce =k+a.. +foy + Bye 
Cette relation étant vérifiée aussi ent -L. il vient : 
(0) ee = 2)+ Poe 3 + Afin — Ve) 
Dans le cadre d'un régime de croissance équilibrée, J4 taux de croissance de la 
consommation est constant ct il est égal au taux de croissance du revenu (l). 


( ün parvient au (oêrme résultat co posant que l'équilibre Statique de long cerme impose 
que les taux de croissance de © et de Ÿ soicnt ouls et que Jes valeurs de Ÿ el C soient 
constantes, à partir dé : 
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d'où : 

Er TE A TEA Ere2 Ve 1 Je F2 TE 
soit g= 8 + Ho + P18 

AD 1=@ +f0+fA 

Ce que l'on pout encore écrire sous la forme : 

O2) 1-m = Bo+A TT 


et donc 
(3) æ=1-7Y 
(4) B=Y-f 


Les valeurs obtenues dans {13} et (14) peuvent être reportées dans (9) 

US) 6 =&+ Pet + Boyer +CY— Bois 4 

(6h cc. = + Batre —312+ FO © 1) 

(6) cc =k+ for - y) Yes — vs) 

Cette équation constitue la forme la plus simple d'un modèle à correction 
d'erreur, Le taux de croissance de la consommation dépend du taux de croissance 
du revenu et d'un terme de rappel qui est mesuré par l'écart entre log c, et log y, à 
la période précédente. 


=1 


Ce terme peut chcorc s'intérpréter comme log et se compare dunc à la 


1-1 


€, 
valeur de long 1erme du rapport n . Si la consommation a augmenté moins 
J 


vitg que le revenu à [a période précédente, le terme de rappel accélère la 
censommation à Ja période présente. Dans le cas contraire, 11 freine la 
consommation. 

Un modèle à correclion d'erreur peut être complété par diverses variables 
explicatives. 


(9 cr = =K+iæ] _- bee] + fo Y— _1)+ CB + fi 1x _I cf Cuthbertsori, Hal) 
et Taylor “Apptied Économétnic Techniques” Philip Allan 1992 p 103. 
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Un exemple de fonction de consommation de type MCE peut être trouvé dans le 
modèle DMS 4 {1j dans lequel ont € ajoutés Le taux d'inflation et Le taux de 
chômage anticipé. Là formulation est un peu différente de l'équation (16) dans 
la mesure où elle est résolue par rapport à l'épargne ei non par rapport à la 
consomiMmalivn. 

On passe nalurellement de l'un£ à l'autre par le fait que E, = Y,-C4 Nous 
reviendrons sur les raisons qui peuvent conduire à estimer plutôt une foaction 
d'épargne dans les modèles multisectorieis, à [a fin de ce chapitre. 

Le modèle à correction d'erreur présente une propriété remarquable qui à été 
démontrée par Granger en 1983 ‘2, Un ensemble de variables co intégrées 
d'ordre {1,1} peut être mis sous la forme d'un modèle à correction d'erreur don! 
toutes les variables sont stationnaires ct dont les coefficicnls peuvent être 
estitués par les méthodes de l'éconvmétrie classique sans risque de corrélations 
fortuites. Le résultat, connu sous fs nom de théorème de représentation de 
Granger, valide de façon générale La démarche du MCÉ pour une classe 
unportante de variables. 

Les modèles à carrection d'erreur peuvent être urilisés paur de nombreuses 
varfables économiques. Il est passible de s'affranchir de l'égalité du taux de 
croissance de la variable explicative et de Ia variable expliquée &t il est 
également possibie d'appliquer le modèle à plusieurs variables explicatives 
figurant dans une même équation. 


Supposons que la variable X,, dépendant de 24, ait Un taux de croissance à long 
terme égal à # fois celui de Z, : 
DUR € Ki. | = Os 


L'équation du MCE est : 


Q} ps 4 modèle dynamique multisectoriel - Collectins de lTNSEE série C 0° 129 - [987 
Chap. IV - p 146 - 150. 


{2} voir aussi Engte et Granger “Co intégration und Error Correction Représenlation, 
Estimation and Testing” Econormétrica 1987. 
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(7) x 24 + ox, + Po + Âizs_) 
La condition de long terme {11 s'écrit désormais 
(8) = 00; + fo + f 
Soit 
_ Both 
Ÿ 
d'où a =i-7 « B=n- fo 
en remplaçant dans (17}: 
(9) x = 4#+il— 7x, + Bo + (YO — Poe 
(20) x Et Por 21) Ye à sr) 
Cette équation peut être estimée par les MC sans restrictions sous la forme : 


(8) I-@ =} 


(20 bis) At, = ay ré Az, +aox_] +é425 +4, 
Les valeurs des coefficients de (20) sont données par : 


- 63 


42 


k=ag fo=a 7=-4 À 
L'équation (20 his) conduit à régressur sur des variables dont l'ordre 
d'intégration est différent : Az, est liojet y, conmme z, , som Zi) 


Engle et Granger proposent une dématche en 2? étapes permettant d'estimer 
directement l'équation {20} sans avoir à passer par {20 bis). Dans une première 
étape, la relation de coïntégration, donc de long terme, entre x, et z, esl 
estimée. 


Li 
Elle fournit la valeur du coefficient Ÿ, sait ©. et étabtit en même lemps que 


h Di 
x, 02, comme t,_1 492,4 Sont stanonnaires par définition de la 
cointégration, 


D 
L'équation (20) peut alors être estimée en utihisant la valeur dé Ô : 


Ax,; = k + Bo — YÉxs 1 — ÿ 41) +u, permettant d'obtenir dos valeurs de k, Bo 
et y. 
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Si maintenant Y, dépend de Z, et de W, avec des réponses à long terme, telles 
que le taux de croissance de X soit égal à 4 fois celui de Zet & fois celui de W, 


le modèle s'écrit : 
QI) x =k rot + Por + Him + bo +811 


et l'équation {18} devient : 


Q2 1-a = 0h _éoté, 
ê E 

d'où: 

&| =1-7 

8 = %Y— Po 

Gi =E7-60 


Soit en remplaçant dans (21) : 
(23) x, =k+Q— pas + Bon (07 fo)z-1 + or +(E7— Én)w, 
(24) x, xs = k+ fol 2-1) + 600m, wi 1) 7x1 — Da, — 061) 
Équation qui peut être estimée par Les MCO sans contrainte sous la farme 
(25) x, = Go +@ 4, +43 AW, + ax, ass. 1 +4sW_1 +U 
Les coefficicnts de 42%) sont : 
k=ay Ho=a Go=e Y=-0 De F2 

a; 43 
Les réponses à long terme de X, les élasticités à long terme par rapport à Z et 


2 a a a ; 
W sont donc données par + à -+ respectivement, tandis que Îes 


&2 43 
élasticités de court terme somt respectivèment a] et 49. 
De tels modèles sont couramment utilisés pour estimer des functions 
d'investissement, d'importations ou d'exportarions (let sont intégrés dans 


{1 parmi de nombreux travaux cf : J.P, Urbain “modèles à correction d'erreur et fonctions 
de detoande d'inportations agrégées", Économie et Prévision 1950 : $. Capet ct Pb. Gudin 
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certains modèles macroéconomiques comme MEÉTRIC. La version 
METRICX 1) réestimée en 1988 fait une large place aux fonctions de lypé 
MCE. Certains auteurs utilisent aussi l'expression “modèle à correction 
d'équilibre” qui peut sembler plus approprié que "correction d'erreur” ct met 
bien cn évidence l'idée d'évolution de long terme ou de situation d'équilibre de 
longue période par rapport auxquellts les écarts déclencheraient des réactions 
de correction {*}, 


SECTION 2 : L'APPROCHE SEMI-GLORALE 


4 l'intérieur de ka variable consommation, il est parfois utile de distinguer les 
bien durables et les biens non durables. 


Les bien nôn durables tels que l'alimentation ou les services peuvent être 
correctement représentés par un modèle du type de ceux qui ont été présentés 
dans 18 section précédente. [ls ne posent pas de problèmes particuliers. 


Les biens durables, en raison de leur durée de vie précisément et accessuñrement 
de teur valeur unitaire, présentent des caractéristiques particulières. 


Comme le notait M.K. Evans dès 1969 ‘#), il est possible de déplacer dans le 
temps la décision d'achar en la reportant ou au contrairc en l'anlicipant, ce qui 
confère à la demande une plus grande variabilité et d'autre part le stock existant 
de ces mêmes biens exerce un effet de frein sur l'atervissement du parc. De plus, 
compte tenu de leur mode d'acquisition, qui fait sauvent appel à la constitution 
d'une épargné préalable ou au crédit, la demande pour cé type de biens est 
beaucoup plus sensible aux variables monétaires et financières que la demande 
dé biens de consommation en général. En somme les biens durables sc situent 


de Vallerin “Fonctions d'importations el d'exporlations : f'apport de la théorie 
écunométrique récent” Économie ct Prévision 1993 ; J. Glachant et J.F. Nivet “Deux 
études macroéconomiques de l'investissement”. Écunonie cl Prévision HJR9. 


Ce modèle METRICX et l'étude de la variante dévaluation”, Économie et Prévision, 
538. 


{4 C'est l'expression vtilisée par Hendry (Equilibuiutn Correction Model) - cf D. F. Hendry 
“Dynamic Econometrics"” Oxford University Press. 1995. 


($} MK. Evans "Macro économie Activity", Harper ct Row, 1969. 
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dans une position intermédiaire entre les investissements et la consommation, 
comme le traduit l'expression “biens d'équipement des ménages” wrilisée par 
J'iINSEE. 


La modélisation de la consommation de biens durables présente certaines 
similitudes avec celle de l'investissement, mais S'en distingue par certains 
aspects tels que la structure des retards de décision et les délais de réalisation. 


Les modèles de consommation de biens durables sont à la fois nombreux et 
divers. C'est depuis longtemps un des domaines principaux d'application dé 
l'économétrie. Faute de pouvoir rendre compte de cette diversité, on se limitera 
à une présentation des principes utilisés pour la construction de çes fonctions, 


Comme pour la FBCF, on peut distinguer un stock désiré et une fonction de 
réalisation qui décrit l'ajustement à court ttrme du stock effectif au stock désiré. 
Le stock désiré, qui peut aussi être un stock d'équilibre dépend uniquement du 
revenu des ménages et du prix relatif de ces biens par rapport à l'ensemble de Ja 
consommation où par rapport aux setvices qu'ils remplacent. 

(D S*= f{y.Pr) 

Cette fonction peut comporter un mécanisme de retard décrivant l'ajustement du 
stock désiré aux variations du revenu. L'utilisation du revenu permanent ou du 
reévénu anticipé à moyen terne permet d'éviter ce détour dans certains cas. Un 
modèle à correction d'équilibre peut également être utilisé pour déterminer le 
stock désiré. La fonction de réalisation décrit l'ajustement du stock effectif au 
stock désiré à partir d'uné distribution de retards sur les variables qui 
déterminent là décision d'achat des ménages, c'est-à-dire par exemple la 
variation de leur revenu, les conditions de crédit ou encore l'évolution 
anticipée des prix des biens durables. 


$, S* 
€2} Sy —S ou SE endse- 7-1 64 Et où X) g1(b(L)Ay,,i,pf) 
1-1] 1-1 


| “ : 
La formulation en —— présente l'avantage de faire apparaître le taux de 
1-1 
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croissance du parc, lorsqu'on passe en logarithme. La fonction g, peut être un 


+ 7 
simple mécanisme d'ajustement partiel du type À(S —5_1) ou Es . 


DE 
La fonction g2est additive ou multiplicative selon la spécification retenue. 


La combinaison des équations (1} et (2; permet d'exprimer la fonction 
décrivant l'accraissement du stock de biens durables ou son taux 
d'accraissement. 


La demande de biens durables est égale à l'accroissement du parc auquel s'ajoute 
ta demande de renouvellement. Cette dernière peut être obtenue en appliquant 
un taux de dépréciation # au stock existant : 


D,, = 05,1 où # => T étant la durée de vie moyenne des biens. Mais cette 


solution ne rend pas compte du rythme selon lequel Île stock à été cunstitué, 
Ainsi un parc automobile peut être très récent et la demande de renouvellement 
ainsi calculée surestimera grossièrement ja demande de renouvellement 
effective. De même cette méthode ne peut rendre compte des éventuels à-coups 
dans la constitution du stock, du iype crises pétrolières pour un parc 
automobile. 


il est donc préférable, quand l'information est disponible, d'appliquer une loi 
de déclassement aux différentes classes d'âge constituant le stock, Pour des 
biens tels que ies automobiles qui font l'objet d'un enregistrement et d'un suivi 
administratif, cela ne pose guère de probièmes. Pour Îles autres catégorics de 
biens durables c'est généralement beaucoup plus compliqué puisqu'il faut 
commencer par construire le stock existant à partir des achats réalisés dans lé 
passé. Rappelons que pour la comptabilité nationale, les biens durables sont 
réputés consommés dès leur entrée dans le patrimoine des ménages. ce qui 
justifie qu'il ne suit pas tenu de compte des stocks détenus par les ménages. 


Parce qu'elles correspondent à des comportements différents, il est donc 
intéressant de pouvoir distinguer Ja consommation de biens durables et celles 
de biens non durables dans un modèle non sectorialisé, On peut ainsi améliorer 
ia précision avec laquelle l'ensemble de la consommation est décrit. 
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Dans les grands modèles sccrorialisés, la consommation est décrite par 
produits, comme on va le voir dans la section suivante. 


SECTION 3 : L'APPROCHE PAR CATÉGORIES DE PRODUITS 


Dans les études détaillées comme dans les grands modèles matroéconomiques, 
une projection globale de la consommation est insutfisante. D'une part il se 
produit des déformations continues de la structure de la consommation sous 
l'influence des élasticités revenus, déformations dont il faut pouvoir rendre 
compte. D'autre part, l'équilihse offre - demande doit être assuré par catégories 
de produits et non plus au niveau de l'ensemble des catépories. Dès lors que des 
catégories ont été identifiées, i] n'est plus évident que les excédents d'offre 
d'une catégorie puissent combler les déficits d'une autre. Or les conséquences 
des situations d'excédent d'offre ou de pénuries ne sont pas symétriques. Une 
pénurie est comblée par recours à des importations en économie ouverte, 
tandis qu'un excédent conduit à une augmentation non désirée des stocks de 
produits finis, puis à un ajustement du niveau d'activité. Qn voit que l'équilibre 
macroéconomique, qui pourrait être réalisé au niveau de l'ensemble des 
catégories de produits, peut être remis en causé par des désajustements 
sectoriels. 


L'analyse de la consomation par produits se réalise généralement en deux 
étapes. 


On pourrait imaginer de déterminer simultanément toutes les fonctiuns de 
demande à partir des conditions dé maximisation d'une fonction d'utilité, sous 
contrainte de revenu. La consommation globale serait alors obtenue par 
sommation des consutnmalions par produits Cette dérnarche est impossible 
pour plusieurs raisons. Tout d'abord les quantités des différents produits, qui 
ont été consommées, ne correspondent pas à des quantités demandées, mais à 
des quantités d'équilibre reflétant le jeu de l'offre et de la demande sur les 
marchés. li serait donc incorrect de vouloir en rendre compte uniquement à 
partir des seules fonctions de demande qui sont celles qui se déduisent de Ja 
maximisation sous contrainte de la fonction d'utilité. I faut rappeler que les 
fonctions de consommation par produit ne sont pas des fonctions de demande 
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(on peut se reporter à ce sujet au chapitre II seclion 2 du présent ouvrage). Une 
deuxième raison est qu'une telle procédure aurait pour effet d'additionner les 
erreurs faites sur les consommations par produils au niveau de la 
consommation globale puisque : 


€; = EC; 
1 


donc si 
Cie = hi Ve Pins Hu 
Ce = Eh Pir...)+ Ets 


Or il n'y a guère de raison à priori pour que les erreurs faites sur les 
consommations par produits se compensent. 

De ce fait la consommation globale serait déterminéc avec une assez grande 
marge d'erreur. Celte marge se retrouverait cneurc amplifiée au niveau du calcul 
de l'épargne que l'on obtient par diflérence avec Le revenu. 

En effet, compte tenu des proportiuns respectives de l'épargne et de la 
consommation dans lé revenu, une erreur laible en pourcentage sur la 
consommation se traduit par une erreur forte en pourcentage sur l'épargne. 


Supposons que la consummation représente 90 % du revenu et l'épargne 10 % 
Une erreur de 1 % sur la consommation entraîne une erreur de 9 % sur l'épargne. 


C'est pourquoi on commence par déterminer l'épargne à l'aide d'une fonction de 
comportement Même si cétte foncliun est de qualité inférieure en terme 
d'erreur quadrätique moyenne rapportée à la variable expliquée 0), elle reste de 
meilleure qualité qu'un résultat déterminé par solde de la consommation sur le 
revenu. 


| 
(L) E 1 pr TX, —$ e 
ÿr 


192 Économéirie appliquée 


Ainsi, par exemple, a fancrion d'épargne ajustée” construite dans lé modèle 
DMS d à une précision de 3,2 % {1} Cette fonction est déduite d'un MCE pour la 
function dt consommation. Pour obtenir le même degré de précision, il 
faudrait une fonction de consommation ayant une erreur quadratique moyenne 
inférieure à 0,5 %. 


Une fois estimée l'épargne, on obtient la consommation globale par 
différence. I laut ensuite répartir celle-ci entre lés différentes catégories de 
produits. La répartition se fait par un système d'équations inspiré de la 
microéconomie. Le plus connu de ces systèmes est celui de Houthakker et 
Taylor. [l conduit à des spécifications dans lesquelles ia consommation d'ur 
produit à la période t dépend de la consommation globale, du prix relatif du 
produit et de La consommation passée du produit : 

(AY 6j = ag +4, +a2pry +axcg-5 + 

La consommation décalée peut imervenir avec le signe positif, auquel cas elle 
traduit un effet d'habimde, ou avec le signe négatif traduisant un effet de stock. 
Un coefficient à3 positif sera associé aux fonctions de consommation de biens 


non durables et un coefficient négalif aux fonctions de consomation de biens 
durables. 


Ce type de fonctiun est compatible avec une spécification litéaire comme dans 
l'équation (1) ou log linéaire comme (2) ou serni log comme {3) et (4). 

(2) Jogci = au +a loge, +aslogp,, +aalogeis_i +4, 

(3) €j = ag ta loge, +a;logp,, +aslogc,_1 +4, 

(4) Jogcs = a6 +0, +pp,, +dati_t + 

Une catégurie de produit dont le coefficient budgétaire est à peu près stable 


pourra être représentée par un modèle log linéaire, Par contre la consommation 
d'une catégone de produits dont la part régresse dans La consommation globale 


(L} DMS 4 modèle dynsunique multisccturiel, Collections de l'INSEE séric C n° 139 1987, p 
148 et suivantés. 
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comme la consommation alimentaire par exemple à des chances d'être mieux 
décrite par uné fonction de type (3). 

Exemple ; analyse de la consommation de carburant en France de 1983 T à 
1993 11. On commence par vérifier que la consummation de carburant (ec,} la 
consommation finale {c,) et le prix relatif des carburants (Pre;) exprimés en 
logarithmes, sont coïntégrés. 

Test de Dickey - Fuller : 


Variable Lec, Ec, EPre, | Acc AEC, APrt 
Valeurs du test 0,919 14,7 -(,GD08 | -7,246 -2,059  -5,526 
Valeur critique à 5 % 1,949 -1,949 -1,949 | -1949  -(,949 -1,949 
Valeur critique à | % -2,619 -2,619 -2,619 | 2,621 2,621  -2,621 


D'où il ressort que les variables sont routes imégrées d'ordre 1. 

Test de Engle et Granger sur Lec,, Le, et LPree 

Valeur du test : - 5,80 

Valeur critique à 5 % : - 3,954 

Valeur critique à 1 % : - 4,662 

résultat qui établit que les variables sont bien cointégrécs. 

Après diverses tentatives on constate que la consommation de carburant 
décalée n'est jamais significative. D'autre part, la comparaison graphique entre 
valeurs calculées et valeurs observées montte un léger déplacement après le 
2 trimestre 1987. L'introduction d'une variable muette (D} prenant la valeur Q 


jusqu'au 2 trimestre 1987 et la valeur 1 après se révèle significative. Un tesi de 
restriction a priori confirme la pertinence de cette variable. 


{5) Lcc, = ag t+a;lc, ta;LPrc, +u, 
(6) Lec, =ap+alc,+a;LPre, +a,D+4, 


Zwr de l'équation (51=0.0171 


Zér? de l'équation (6) = 0,0143 


194 Économétrie appliquée 


Ehr-Ele  0,0171-0,0143 


8 1 : 
— = ——  —— = 6,65 
A 0,014 L 
EG? 0,0143 


TL 34 
La valeur critique de la table de Fisher est : 
Frt34)= 4,17 par cacès 
L'ajustement définitif est : 


€7) Pro Is 2e +0,002D-0,017£Pre, Fe R? = 0,968 
(283 (3.501 


estimation obtenue sous AR{4} pour tenir éompie d'un fort mouvement 
saisonn)er. 

Incidemment, ce résultat montre que contrairement À ec qui est souvent avancé, 
Félasticité de la consommation de carburant est supérieure à l'unité et que la 
consommation à même eu tendunve à augtnenter récemment. {ne 
décomposition plus fine montre que celle cruissance est essentiellement due au 
carburant diesel. 


Dans un excreice de modélisation de la consommation par produits, rien 
n'assure que la somme des consommations par produits soit égale à La 
consommauon globale calculée séparément. 


Dans le cas général : 
Zcysy,-E, 
1 
I y a deux solutions qui permettent d'assurer la compatibilité entre ces deux 
variables. 
La première consiste à laisser libre une catégorie de produits et À la déterminer 
pur sulde. 
Cr=h-k-EcG 
de Tes 4) it 
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Compte tenu des phénomènes d'accurmulation d'erreurs qui ont été évoqués 
précédemment, il est souhaitable que Ha catégorie détérminée par solde 
représente une fraction assez importante de la consommation totale. 
C'est pourquoi dans certains exercices de modélisation de la consommation par 
produits, on soîde sur le poste services qui représente près du ticrs de la 
consommation totale. 
L'autre solution consiste à répartir l'écart constaté cntre les catégories de 
produits en proportion de leurs valeurs. 
si CE 

Ec; 

î 


Chaque poste est réévalué par le coefficient m. 


Chapitre VI : Application aux échanges 
extérieurs 


Les échanges cxtérieurs constituent un autre domaine d'application intéressant 
des instruments de l'économétrie. On sait que les relations économiques avec 
l'extérieur sont constituées de différentes catégories d'opérations 
mouvements de biens ct services, mouvements de capitaux, opérations de 
répartition etc. On s'intéressera exclusivement aux opérations commerciales 
qui sont les plus simples à modéliser. Certaines opérations peuvent leur Être 
rattachées comme par exemple les crédit à l'exportation ou les investissements 
directs qui ont tendance à dépendre des flux de marchandises. 


La modélisation des échanges commerciaux présente l'avantage de pouvoir se 
faire avec des spécifications relativement simples, elle aussi, et élle permet 
d'aborder quelques problèmes intéressants. 


On commencera par les importations, puis dans une deuxième section, on 
étudiera le cas des exportations. 


SECTION 1 : LES IMPORTATIONS 


Elles dépendent à titre principal de la demande intérieure qui peut être 
représentée par l'activité économique ou le PIB ct de la compétitivité des 
importations par rapport aux productions nationales pour les produils qui sont 
substituahles à la production intérieure, Pour de nombreux pays engagés dans 
une démarche d'intégration régionale ou d'ouverture extérieure, il existe une 
poussée autonome des importations qui traduit cette intégration ct l'évolution 
des spécialisations qui en résulte on) 


(D Sur l'évolution des spécialisations internationales cf. H, F. Henner “La spécialisation 
imemationale de l'économie française" Revue Économique 1976 et K. Abdel Rahman 
“Réexaren de la définition de la mesure des échanges croisés de produits similaires entré 
les nations‘, Revue Économique, 1986. 
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Ces influences peuvent être traduites dans Ia formule : 


Im = f(Y, Pr, À 
Connue sous le nom de “Formule de Luxembourg” (1. 


Dans la mesure où des tensions sur les capacités de production peuvent 
entsaîner des retards de livraison de là part dés entréprises intérieures ôu 
conduire à une certaine saturation des capacités de réponse de l'offre intérieure, 
on peut compléter cette formule par ut indicateur dés taux d'utilisation de 
capacités. Lorsque les capacités sont utilisées plus intensément, Îles 
importations augmentent. 


La réponse des importations aux variations d'activité ou de compétitivité est 
généralement rapide. En données annuetles, elle se produit entièrement au 
cours de la période. En données trimestrielles, les retards ne dépassent guère 
2 trimestres et il n'est pas nécessaire de construire de structure de retard dans la 
plupart des cas. 


La fonction d'importation est compatible avec une spécification linéaire, ou 
avec une spécification log linéaire à élasticité constante : 


Imi=4l +4 Pr, +aaFu,. +us +4, 

à, représente donc le contenu en importations d'un accruissement de revenu ou 
d'activité. 

On peut aussi distinguer les différents postes de la demande finale. 

Ce qui conduit à estimer des contenus en importations différents pour ces 
catégories : 

Imt=al,+aG, +a3X, +aal, +as Pr +asTu,. +47 +u, 

Toutefois, il est quasiment impossihlt dans cette forme d'éviter les problèmes 
de colinéarité généralisée. 

Avec des élasticités constantes le modèle s'écrit : 


«a Cf B. Bobe et F.H. Derycke ‘Projections des échanges extérieurs et balances des 
paiements", Econotmca 1975. 
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} Jogimf = 0; log ÿ, +0 log Pr, + 03 log Tr. + 04 +4, 
œ, est l'élasticité revenu des importations 
a, est l'élasticité prix des impartations. 


De nombreuses études éconcmétriques ont été réalisées dans le passé à l'aide de 
la relation (1} ou de formes voisines. 


Exemple : importations italiennes on données trimestrielles du chapitre ! 
section 3. 


Avec une spécification à élasticité constante, on obtient pour les importations 
italiennes de produits manufacturés sur la période 1978 I - 1492 IT: 


logImr = 1,916 log Y, — 0,147 log Pi, + 0,271 D90— 11,89 
{22.34) ‘4,913 42,81 (10, 14; 


R2=0,976 DW=2,02 p=0417 
Im,= importations italicnnes de produits manufacturés, déflatécs 
Y, = PIB déflaté 

= prix des importations de produits manufacturés 


D98 est une variable muette qui vaut L pour je premier trimestre de 1999 et 
permet d'éliminer un écart important invxpliqué pour cette période. 


On vérifie que og Im, et log Ÿ, sont cointégrés : 
Vecteur cuïntégrant : 

Log Im: 1 

Log Y::-1,573 

Test de Engle et Granger : - 5,016 

Valeur critique à 5 % : - 3,44 


Dans la plupart des pays, les élasticités revenu sont bien estimées et varient 
entre un peu plus de 1 et 1,8 selon les pays et les périodes. Le Japon est un des 
râres pays de l'OCDE à avoir eu pendant longtemps une élasticité inférieure à 1. 
L'estimation de l'élasticité prix n'est pus aussi satisfaisante. Il est souvent 
difficile d'obtenir des coefficients significatifs. ét du bon signe. La raison de 


200 Économéirie appliquée 


ces difficultés a été exposée dans le chapitre 4. En effet l'équation (1) contient 
des variables qui seront nurmalernent 1 (1), Im, et Y, et une variable qui ne peut 
pas être I (1) à savoir Pr. La logique économique veut que !es autorités 
monétaires veillent à maintenir le niveau de compétitivité des produits 
nationaux en ajustant le taux de change lorsque les prix intérieurs ont évolué 
durablement dé munière différente des prix des principaux partenaires. Par 
conséquent Pr, doit être normalement stationnaire. 

La relation (1) contient des variables d'ordre d'intégration différents et ne peut 
constituer une relation de cuïntégration, c'ést-à-dire une relation de long 
terme. 


Par contre 


(2) logim, = logX +4: peul être une telle rélalion si on suppose que les 
autres vatiahies n'ont pas d'influence en longue période. 


Dans fa mesure où &y > 1 la croissance des importations déplace des 
productions intérieures. En effet la production intérieure ost calculée par solde 
dans l'équation d'équilibre ëntre ressources cl emploi : 
(DNim,+k=C,+1+G.+X, 

CL, X, et Im, sont déterminés par des relations de comportement et G, est 
cxogtne. 

On peut cependant juger que ce mécanisme est insuffisamment contraignant, 
dans la mesure où, en particulier, la croissance des OHEUON prend appui 
sur la croissance du revenu ou de l'activité. 


On est alors amené à remplacer la fonction d'importation par une fonction de 
partage du marché intérieur. Un des premiers modèles à contenir ce type de 
spécification est le modèle FiFi, construit en France à la fin des années 60 et 
utilisé pour ta préparation du 6° Plan. Ce modèlé inspiré de la théorie des 
économies concurrencées, élle même inspirée des travaux menés én Norvège 
par Aukrust, présente une version extrémiste du mécanisme de partage. 


Dans le court terme les prix sont Üxés par lc marché mondial. La France est un 
pays non dominant, dans une situation de price taker. 
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Au niveau de prix P,, les producteurs nationaux offrent une quantité Y. La 
différence entre la demande intérieure Di et Ÿ, est comblée par les 
importations (1), 

À moyen terme, les importateurs relèvent leurs marges lorsque leur part de 
marché augmente, cé qui fait que le niveau de prix n'est plus horizontal et les 
entreprises se modernisent et réduisent leut coûts, ce qui déplace vers le bas la 
fonction d'offre de O1 en O3. Les importations s'établissent à DT’ - Y'. 


(D M, Aglieta, R. Courbis, C. Seibel “Le modèle FiFi” Collection de l'INSEE C22, 1973, 
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Cette approche consacre là perte quasi totale du contrôle de l'activité 
économique par les moyens traditionne]s de la politique économique qui 
agissent sur la demande globale. En effer un déplacement vers la droite de La 
courbe de demande se répercute presque intégralement sur Les importations et 
laisse inchangée ja production intérieure et l'emploi. Face uux traditionnels 
multiplicateurs de demande, il existe donc des multiplicateurs d'offre qui 
mesurent l'effet du desserrement de la contrainte de compétitivité sur la 
production et sur l'emploi. 
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On peut regretter que cette analyse n'ait pus reçu les prolongements qu'elle 
mérite {1}, Quoi qu'il en soit de nombreux aultes travaux Ont &é faits avec des 
mécanismes de partage du marché dont Ja forme générale est : 


(4) e = a(LPr,+ f{#)+dTuc, +u, 


Le 


(5) Le = at L)log Pr, + fi) duc, 0, 


? 

a(L) et œil} étant des polynômes de retard ct t désignant le temps. 

Lc partage du marché dépend de Ia compétitivité prix des importations avec une 
distabution de retards de T'ardre de 6 rrimestres, d'un terme tendançiel et d'un 
indicateur conjoncture] qui peut être représenté par le taux d'utilisation des 
capacités. 

On remarquera que ces équations ne satisfont pas les conditions posées au 
chapitre à. 


Dans certains pays, il existe des restrictions aux importations sous la forme de 
licences uù d'allocations de devises ou de limitations quantitatives, Les 
importations réalisées sont alors inférieures à la demande d'importation. 1] faut 
dans ce cas distinguer les importations contraintes et les importatious 
potentielles. 


Cette distinction est fondamentale lorsque l'on envisage d'assouplir le régime 
des importations, Sc rélérer aux valeurs des élasticités revenu et prix de la 
fonction d'inportätions réalisées conduit à sous-estimer gravememt Les effets 
d'une éventuelle libéralisation. Si on connaît la fonction d'utilité des autorités, 
on peut en déduire le mécanisme de rationnement qu'elles appliquent entre les 
importations demandées ct les importations réalisées. En sens inverse, en 


() Bien qu'elle ait incontestablement un caractère dérangeant par rapport au mouvement 
de mondialisation de l'économie, Mais peut-il y avoir un équivalent du “politically correct” 
én économie ? 
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connaissant le mécanisme de rationnement et les importations réalisées, 11 est 
possible, dans certains ças, de remonter à la fonction de demande 
d'importations. 


Différents modèles théoriques permettent de faire cette opération qui n'est pas 
sans rappeler l'identificalion présentée au chapitre 2. Nous présentons ci- 
dessous un modèle dû à Goldstein et Khan {1}. 


La fonction de demande d'impartations émane du secteur privé et s'écrit : 
{6) mi = Xp + di} + da Pre 


Cette fonction n'est pas entièrement satisfaite faute de moyens de paiement 
suffisants. Les autorités offrent un montant de dévises qui dépend de leurs 
propres objectifs. Ces objectifs sont à la fois des objectifs de court terme et 
des objectifs dé long terme : 


- à court terme, les autorités cherchent à satisfaire le muximom de demandes 
d'importations et à éviter de trop brusques variations du niveau d'importations. 


- à long térme, elles cherchent à maintenir Îes importations en ligne avec (a 
capacité de paiement, f,. 

Elles désirent également maintenir un certain niveau de réserves en devises, r° 
Les réserves en t sont données par : 

(D n=f-mt+r 

f,= flux de recettes provenant dés exportations et des entrées de capitaux. 

m,= importations réalisées en t. 

Les autorités cherchent à minimiser le coût de déviation par rapport à une 
fonction de coût quadratique, 


min C,= BCm, mé} + 80m, -m,) + Bin -rf) +Batm, mi) (8) 


(D) M. Goldstein et M.S. Khan “The Supply and Demand for Exports : a Simultançous 
Approach” The Review of Economics and Statistics, 1978, 
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où m, représente Les importations compatibles avec la capacité dc paiement à 
long terme : C7 =}. 

Le coefficient B; représente Le coût associé à un écart par rapport à l'objectif 
correspondant. 

La valeur optimale des importations est celle qui minitnise le coûLl total de 
déviation mesuré par la fonction €.. Elle s'obtient en annulant la dérivée par 
rappott à Me: 

FC, 


Fr sachant que = f—m,+r 1 (2) 


ue = 2h Lin —2 Bim£ +2B;m, — 2A5m,. + 2Bzm, _ 2B:f, +2 Br 
{ 

—2Bar,_| +2Bam,-2Bym =0 (9) 

Soit 


(10) (Bi + Bo + Ba + Bad, = Ram + Bin, + Baf, — Bar + Bars + Ban =0 
En norant les coefficients, il vient : 


5. 
bi et Eh oil 
CR 


£ 


et en remplaçant dans (10) : 

QU) om = him + bams_n + ba + batré ri 1)+ ame 

2} mu sb +ba)f +batr =) +bm,_ +hym 

Soit en remplaçant md par sa valcur dans {6) 

my = Cha ba fs + bar 6) + homes + bei + bia Pr, + bic 
Ce qui conduit à estimer : 

(13) om =ap+eh + a Pr+asf tar nr 1 )+ sm +, 


avec : 
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b3+h1=0 
b;=as 
b}=as 
h=1-a;-as 
& a] 
M=—= 
h l-u3-0s 
ä a 
a =? e 


Si les variables figurant dans l'équation (13) sont les logarithmes des 
variables, les cocfficicnis a et a3 représentent les élasticités apparentes des 
importations par rapport aux revenus et aux prix, tandis Que &, et &; sont les 
élasticités réclics de la demande d'importation par rapport aux mêmes 
variables, Même si les autorités accordent un poids prépondérant au premier 
objectif ct s'efforcent de satisfaire un maximum d'impuriations, te résultat 
montré {a possibilité que les “vraies” élasticités soient substantiellement 
supéricures aux élaslicités apparentes qui résulteraient d'une régression directe 
sur les volumes d'importations réalisés. 


Comme indiqué dans le chapitre Ÿ, les fonctions d'importations peuvent 
Également Glre estimées en utilisant une spécitication à correction d'erreur. 

La formulation la plus générale est : 

(14) AlogiM, =k+ foñlog Y + ÉpAlogPr - yilog FM. log -logPr, ) 
qui peul être estimée par ta méthode de Engle et Granger en 2 étapes (1} sous 
celte forme où plus simplement par les MCO 3 partir de l'équation suivante : 
(S)alogiM,=k+e Alout +43AlogPr, +aslog M, _] + @1 log +as log Pr, +4, 
avec les propriétés suivanes : 


LE) Qn sait que dans ce cas. la relation de cointégration entre log [M, log Y,ut log Pr n'est 
pas unique écf, ch. 4}. 
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Élasticités à court terme des importations : 
— Par fapport Au FEVENU : Fo = 
— par rapport aux prix : Éo = 
Élasticités à long terme des importations : 


as 
— par rapport au revenu:  Ê=—< 
a 
ï —âs 
— par rapport aux pfix : E=— 
43 


SECTION 2 : LES EXPORTATIONS 


Les exportalions d'un pays sont par définition les importations du Reste du 
Munde. Toutefois Le Reste du Monde n'est pas, dans la réalilé, un ensemble 
homogène ct facilement identifiable comme on le suppose dans les modèles 
simples à deux pays, mais un ensombic hétérogène de pays à Fa fois clients ct 
concurrents. 


De là il vient que si la spécification d'ensemble des fonctions d'exportaions 
est similaire de celle des fonctions d'importations, la signification des 
variables qui y figurent es! différente. 


La variabie d'échelle, correspondant an revenu ou à l'activité dans la fonction 
d'importation, est donnée ici par le revenu où les importations des pays 
clients. Les imporrations présentent l'avantage d'être plus près de la démande 
adressée par les pays clients au Reste du Monde. Mais parmi les pays clients, 
tous ne pèsent pas du même poids et par exemple Les conséquences d'une 
augmentation du Revenu qu du Produit Intérieur Brut de l'Allemagne sur les 
exportations de la France ne sont pas les mêmes que celles d'une augmentation 
identique du Revenu où du Produit Intérieur Brut de ta Malaisie. C'est pourquoi 
on construit l'indicateur de demande étrangère en pondérant les Révenus ou [és 
PIB ou les Importations des pays clients exprimés en unité commune, par leur 
part dans les exportations du pays étudié. 


Ainsi pour la France, l'INSEE construit deux Indicateurs : 
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PIBET = Z PIB des pays clients de la France pondérés par leur part dans les 
exportations françaises. 


IMPET = ZX Importations des pays clients pondétées par leur part dans les 
exportations françaises. 


Ce dernier indicateur est voisin de ce qu'on désigne parfois en économie 
intemationale sous le nom de ‘emande adressée”. 


La compétitivité - prix des exportations peut être mesurée par le rapport des 
prix à l'exportation sur les prix des produits concurrents, le tout exprimé en 
une unité monétaire commune, Les différents pays concurrents ne pèsent pas le 
même poids dans le commerce mondial ét lés conséquences d'une variation de 
la compétitivité prix des exportations françaises par rapport aux prix italiena 
ne sont pas les mêmes que la même variation par rapport aux prix argentins. Il 
faut donc pondérer les indices de prix des pays concurrents par leurs parts de 
marché respectives. On peut sur cette base construire un premier indice de 
compétitivité - prix. Mais les exportations d'un pays ne se réalisent pas sur ut 
vaste marché homogène qui serait constitué par le Reste du Monde. Elles ont 
lieu sur différénts marchés nationaux qui représentent chacun une certaine 
partie du total ét sur ces différents marchés, la concurrence des autres pays est 
plus ou moins marquée. 


Un bon indicateur de compélitivité-prix doit établir pour chaque pays client un 
indice du prix des exportations par rapport aux prix du pays client et à ceux des 
principaux concurrents pondérés par leur part du marché dans le pays client. 
Puis la pondération de ces indices élémentaires par Ja part des différents 
marchés dans les exportations lolales permet d'obtenir l'indicateur de 
compétilivité-prix globale. 

Suppusons que l'on ait identifié m marchés clients désignés par i = 1. m dant 
la part dans les exportations est e; à e,, el qu'il existe n pays concurrents 
désignés par ÿ = 1 .…. n dont les parts de marchés sont respectivement e; Y. 


&Y= part du pays y dans le marché i. 
Ze;y=1 et Le; =1 
à f 
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L'indicateur de compétitivité à l'exportation sur Le marché i s'écrit : 
comp; = SRE 
TP) 
Y 
P étant l'indice des prix à l'exportations en unité monétaire de référence. 
P;7 : l'indice de prix, en unité monétaire de référence, du pays y sur le marché i. 
L'indicateur global est alors : 
Com p = T{com pr, )* 
t 
Ces variables peuvent être incluses dans une formule de Luxembourg 
symétrique de celle posée pour les importations : 
x, = get, Comp,, t} 
où Det désigne là dernande étrangère. La variable Lemps est souvent omise, 
comme dans la suite de cette section. 
Cette formulation peut être spécifiée de façon linéaire : 
(t) X,= her +h, Comp, + b3 +u, 
dans ce cas h1 représente la part du pays dans un accroissement du commerce 
international. 
ou Tag linéaire : 
{2} log X,= Blog Der + B; logComp, + PB +u, 
8, est alors l'élasticité des exportations par rapport à la demande ou par 
rappott au commerce mondial et 8 est l'élasticité - prix. 
B, se compare à la valeur de rétérence qui est l'unité : 


B,=1 les exportations du pays augmentent au même rythme que la demande 
étrangère. La part du pays sur ses principaux marchés est stable. 


B,>1 Les exportations augmentent plus vite que les principaux marchés, le 
pays gagne dés parts de marché. 


B<1 Les exportations augmentent moins vite que {s demande étrangère, le 
pays perd des parts de marché, 
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Les valeurs de B; sont cn général légèrement supérieures à 1 pour la plupart des 
pays développés à l'exception du Royaume-Uni où elle est nettement 
inféricurc, Ces valeurs traduisent la qualité des spécialisations et la capacité 
des économies à offrir &s produits qui font l'objet d'une forte demande sur le 
marché mondial. 


La comparaison des valeurs de long terme des élasticités revenu des 
importations, &, et des exportations par rapport à à demande étrangère, B;, 
donne des indications sur la capacité d‘une économie à soutenir un laux de 
croissance égal ou supérieur à celui du Reste du Monde. 


Un pays pour lequel &, < B, dispase de marges de croissance, qui toutes choses 
égales par aiïlieurs lui pérmettent de connaître un taux de croissance 
durablement supérieur à celui du reste du mande sans mettre en péril son 
équilibre extérieur. 

Le Japon a pendant longtemps bénéficié d'une telle situation. 


Âu contraire, un pays pour Lequel &, êst nettement supérieur à &; viendra huter 
sur là contrainte extérieure dès que sa craissance interne aura tendance à 
s'accélérer, Un lei pays ne peut suivre le rythme de croissance du reste du 
monde, qu'en compensant par des gains de compétitivité-prix obtenus par des 
dépréciations de sa monnaie. 


En fait, la symétrie qui existe entre Ics équations d'importations et 
d'exportations est plus formelle que réelle. D'une part, il existe des délais de 
réaction beaucoup plus longs dans Le cas des exportations et d'autre part ia 
distinction entre exportations demandées et réalisées s'impose avec davantage 
d'évidence. 


La réaction des cxportations cst beaucoup plus lente que celle des 
importations. Ce qui impose presque toujours l'utilisation de retards en 
données trimestrielles et parfois en donnécs annuelles. Ces retards 
s'appliquent à la demande et aux prix. Une solution simple consiste à 
introduire la variable endogène décalée pour traduire la présence d'un 
mécanisme de retard à poids géumétriquement décroissants : 
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(33 log X,= fhilog Dei+ P, log Comp, + Blog X,_1 + Ba +u, 
ce qui permet de faire apparaître des élasticités de court terme : 
CT : élasticité à (a demande B; 
élasticité aux prix B> 
et des élasticités de long terme : 


LT : élasticité par rapport à la demande : À 
T3 
Free #3 
élasticité prix : 
1-8 


Exemple : 


L'équation (3) à été utilisée pour estimer [es exportations de produits 
manufacturés du Marac vers la CEE. 


On a vérifié dans un premier temps que toutes les variables sont intégrées 
d'ordre 1 y compris l'indicateur de compétitivité caléulé par rapport à un 
échantillon de pays concurrents. 


L'ajustement réalisé donne sur 1985 - 1] - 1992 - 4 


LX, = 0,36 LIMPET + O,12 I. Pr, + 0,30 Abus L 47 set Aa ji 1,25 
(204) {0,47} (D.% (0,6} 


k? =0,68 DW=I,89 
{toutes les variables sont en log} 


ou LX, = indice de volume des exportations désaisonnalisées de produits 
manufacturés vers La CEE. 


LIMPET : indice désaisonnalisées des importations des pays dé la CEE 
pondérées par leur part dans les exportations marocaines. 


LPr, : indice de prix relatif du Maroc / CEE 


Lcomp, : indice de compétitivité du Maroc par rapport à un échantillon de pays 
concurrents. 


Les éfaslicités prix ne sont pas significatives, 
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L'élasticité à la demande CEE est de 0,36 à court terme et de 0,705 à long 
térme, 


On peut également utiliser des retards €chelonnés sur Det et des retards 
polynomiaux sur Comp. 


(4) log X, = filog Det + fl log De,_] + flog De,_5 + diL)Comp, + Ba +u, 


Les exportations qui sont déterminées par l'équation (2) sont les exportations 
demandées. Elles peuvent différer des exportalions réalisées. Dans le cas des 
importations, l'offre venait du Reste du Monde et pouvait être supposée 
infiniment élastique. Dans le cas des exportations, l'offre provient d'un seut 
pays et cette hypothèse ne tient plus. 


IL faut donc distinguer l'offre d'exportations et la demande d'exportations. En 
tour lagique ceci impose d'écrire un modèle à plusieurs équations dans lequel 
les exportations réalisées soit égalisent l'offre et la demande d'exportation, 
soil se fixent au minirmum de l'offre et de la dernande. 

La fonction d'offre dépend de la production intérieure, de la rentabilité des 
exportations et des tensions sur les capacités de productions. 


{5) log x = À log Det + B; log Comp, + fi +u, 
(G)_X} = h(y, Tue) 


(D) X,=xf=x 

ou (7 bis)X, =min(X#, X,) 

Une solution simple cansiste à introduire dans l'équation (2) des variables 
additionnelles qui traduisent d'éventuelles tensions sur l'offre, L'équation (2} 
devient alors une équation décrivant les exportations réalisées. Cette 
transformation est adaptée pour rendre compte d'écarts temporaires résultants 
d'incapacités poncruclies, mais conduirait à une sous-estimation de l'élasticité 
des exportations pat rapport à la demande étrangère si les exportations 
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réalisées étaient de façon permanente limitées par la capacité d'offre. Le 
recours à une forme structurelle identifiable serait indispensable. 


Comme cela a été mentionné dans Ja section !, les équations d'exportations 
posent un problème d'ordre d'intégration des variables, Les variables telles que 
les exportations au la demande étrangère seront intégrées d'ordre 1 tandis que 
la compétitivité doit logiquement être stationnaire. Les mêmes raisons 
conduisent en effet les autorités monétaires à lisser Les écarts importants qui 
poucraient apparaître au niveau de la compétitivité prix, De cc fail les 
élasticités prix sant de toute Façon délicates à estimer. L'économétrie mesure 
d'autant mieux que les variations sont plus marquées, Toutefois des travaux 
récents tendent à infirmer cette hypothèse de stationnarité de la compétitivité 
et concluent à une intégration d'ordre 1 pour cette variable dans Le cas de 
l'Allemagne, des États-Unis ou de la France. 


Les modèles à correction d'erreur peuvent être utilisés pour les fonctions 
d'exportations, soit que l'un se réfère aux propriétés du théorème de 
représentation de Granger lorsque toutes les variables sont I {1}, soit que l'on 
introduise une variable stationnaire dans une spécification dont toutes les 
variables sont stationnaires : 


(8) Alog X,=k+ fyAlog Det - y(logX,_1 —-Ë log De,_))+u, 


dans ia mesure où X,et De, sunt coïntégrés (1,1) et si Ÿ est le paramètre du 
vecteur de coïntégration, il est légitime de rajouter dans cette régression un 
(erme également stationnaire : 

(9) AlogX, =k+ BoAlog Det + Blog Comp, - y(log X,_: -ÉlogDe,_1)+4, 
Le terme entre parenthèses n'est autre que l'écart en 1-1 par rapport à la relation 
de long terme. D'où une autre formulation de La fonction d'exportations : 

(10} AlogX, =k&+ foÂlog Det + 8: log Comp, -YECMr-1 +0, 


qui peut être estiméc par les MCO sous cite forme. 


Tables statistiques 


Loi normale centrée réduite, Z 
(oi normale de moyenne 0 et de 4 é 
variance 1) 
Zœ 
La table donne {a valeur de zg tétle que © = MZ>z2a} 
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Table de [a loi de khi-2 (x2) 


La table donne la valeur # 2 telle ques Pix2>y 5 
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Table de distribution de F EL de Fishce Snedecur) 
Valeurs de £'ayant ln protabstité d'être dépassées (F =? f3À) 


#p 


=] x, =2 
vs F=AUS : P=0or PU P=Q 
[l lér 4 . +052 | 41:00 
2 1851 : 200 
3 16,13 30,61 
4 LA J#4x 
L fl 1427 
| ñ EX] 10:91 
! 7 4,59 9355 
ï 8 su h4£ 
| { S,L2 #97 
Li 4% 756 
il +,44 J5$ 14 770 
[E 4,75 4,53 ZH4 6.92 
13 EAU #.70 4 
art tSL H 
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343 #24 24 
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3.53 &ûL de 
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335 549 6 LÉ 
134 s4S } 
139 342 373 
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si 64 la phos grandir les deux varfances climées avec v degrés de libertés. 


Tabke de distribution de F {Loi de Fisher.Snadecuri 
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Tables du d de Durbin-Watson 


Risque à 5 % 

T [ken [ Kk=9 k=3 | K=4 [|  «-5 

[di  d2 [di 42 dr, _d2 [#1 42 | a ar | 
15 108 1,36 [0,95 1,54 |0.82 1,75 | 0,69 1.97 |O0.£8 2,21 
16 [1,18 1.37 | 0.98 1.54 |ü,g&s 1,73 | 0,74 1,93 | 0,62 2,14 
17 [143 1.381102 1.54 |0.90 1.71 [078 1,90 |0.67 2.16 | 
18 |1,1# 1,39 11,05 1,52 |0,94 1.69 | 0,82 1.87 |OQ.7L 2.06 
19 1,8 1,40 ! 1.408 1,53 |0,97 1,68 | 0,86 1,45 | 0,75 2.02 
20 |1,20 L,41 1,10 1.54 | 1,00 1,68 | 0,90  i,k2 | 0,79 1,99 
21 1,22 1,4% 1,34 1,54 | 1,07 1,67 : 0,94 1,81 | O8 1,96 
22 |1.24 1.43 : 115 1,54 | 1.05 1,66 , 0,Jé 1,80 | 0,%6 1,94 ! 
23 1,26 L44 ji 1.17 1,54 | 1,08 1,66 0,99  i,79 | 0,90 1.92 | 
24 |E27 1,45 1119 1,55 | 1,10 1,66 1,01 1,78 | 0,92 1,90 : 
25 |1429 1,45 11.21 155 | 1412 166 104 1,77 | 0,95 1,89 
26 | 4130 1.46 122 LSS | 1,14 1,65 1.06 1,76 | 0,98 1.88 
27 | 132 1,47 11.34 1,56 | 1,18 1,65 1,08 1,76 | 1,01 1,86 
28 |1,24 1,48 | 1.28 1,56 | 1,18 1,65 1.10 1.75 | 1,03 1,85 
29 | 1,34 148 | 1.27 1.56 |1,20 1,69 1.12 1,74 | 1.05 1.84 
39 |1,35 1,49 [1,28 1,47 | 1.21 165 [1,14 1,74 | 107 1.82 
ñ1 |1,36 1.50 | 1,30 1.57 |1,23 1,65 | 1,16 1,74 | 1,09 1,83 
32 |1.37 1.50 | 1,31 1.57 |i,24 1,65 | 1,18 1,73 | 111 1.8? 
33 1.38 1,51 1,342 1,58 1,26 1,65 1,19 1,73 1,13 [ee:1 
14 1,39 1.5! 1,34 1.54 1,27 1,6$ 1,21 1,73 1,15 1,#i 
35 |1,40 1.52 | 1.54 1,58 | 1,24 1.65 | 1,22 1,73 | 1.16 1.80 
36 1,41 1,52 | 1.35 1,39 1,29 L65 1,24 1,73 t,18 t,80 
37 |142 1,43 ! 1,34 1,59 | 1.31 1.66 | 1,25 1,72 | 1.19  E.80 : 
38 |1,41 1,54 L,27 1,29 | 1,32 66 ;, 1,26 1,72 | 1,21 1,79 
39 |L41 154 138 1.60 | 1,52 1,66 | 1,27 1,72 | 1,22 1,79 
49 |Hd4 1,54 1.39 1.60 | 1,34 166 : 1,29 1,72 | 1,23 1.79 
45 148 [ST 1.47 [62 |1,38 1,67 1,34 1,72 | 1,29 1,78 i 
sQ ÀE,50 1,49 1.46 1.63 | 1,42 1,67 1,38 1,72 | 1,è4 1,77 : 
85 a,53 160 11.49 1,64 | 1,4% 1,68 | 1,41 1,72 | 1,38 1.77 : 
69 |1,55 162 | 1.51 165 | 1.48 1,69 | 1,44 3,73 | 14h 1,77 : 
65 |1,57 1.63 | 1,54 Léé | 1.50 1,70 | 1,47 1,73 | 1,44 1,77 
70 1,58 1,64 1,55 1,47 1,5% 1,70 1,49 1,74 1,46 1,77 
45 1,60 1.65 | 1.57 168 | 1.54 1,71 | 1,51 1,74 | 1,49 1,77 
Ro |1,61 1,66 . 1,59 1,69 | 1,56 1,72 1,53 1,#4 | 1,51 1,77 ! 
85 |1,62 1,67 1.60 1,70 |1,57 1,72 1,5% 1,75 | 1,52 1.77 : 
9Q ÀE.61 1.68 1.61 1.70 | 1,59 1,73 1,87 1,75 | 1.54 E.78 
95 |1.64 1.69 1.62 1,71 |1,60 1,73 1,58 1,7S | 1.56 1,78 
100 | 1,65 69 L61 1,72 | 1.61 1,74 1,89 1,74 | 15% 174 


kest le nombre de variables excgènes (constante exclue}. 
T est la taille de l'échantillon 
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